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Mercury Photosensitized Reactions of Propane* 
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(Received June 13, 1940) 


An investigation has been made of the mercury photosensitized reactions of propane, alone 
and in the presence of hydrogen at temperatures from 25 to 323°C. The products of the reac- 
tion are predominantly hydrogen, hexanes (mainly 2,3-di-Me-butane and 2-Me-pentane), and 
a small amount of higher paraffins. The postulated mechanism is: 


Hg (@P;) +C3Hs—Hg (59) +C;H;+H 
Hg(®P:) +H2—Hg('So) +2H 
H +C:H s—>C;H; +H, 
2C;H7>CeHu 
2H—H, 
C;H;+C3Hs>CeHut+H. 


The quantum yield of hydrogen is approximately 0.15 at 25°C and rises somewhat with in- 


creasing temperature. 





INTRODUCTION 


UCH valuable information concerning the 

kinetics of elementary processes has been 
been obtained from the investigation of photo- 
sensitized reactions. In the case of the mercury 
photosensitized decomposition: of ethane!~* there 
seems to be no doubt that the primary process 
involves a C—H bond split, 


C.He+Hg (Pi) =C2H;+H+Hg(' So). (1) 


It also seems to be well established that the main 





_ *A preliminary account of some of this work was given 
Na paper presented at the Symposium on Photochemistry 
of the New York Academy of Science, May 3, 1940. 

+ Present address, National Research Laboratories, 
Ottawa, Canada. 

‘E. W. R. Steacie and N. W. F. Phillips, J. Chem. 
Phys. 6, 179 (1938). 

*E. W. R. Steacie and N. W. F. Phillips, Can. J. Re- 
search B16, 303 (1938). 

*E. W. R. Steacie, W. A. Alexander, and N. W. F. 
Phillips, Can. J. Research B16, 314 (1938). 
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secondary reactions are*~* 


H+C:He=C2H;+H,, (2) 
H+C.H;=2CHs, (3) 
CH3+C.H;=CsHs, (4) 
2C.H;=C.Hw, (5) 
2H=H,, (6) 
H2+Hg(*P1) =2H+Hg('So). (7) 


In view of the interest attaching to the reactions 
of the aliphatic free radicals, it seemed desirable 
to investigate the reactions of propane as well. 

An investigation of the reaction of hydrogen 
atoms, produced by the Wood-Bonhoeffer meth- 
od, with propane has already been reported.® 


4H. S. Taylor, J. Phys. Chem. 42, 763 (1938). 


5 E. W. R. Steacie and N. A. D. Parlee, Trans. Faraday 
Soc. 35, 854 (1939); Can. J. Research B17, 371 (1939). 


E. W. R. 


Fic. 1. Apparatus. 


It was found that the products of the reaction 
were solely methane at low temperatures, and 
methane, ethylene, and ethane at higher tem- 
peratures. It was concluded that the results 
furnished very strong confirmation of the im- 
portance of reaction (3) in free radical processes. 
The main steps in the postulated mechanism 
were: 

Primary reaction 


H+C3Hs=C3H;+H:. (8) 

Secondary reactions at low temperatures 
H+Cs3H;=CHs+C:H;, (9) 
H+C;H;=2CHs, (3) 
CH;+H=CH,. (10) 
Additional secondary reactions at high temperatures 
C3H; =C2Hi+CHs, (11) 
2CHs=CsH., (12) 
H+C,H4=CoH;, (13) 
H+C.H;=C2He, (14) 
H2+C2H; =C2He+H. (15) 


As far as we are aware, no previous work has 
been reported on the mercury photosensitized 
reactions of propane. 


EXPERIMENTAL 


A few preliminary experiments were first car- 
ried out in a “single pass’’ flow system of a type 


described elsewhere. Such a system reduces 


SE. W. R. Steacie and R. L. Cunningham, J. Chem. 
Phys. in press. 


STEACIE AND D. J. 


DEWAR 


complications due to secondary reactions to a 
minimum, but is very tedious in operation. The 
major part of the work was therefore carried out 
in a circulating system as described below. The 
two types of system yielded substantially the 
same results. 

The apparatus employed is shown diagram- 
matically in Fig. 1. The mercury saturator S sup- 
plied the necessary mercury vapor to saturate 
the reactant gas. It was heated electrically to 
about 60°C. A desaturating trap D, filled with 
iron pellets, and maintained at 25°C, removed 
excess mercury vapor and maintained the mer- 
cury pressure at that corresponding to its vapor 
pressure at 25°C. The circulating pump P has 
been previously described.' It was connected to 
a series of mercury valves V, so arranged that 
the gas in the line always flowed in the one direc- 
tion in spite of the reciprocating motion of the 
pump. 

The condensable reaction products were col- 
lected in the trap T which was filled with iron 
pellets and maintained at — 80°C. 

The fused quartz reaction vessel C was annular 
in form, 10 cm long, 8 cm in diameter, with a 4-cm 
diameter hole through the middle. Its volume 
was 320 cc. The reaction vessel was illuminated 
by a Hanovia Sc2537 high voltage mercury lamp. 
The lamp was in the form of a double-U and was 
run at 100 milliamperes from a 5000-volt sign 
transformer. 

As the Hanovia lamp was unsuitable for high 
temperature investigations, the runs at higher 
temperatures were made with a combined lamp 
reaction vessel of the type previously described.’ 

The reaction vessel was provided with a 
thermocouple well, and temperatures were 
measured with a Chromel-Alumel thermocouple. 
For high temperature runs the reaction vessel 
was enclosed in an electrically heated furnace, 
with the exception of the electrodes which pro- 
jected from the furnace and were water or air- 
cooled. 

The main proportion of the apparatus used for 
low temperature runs had a volume of 1680 cc. 
The high temperature apparatus had a volume of 
2220 cc. The usual pumping system, gas reser- 
voirs, Toepler pump, gas-holders, etc., were used. 


7E. W. R. Steacie and N. W. F. Phillips, Can. J. Re- 
search B16, 219 (1938). 
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The intensity of the radiation emitted by the 
lamps was determined by measuring the rate of 
hydrolysis of monochloracetic acid. The acid 
used was 0.5 N, and the procedure followed that 
of Steacie and Phillips.' Recent work® has indi- 
cated lower values of the quantum yield than 
those originally given by Rudberg. However, 
there are some indications from work in this 
laboratory that the quantum yield may be higher 
for high light intensities. In any case the absolute 
value of the quantum yield is not of great im- 
portance for the present purpose, and Rudberg’s 
values have been used for the present to facilitate 
comparison with the ethane decomposition. 

All runs were carried out at the vapor pressure 
of propane at —80°C (about 14 cm). Excess 
propane was admitted to the system, the trap 
was cooled to —80°C and part of the propane 
condensed in it. This acted as a reservoir to keep 
the partial pressure of propane constant during 
the experiment. 

At this pressure the quenching by propane is 
sufficient to ensure practically 100 percent ab- 
sorption of the resonance radiation. 


ANALYSES 

At the end of a run the gaseous products were 
pumped by a Toepler pump into a portable 
mercury gas-holder. The gas was then analyzed 
in a low temperature fractional distillation ap- 
paratus of the Podbielniak type. The hydrogen 
and methane were taken off together in the distil- 
lation and this fraction was analyzed by com- 
bustion as previously described.” Occasional tests 
on fractions for unsaturated hydrocarbons were 
made in a Burrell gas analysis apparatus. 

The condensate in the trap, which consisted 
of higher boiling hydrocarbons, was weighed 
after each run and saved for future analysis. The 
liquid products from several runs were combined 
and separated by distillation in a semi-micro- 
fractionating apparatus." Cuts from the micro- 
distillation apparatus were analyzed by combus- 
tion and tested for unsaturates by bromine 
absorption. Determinations were also made of 
their melting points and refractive indices. 

*E. Rudberg, Zeits. f. Physik 24, 247 (1924). 
*W. G. Leighton, R. N. Smith, and P. A. Leighton, J. 
Am. Chem. Soc. 61, 2299 (1939); 60, 2566 (1938). 


'°E. W. R. Steacie, D. J. Dewar, and R. Potvin, forth- 
coming paper. 


PROPANE 


MATERIALS 


Propane was obtained in cylinders from the 
Ohio Chemical and Manufacturing Co. It was 
stated to be 99.8 percent pure and no appreciable 
impurities were detected in Podbielniak analyses. 
The gas was therefore used without further 
purification. 

Hydrogen was obtained from commercial 
cylinders. It was purified by passage over platin- 
ized asbestos at 600°C and through a liquid-air 
trap. 


RESULTS 


Some 20 preliminary runs were made in the 
““single-pass’’ system. These indicated a quantum 
efficiency of hydrogen production of approxi- 
mately 0.18. The gaseous products consisted of 
hydrogen and a trace of methane. The liquid 
products are reported in Table IV. 

The results of experiments with the circulating 
system are given in Tables I-IV. The initial and 
final pressures given in columns 3 and 4 of 
Table I are those with the trap at room tem- 
perature. The actual partial pressure of propane 
during a run was the vapor pressure of propane at 
— 80°C, viz. approximately 14 cm. 

The liquid condensate is expressed in column 
9 as ‘‘CsH,4”’ since the main fraction of the distil- 
late had this composition. The condensate also 
contained some hydrocarbons of higher boiling 
point and traces of lower boiling hydrocarbons, 
so that this method of reporting the results is 
slightly in error. 

TABLE I. The mercury photosensitized decomposition 

of propane. 

Volume of system = 1680 cc. Propane vapor pressure= 14 cm. 

Volume of reaction cell= 320 cc. Trap temperature =—80°C. 

Are current =0.100 amp. Reaction temperature = 24-27°C 


Resonance radiation absorbed =4.0X 10-6 Einstein/sec. 
Circulation rate= 1900 cc/min. 
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TABLE II. The mercury photosensitized decomposition of 


propane at higher temperatures. 

Volume of system = 2215cec. Circulation rate = or, 

Volume ofreactioncell= 870cc. Trap temperature =-8 
Are current =0.090 amp. Propane vapor pressure = 

Resonance radiation absorbed =9.5 X 10-6 Einstein/sec. 


id i: 








Propvwcts, 


Int- | Decompo- Mote Percent 


TIAL SITION 
PrEs- Rate, 
SURE, | MOLES PER 
SEC. X 106 H2 


QUANTUM 
YIELD 





CHa | “CeHu” 
57.9| 0.0 42.1 
53. 46.2 
53.9! 0.0 46.1 
49.0} 0.2 50.8 
52.6} 1.9 45.5 


519} 3. 44.6 
57.1] 2. 40.6 


553 | 0: 44.7 

55. 44.4 
} 61.4 Ap 36.8 
0.33 45.2 


a 
Blank run lamp off 
1.7 43.8 














(0.24)* | 54. 
(0.26)* | 55.0] 03 
(0.27)* |51.3| 1.7 
(0.28)* | “53.7 


(0.15)* BLA 
































* In these runs at higher temperatures, difficulty was encountered in keeping the 
electrodes sufficiently cool (owing to the fact that the electrodes were too close to 
the reaction vessel). Under these circumstances the mercury vapor pressure was 
too high, and \2537 was partially reversed and therefore less intense than before. 
As a result these quantum yields have no real significance. At 175°C, however, the 
intensity was as stated, and the rise in quantum yield on going from room tempera- 
ture to 175°C is definitely established. 


The results of experiments in the presence of 
added hydrogen are given in Table III. In this 
table the quantity of hydrogen in the products is 
that formed during the reaction, i.e. the amount 
of hydrogen added initially has been subtracted 
from the amount found by analysis. 

The results of analyses, etc., of the liquid prod- 
ucts are given in Table IV. In these tests the 
products from several runs were combined in each 
case so as to provide a sample of 4 to 10 cc. In 
the table the high boiling products are expressed 
in volume percent of the total condensate. 

It will be seen that, except for a trace of un- 
saturates, the main fraction of the liquid prod- 
ucts consists of hexanes. Furthermore, within the 
accuracy with which distillations on very small 
amounts of liquid can be carried out, there is no 
appreciable change in the products with changing 
temperature, or in the presence or absence of 
added hydrogen. 

In Table V the physical properties of the 
isomeric hexanes are summarized. The boiling 
point of the product definitely rules out the 
possibility of -hexane, or of 2,2-di-Me-butane. 
The melting point data on Table IV are uncer- 
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tain on account of the ease of supercooling and 
the tendency to the formation of glasses. How- 
ever, the data are sufficiently reliable to rule 
out n-hexane and 2,2-di-Me-butane also. 

Of the remaining three compounds, 3-Me- 
pentane may also be ruled out. In the first place 
it is intrinsically unlikely, since neither it nor 
2,2-di-Me-butane can be formed directly by any 
possible combination of n-propyl or isopropyl 
radicals. Furthermore, its boiling point is rather 
TABLE III. The mercury photosensitized decomposition of 

propane in the presence of added hydrogen. 

Volume of system = 1680 ce. Reaction temperature = 24-27°C. 

Volume of reactioncell= 320cc. Trap temperature =—80°C. 

Are current =0.100amp. Propane vapor pressure = 14cm. 


Resonanceradiation absorbed =4.0X 10-6 Einstein/sec. 
Circulation rate= 1900 cc/min. 








Propvcts, 


PRESSURES Moe Percent 








INITIAL 


Decom- 
PRESSURE 


POSITION 
Rate, | QUANTUM 
Mo.es YIELD 

PER SEC. | (or He 
106 +CH,) “CeHu” 


0.64 0.16 51. mi 47.3 
0.64 0.16 54. 45.3 
0.67 0.17 52. J 47.8 
0.67 0.17 , 
0.70 0.17 = 






































high, its melting point low, and its refractive 
index high. Hence all the data combine to make 
its presence unlikely. 

We are therefore left with the possibility of 
2-Me-pentane and 2,3-di-Me-butane. 2-Me-pen- 
tane could be formed by the combination of a 
n-propyl and an isopropyl radical, while 2,3-di- 
Me-butane would result from the combination of 
2 isopropyl radicals. As far as the data of Table 
IV are concerned, the distillation data are not 
reliable enough to enable any decision to be 
made between the two isomers. The refractive 
indices slightly favor 2,3-di-Me-butane, while 
the melting points somewhat favor 2-Me-pen- 
tane. 

It may be observed that if, as appears prob- 
able, the hexanes result from the combination of 
two propyl radicals, then ignoring steric factors 
we may calculate the relative amounts of 2,3-di- 
Me-butane, 2-Me-pentane, and n-hexane which 
would result from any given relative amount of 
n-propyl and isopropyl radicals. Thus if we had 
equal amounts of the two radicals we would ex- 
pect equal quantities of the three isomers. Now 
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n-hexane is not detected. If we had say 5 iso- to 
1 n-propyl radical, we would expect 


25 2,3-di-Me-butane 
5 2-Me-pentane 
1 n-hexane. 


It is certain that this amount of n-hexane would 
have escaped detection. It may be concluded 
therefore that on the basis of chance recombina- 
tion of radicals, the absence of n-hexane indi- 
cates that the condensable products are probably 
preponderantly 2,3-di-Me-butane rather than 
2-Me-pentane. This argument is, of course, 
weakened by steric considerations. However 
it receives support from the generally accepted 
idea that secondary hydrogen atoms are more 
weakly bound than primary. 

The physical properties, combustion data, etc., 
indicate that the high boiling fraction of the 
products is probably predominantly nonanes. 
These are presumably formed by secondary re- 
actions of propyl radicals with hexane or with 
hexyl radicals. 

It may be remarked that since the ‘‘single- 
pass’”’ runs and those in the circulating system 
yield the same products, it may be concluded 
that the products formed are those actually pro- 
duced in the reaction, and that secondary reac- 
tions have little effect on the results. 


TABLE V. Physical constants of the hexanes.* 








CARBON M.P., | B.P., 
SUBSTANCE SKELETON *~< — 


C-C-C-C-C-C| —94 | 68.8 | 1.37506 





n-hexane 


Cc 


2-Me-pentane 60.2 | 1.37149 


2,3-di-Me-butane 1.37499 


. 


| 
c-C— 

Cc 

| 
a on 


>=C=-C 
Cc -C 


G 
3-Me-pentane 


2,2-di-Me-butane 1.36887 




















. From Egloff, Physical Constants of Hydrocarbons (New York, 1939), 
Jol. 1. 


DISCUSSION 


The fact that the products are almost exclus- 
ively hydrogen and hexanes, indicates that, as 
with ethane, the primary step is a C—H bond 
split, 

C3Hs+Heg(@P1) -C3H;+H+Heg('Sp) 
or —C3;H;+HgH. (16) 


Furthermore the following secondary reactions 
are obviously reasonable, 


2C3H;-C,Hu, (17) 
2H—H,g, (18) 
H+C3Hs—-C3H;+ H2. (19) 


TABLE IV. Liquid products of the reaction. 








MAIN FRACTION 


Hi1GH BOILING FRACTION 





RUNS — 


B.P. COMPOSITION 
(BY COMBUSTION) 


AMOUNT 
(AS VOLUME 
PERCENT OF 
TOTAL LIQUID B.P. 

PRODUCTS) » 


REMARKS REMARKS 





Preliminary 


single-pass 57.5+1 C5.90.1 


runs 13.70.5 


Unsaturates, ca. 10% 142-149 Cs. 90.5 
<5% H20.7+1.0 





1-7 58.0-58.5 Co.0+0.1 


—_ present — 





8-12 
32-36 } —154 58.8-59.4 |1.3735 


present 





(a) 59.2-59.7 |1.3736 
(b) 60.2-61.0 |1.3733 


13-18 — 169 


Trace of lower 
boiling fraction 10 
also present 





(a) —167 58.6-59.9 |1.3736 
(b) —169 61.0-61.2 |1.3734 


Trace of lower np? = 1.3997 
boiling fraction 


to 
also present 1.4003 





— 167 58.8-59.4 |1.3725 



































576 E. W. R. STEACIE 
Since the products are unchanged when hydrogen 
is added initially, we obviously have in this case 
the same secondary reactions as above, together 
with the primary step 


H.+ Hg(@P,)-2H+ Hg('So) 


or —HgH+H. (20) 


At first sight these results appear to conflict 
with those of Steacie and Parlee on the reaction 
of hydrogen atoms with propane. They suggested 
reaction (19) as the primary step in their investi- 
gation. They accounted for the formation of 
large amounts of methane (which are not found 
here) by reaction (9). 


H+C3;H;-C.H;+CH 3. (9) 


At higher temperatures these results also sug- 
gested strongly that the propyl radical was de- 
composing by 


C3H;-C.H.+CHs, (11) 


while here the products do not change appre- 
ciably with rising temperature. 

However, in the case of mercury photosensiti- 
zation the concentration of hydrogen atoms is 
very low, and that of propane relatively high, 
thus making (9) unimportant compared with 
(19). Furthermore, recombination reactions of 
radicals will be favored by the higher pressure, 
thus diminishing the chance of (11). There is 
also the possibility that hexane may occur here by 


C3H;+C3Hs—>CsHut+H, (21) 


on account of the relatively high concentration 
of propane. The results of the investigation of the 
reaction of hydrogen atoms with propane are thus 
not in any way incompatible with the present 
work. 


AND D. J. DEWAR 

The assumption that reaction (21) plays an 
important part appears to be unavoidable. It 
seems to be very unlikely that propyl radicals 
could disappear fast enough by recombination 
alone at high temperatures to avoid some decom- 
position of the propyl radical by (11). If this oc- 
curred, methane and ethane would ultimately be 
formed. Since they are not present in appreciable 
quantities, propyl radicals must disappear by 
some reaction such as (21). 

It may be added that the absence of n-hexane 
as an important product favors considerable 
hexane formation by (21). We may expect that 
the relatively weaker binding of secondary hy- 
drogen atoms will lead to a preponderance of iso- 
propyl radicals in reactions (16), (19), and (21). 
If this is the case, then (21) will lead almost exclu- 
sively to 2,3-di-Me-butane and 2-Me-pentane. 
The occurrence of reaction (21) would also im- 
part chain characteristics to the reaction, and 
would explain the increase in quantum yield at 
higher temperatures. 

We may therefore conclude that the main reac- 
tions of importance in the mechanism are 


Hg(®P1) +CsHs—Hg (So) +CsH7+H 
Hg(?P;) +H2—Hg(!S0) +2H 
H+C;H;—-C;H:;+ H2 
2C3H;—>CeH 4 
2H—H, 
C3H;+C3Hs-CeHu+H. 


In a forthcoming paper the present results will 
be compared with the results of an investigation 
of the cadmium photosensitized reaction. 

The award of a Studentship to one of us (D. 
J. D.) by the National Research Council of 
Canada is hereby gratefully acknowledged. 
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The Infra-Red Absorptions of Organic Hydroxyl Compounds in the Solid State 
in the 3u Region 


MANSEL M. DAvIiEs 
The Department of Physical Chemistry, Cambridge, England 


(Received May 30, 1940) 


The absorptions in the range 2.6u to 3.54, covering the OH and CH fundamentals, have been 
measured for a variety of solid organic hydroxyl compounds in the form of thin films. For 
simple alcoholic and carboxylic structures the results, as expected, show complete association 
in the condensed state. In a number of other cases a sharp band with a narrow doublet structure 
(Av~40 cm!) is found near 2.9u, and this is ascribed to relatively ‘‘free’’ OH groups in the 
solid resulting either from steric protection, as in triphenylmethyl carbinol, or weak intra- 
molecular binding as in ethyl mandelate. Two very general features are the shift of the absorp- 
tions to longer wave-lengths in the condensed phases, and the doublet nature of the absorption 
of the ‘‘free’’ hydroxyl groups in the solid as well as in other phases. Both these factors are 
discussed in the light of the new and previous data. 





HE examination of infra-red spectra has 

been extensively used in recent years to 
detect the inter- and intra-molecular interactions, 
more particularly of hydroxylic compounds. If 
only in view of the possible extensions of this 
method to such structural problems as occur in 
the carbohydrates, the proteins, and related 
synthetic polymers it was obviously of interest 
to examine first the absorptions of simple 
hydroxylic compounds in the solid state. A con- 
siderable proportion of the various hydroxylic 
derivatives studied had already been examined 
in dilute solution, usually in carbon tetrachloride. 
Comparison of the results for the two states leads 
to conclusions as to the degree of interaction of 
the hydroxyl groups in the solid and the extent 
to which “hydrogen bond” formation occurs in 
these cases. The results for the compounds 
studied are reported individually and are followed 
by a discussion of the more general features 
observed. 


EXPERIMENTAL 


As in previous work! the Hilger D. 83 spec- 
trometer has been used with the quartz prism. 
Having realuminized the mirrors it was possible 
to work with slit widths of 0.0204 from 2.7 to 
3.5u: generally 5\ was 0.025y. 

The usual form in which the solids have been 
examined is as a film pressed between two very 
thin glass slips. Care was taken to prove that 
the latter had no appreciable absorption in the 


‘ Davies and Sutherland, J. Chem. Phys. 6, 755 (1938). 


region concerned and repeated measurement 
showed that two such glass slips together gave 
80+2 percent transmission from 2.74 to 3.5y, 
the 20 percent loss of energy arising from four 
reflections at air-glass interfaces: the absorption 
curves have been corrected for this loss. A small 
quantity (some 10 or 20 mg) of the substance to 
be examined is gently heated to melting on one 
such glass slip (2.5 by 4 cm) and the second, 
previously warmed, pressed on it. This gives a 
reasonably uniform film whose mean thickness 
has been estimated from the area occupied and 
the loss in weight on cleaning the slips. Such 
films invariably scatter a certain proportion of 
the incident radiation, which accounts for the 
apparent background absorption in most of the 
curves: this scattering is minimized by inserting 
the film as nearly as possible at the focus in the 
absorption path, and in any case is considerably 
less than for similar films deposited from solu- 
tions. Liquids have been examined in the same 
way, the two glass slips holding together with 
no apparent loss during the course of a run. 
Cells of the usual type fitted with fluorite 
windows have been used for solutions. 

The substances examined, unless of sufficient 
initial purity, were first carefully recrystallized 
or fractionally distilled. 


RESULTS 


In solution the simplest type of absorption 
corresponding to the valence O—H vibration is 
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Fic. 1. Absorption of thin films of solids. (a) Cetyl 
alcohol, (b) lauric acid, (c) 1-Pentadecanol 15-acid. Wave- 
length scale (in 4) below, frequency scale (in cm) above. 
For curve (b) add twenty percent to percentage absorption 
scale on left. 


given by the primary alcohols. This is found to 
consist of two components: A very sharp well- 
defined band, showing no structure, is centered 
at 2.75 and certainly arises from the free or 
monomeric molecules, while the intensity of a 
broad band extending from at least 2.8u to 3.15, 
with a maximum near 3.0u, increases with the 
number of associated molecules in solution.? 
Cetyl alcohol, C1gH330H, typical of the straight 
chain compounds, in a film estimated to be 10# 
molecules thick gave the absorption curve shown 
in Fig. 1. The absorption which spreads from 
2.84 to beyond 3.24 corresponds almost exactly 
to the associational band in solution, with the 
difference that the maxima in the absorption are 
now given as 3.02u and 3.10u, whereas in carbon 
tetrachloride they are at 2.87u and 2.99yu. The 
CH absorption, within the limits of its definition 
in this case, appears to be in its usual position 
between 3.44 and 3.5u. This result is entirely 
what would have been expected from the known 
x-ray crystal structure for such alcohols, where 
the hydroxyl groups come together and form a 
layer lattice with spacing twice the molecular 
length.* Obviously the association of the hydroxyl 


2 Fox and Martin, Proc. Roy. Soc. A162, 419 (1937). 
3 Malkin, J. Am. Chem. Soc. 52, 3739 (1930). 
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groups has now proceeded to completion, as may 
be represented by (a). The shift of the peak 
absorption by about 0.1y to longer wave-lengths 
is found to be a frequent characteristic of OH 
bands on passing from solution, or vapor, to the 
solid state. 

In further agreement with expectations, lauric 
acid, Cy,He3-COOH, a simple, long chain, car- 
boxylic acid, shows a rather featureless absorp- 
tion starting below 2.9u and extending beyond 
3.5u, the CH bands being only just discernible 
as peaks around 3.45y (Fig. 1b). The greater 
extent of the association band for the acids has 
already been found in solution. The molecules of 
such acids, of course, crystallize in the same 
head-to-head fashion as the alcohols.* 

The typical appearance of the normal alcoholic 
and carboxylic hydroxyl absorptions in carbon 
tetrachloride solution at 2.764 and 2.83 is 
shown by the curve for 15-pentadecanol 1-acid, 
HOOC: (CH2)14- CH2OH, in Fig. 2. In the solid 
it shows the most featureless absorption, starting 
at 2.84 and continuing into the CH bands at 
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Fic. 2. 0.0189 molar 15-pentadecanol, 1-acid in 3.05-mm 
cell. Temperature 58°. 


* Brill and Meyer, Zeits. f. Krist. 67, 570 (1928). 
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3.4u (Fig. 1, curve c). This conforms to a likely 
supposition that the crystal structure consists of 
alternate double layers of COOH and OH groups: 
—BB AA B B—. Association in pairs 
of the carboxyl heads is known to be pronounced 
even in dilute solutions.® In the solid both groups 
will be completely associated and their corre- 
sponding absorptions smeared into the unrelieved 
structure found. 

Some of the most interesting results for the 
solid state are provided by triphenylmethyl 
carbinol, (CsH;)3-COH, Fig. 3. This compound 
gives in the range 2.6u to 3.5u two relatively 
sharp absorptions whose identification as arising 
from the OH and CH vibrations follows immedi- 
ately on comparison with the absorption of, say 
naphthalene, in the solid state: the hydrocarbon 
shows only the band at 3.3y. For triphenyl- 
methyl carbinol the CH band appears to have 
much the same structure as in solution, being 
correspondingly sharp but apparently shifted 
0.024 to longer wave-length. The OH absorption 
is broader than is usual in solution and now has 
a doublet structure, with peaks at 2.86;u and 
2.90 ou of equal intensity. The displacement from 
the center of the sharp peak in solution (at 2.769y 
in carbon tetrachloride) is considerable: Av = 137 
cm; the frequency difference of the compo- 
nents of the doublet is 42 cm~'. These features 
will be discussed later, but it may be suggested 
here that the shift from 2.774 to 2.884 may be 
related to some directional field in the solid 
(electrostatic) as the macro-physical character- 
istics of the hydrocarbon medium (i.e., provided 
by the C,H; groups) would appear to resemble 
those of the solvents CCl, and CHCl ;. The 
relative sharpness of the band, besides its 
position, is entirely against ascribing it to 
associated molecules, and the absence of any 
appreciable absorption characteristic of the latter 
shows that even in the solid state the steric 
influence of the bulky (CsH;)sC group prevents 
the hydroxyls of different molecules getting 
sufficiently close for their molecular fields to 
interact considerably.? This conclusion is readily 
appreciated from a molecular model for this 
compound. 

Another interesting feature in Fig. 3 is the 








* Davies, Trans. Faraday Soc. 34, 410 (1938), and J. 
Chem. Phys. 6, 770 (1938). 


obviously greater intensity of the OH band 
compared with the CH, despite the great pre- 
dominance of the latter groups. Relative to the 
CH, the OH absorption seems to increase in 
intensity about four times on passing from solu- 
tion to solid. This recalls the increase in intensity 
of the HCI fundamental absorption between the 
gas phase and solution.® The intensity variation 
may prove to be general for strong dipole groups 
on condensation.* 
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Fic. 3. Absorption of thin film of solid triphenyl methyl 
carbinol. d=0.012 mm. 


Two a-hydroxyacids have been examined both 
in carbon tetrachloride solutions and in the 
solid state. These compounds, mandelic acid 
(CsH;-CHOH-COOH) and a-hydroxymyristic 
acid (Cy2H2;-CHOH-COOH) show entirely simi- 
lar results (Fig. 4a and b). In solution the normal 
alcoholic OH absorption is missing, being repre- 
sented at most by a shoulder at 2.78u on the main 
peak which is perfectly sharp on the long-wave- 
length side, and is centered at 2.83y. This result 
is not so unexpected as at first appears. Thus, in 
the a-hydroxyketone, benzoin (CsH;-CHOH-- 
CO-C,H;) the main OH absorption appears at 
2.88 although here the ‘‘normal” band at 2.78u 
is clearly resolved from the latter.? Apparently 
the intramolecular interaction is not so intense 
in the a-hydroxyacids. From the form of the 
curve it would appear that there are only two 


6 West and Edwards, J. Chem. Phys. 5, 14 (1937). 

* An attempt was made to examine the triphenylmethyl 
carbinol absorptions under the high dispersion of a grating 
instrument. With the present arrangement of the instru- 
ment, the higher orders overlapping the 3.0 first-order 
spectrum are eliminated by the filter-shutter method 
(Sutherland, Lee, and Conn, to appear shortly), but this 
proved unsatisfactory for the present purpose owing to 
the scattering of the short-wave radiation by the film. 
Accordingly, for the present we shall have to be content 
with the above results. 

7 Davies, Trans. Faraday Soc. 36, 333 (1940). 
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Fic. 4. Absorption of a-hydroxymyristic acid. (a) 0.0230 
molar solution in carbon tetrachloride with 3.05-mm cell. 
Temperature 64°. (b) Thin film of solid d=0.012 mm. 


components in the absorption : the sharpness and 
symmetry of the 2.834 peak suggests that in the 
molecular orientation (c) both hydroxyl absorp- 
tions coincide, while a structure (b) provides the 
“free,” 2.78u, alcoholic band. 
\ ~ 
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A comparison of the OH absorptions of the 
a-hydroxyacids with say acetic acid under the 
same conditions! shows that the intensity of the 
sharp band is distinctly greater in the former. 
This results from the presence of two OH groups 
in the acid molecule, and the fact that the intra- 
molecular interaction in the a-hydroxyacids 
probably reduces the tendency of the carboxyl 
groups to associate. Unfortunately, owing to their 
sparing solubility, the latter point could not be 
checked by molecular weight determinations in 
carbon tetrachloride 

In the solid both acids show a sharp band 
arising from the a-hydroxyl group. In a-hydroxy- 
myristic acid it is centered at 2.924 with a 
shoulder at 2.93;u, whilst in mandelic acid where 
it is not so sharp, it appears as a doublet peak at 
2.9054, 2.9391 (Av=30 cm). In all cases an 
association band appears on the long wave side. 


DAVIES 


If, as the results for solution suggest, the mole- 
cules are largely in the configuration (c), it is 
quite likely that the alcoholic OH group is little 
disturbed on association of the carboxyl heads 
and, its freedom of rotation both around the 
C—O bond and the C—C chain being removed, 
there is small likelihood of these OH groups in 
different molecules interacting with one another. 
The sharpness of the OH band in the solid then 
follows from the intramolecular ‘‘fixation’’ of the 
hydroxyl group configuration. 

It would be expected that in comparison with 
the acids, the a-hydroxyesters would give even 
clearer results owing to the removal of the free 
OH group of the carboxylic structure. This has 
been found to be the case, and a careful study of 
ethyl mandelate has been made. In dilute carbon 
tetrachloride solution, 0.0232 molar or 0.28 per- 
cent by weight, the ester shows a single sharp 
band at 2.83;u, followed by no association band 
(Fig. 5a). The apparent symmetry of the absorp- 
tion in the most dilute solution suggests that in 
the isolated molecules all the OH groups are in 
one configuration (c, of the acid). The small 
peak at 2.924 (3425 cm~') often appears in this 
neighborhood with hydroxylic derivatives. While 
it is difficult to find conclusive evidence for its 
origin, it may possibly be due to an overtone or 
combination frequency involving the C —O bond. 
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Fic. 5. Ethyl mandelate in carbon tetrachloride solution. 
(a) 0.0232 molar (i.e., 0.28 percent by weight in 3.05-mm 
cell. (b) 2.62 percent by weight in 0.35-mm cell. (c) 47 
percent by weight between two glass strips. 
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Owing to the low melting point of ethyl mande- 
late (34°C), when the thin film of solid was 
placed in the focus of the beam it liquefied. 
Accordingly, by placing the film at some distance 
from the focus it has been possible in this case 
to measure the absorption of the same layer 
both as liquid and solid (Fig. 6a and b). To 
complete the sequence, solutions of 2.62 and 47 
percent by weight in carbon tetrachloride have 
also been examined (Fig. 5b and c). It is first 
clear that ethyl mandelate in the solid gives an 
absorption somewhat similar to triphenylmethyl 
carbinol although broader than the latter: an 
important feature is the absence of an association 
band, in contrast with mandelic acid for which 
this component (arising from the COOH group) 
is very obvious. The free hydroxyl group in 
ethyl mandelate is again protected by its intra- 
molecular interaction, the peak at 2.92 corre- 
sponding to the doublet in triphenylmethyl 
carbinol. The peak at 3.00u is more difficult of 
interpretation. It is clearly present in the liquid 
and probably appears at 2.98u in the concen- 
trated solution. A maximum near 3.00 is found 
in the curves for a considerable number of the 
other compounds studied and it may again be due 
toa common overtone or combination frequency. 

The clearest result of comparing the series of 
curves for ethyl mandelate is the gradual nature 
of the transition from the dissolved to the solid 
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Fic. 6. Absorption of thin film of ethyl mandelate. (a) As 
liquid, (b) as solid, d2=0.016 mm. 
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Fic. 7. Absorption of 6-hydroxyester. (a) 0.0168 molar 
solution in carbon tetrachloride with 3.05-mm cell. (b) Thin 
film of solid. d=0.019 mm. 


state. In the concentrated solution the sharp 
peak characteristic of the free molecules in dilute 
solution is still prominent, although shifted 
slightly to 2.84u: the latter effect is probably 
more apparent than real, being a result of the 
overlapping with the strong absorption on the 
long wave side of the band. At the same time 
the absorptions at 2.94 and 3.0u found in the solid 
are already present. All these features are also 
shown in the liquid, with a change in emphasis 
shown by the increased intensity of the long wave 
components. This continuous change on in- 
creasing condensation of the molecules is in 
agreement with a growing mass of evidence on 
the nature of solutions and the liquid state.® 
The mean degree of aggregation (association) as 
indicated by the appearance of the spectrum 
characteristic of the solid, increases gradually 
through the solutions to the solid, in which 
phase there are no longer any “free” molecules. 
A §-hydroxyester with the structure 8-68-B-- 
phenyl-diphenyl-hydroxyethyl propionate, 


CH; OH 
‘i 


CoH;- CoH, 


8 Trans. Faraday Soc. 33, 1-279 (1937). 


CH.-COOC2H;, 
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Fic. 8. Absorption of thin films of solids. (a) catechol, 
(b) resorcinol. d2=0.074 mm. 


has also been examined in solution and as a 
solid (Fig. 7a and b). In solution the OH absorp- 
tion is shifted to 2.86u, where it is broader than 
most such bands in carbon tetrachloride. That 
the interaction between the OH and C=O is 
greater in this case than for the a-hydroxyester 
may be explained on the basis of observations 
with accurate molecular models for these struc- 
tures. In mandelic ester it is always possible 
easily to get the OH rotating past the O of the 
C=O group of the carboxyl! structure, but this 
is not always possible for the 6-hydroxyester, 
where the interaction can accordingly become 
more pronounced, i.e., the OH and C=O groups 
can get closer to one another. The orientation in 
the y- or higher hydroxyacids can also become 
entirely favorable to the OH and C=O interac- 
tion (but not more so than for the B-hydroxyacid) 
although the a@ priori probability of this con- 
figuration will rapidly decrease. In the solid the 
B-hydroxyester gives the sharpest OH band 
measured, centered at 2.904 with perhaps a 
satellite at 2.944. The exceptional freedom of the 
OH vibration in this solid, as measured by the 
sharpness of the band, agrees with previous 
indications. An important point with the same 
significance is the small shift from dilute solution 
to the solid, only 0.044 or 48 cm compared 
with the mean displacements Av=137 cm~ for 
triphenylmethyl carbinol and Av=108 cm~ for 
the a-hydroxyacids and ester. Clearly the intra- 
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Fic. 9. Absorption of solid films. .(a) hydroquinone, 
d=0.024 mm. (b) Quinhydrone, d=0.016 mm. 


molecular binding and the steric protection of 
the hydroxyl group operate to stabilize the OH 
frequency on passing to the condensed phase 
and both these factors are large for the £- 
hydroxyester. 

Three compounds of interest in the present 
study are the isomeric dihydroxybenzenes. Cate- 
chol, the orthocompound, gives in solution a 
clearly defined absorption with two peaks at 
2.76su and 2.802u.° The results for the solid show 
an absorption starting at 2.8u and rising sharply 
to a main peak at 2.91+0.01y, from which it 
continues, with decreasing intensity however, 
into the region of the CH bands which are just 
discernible at 3.34. By comparison with the 
previous results, the steepness of the curve up to 
and the peak at 2.914 suggest a ‘‘free’’ hydroxy! 
group, the following broad band corresponding to 
associated groups. This is compatible with the 
cis-trans structure (d), apparently the only one 
existing in solution: that is, we deduce that the 


4 


OL. 


(d) 


intramolecular interaction of some of the hy- 
droxyls persists in the solid, the remainder being 
engaged in associative interaction. 


9 Davies, Trans. Faraday Soc. 34, 1427 (1938). 








present 
. Cate- 
tion a 
aks at 
1 show 
harply 
rich it 
wever, 
re just 
h the 
> up to 
droxy!] 
ling to 
th the 
ly one 
at the 


INFRA-RED ABSORPTIONS OF HYDROXYL COMPOUNDS 583 


Resorcinol, or meta-dihydroxybenzene, shows 
a notably different absorption from catechol 
(Fig. 8b). It has been found at 73° in carbon 
tetrachloride, to give a single sharp symmetrical 
peak at 2.77;u, indicating no interaction between 
the hydroxyls. In the solid the absorption starts 
about 2.85u, has an insignificant peak at 2.95yu 
but increases in intensity up to a maximum at 
3.10u. The breadth and distribution of intensity 
in the band correspond to complete association 
of the hydroxyls. This agrees with the x-ray 
results of Robertson and Ubbelohde!® whose 
crystal structure analysis points to “hydrogen 
bridges” throughout the solid. Compared with 
catechol, this is a consequence of the meta- 
position of the two groups preventing any degree 
of internal ‘‘protection” such as occurs for the 
cis-group in that case. 

In hydroquinone, i.e. paradihydroxybenzene, 
the saturated solution in carbon tetrachloride 
shows the same single absorption peak at 2.77u 
as in resorcinol. In the solid the absorption of 
these compounds is also similar, although the 
band appears a little narrower and at slightly 
shorter wave-length in hydroquinone (Fig. 9a). 
There is, however, no doubt that associative 
interaction of the hydroxyls is complete in the 
latter case also. Quinhydrone, C,sH,(OH)>2-- 
CsH,Os, the molecular complex formed by hydro- 
quinone and benzoquinone, has also been ex- 
amined in the solid state (Fig. 9b). In carbon 
tetrachloride it was too sparingly soluble for the 
absorption to be measured. The broad, almost 
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_ Fic. 10. Absorption of 9-10-dihydroxymethy! stearates 
in carbon tetrachloride. 3.05-mm cell. Temperature 51°. 
(a) 0.0173 molar, cis-compound. (b) 0.0219 molar trans- 
compound. 


1938 and Ubbelohde, Proc. Roy. Soc. A167, 122 
( , 


featureless absorption rising to a maximum at 
3.10u resembles that of hydroquinone, or better, 
resorcinol, and the “hydrogen bridges’ must be 
in very much the same state in all these cases. 
With Badger’s form of the relation of inter- 
nuclear distance with force constant," and a 
vibration frequency of 3225 cm (i.e. corre- 
sponding to 3.104) for the associated OH groups, 
the calculated bond distance is 1.02A: this com- 
pares with a distance of 0.97A in the monomeric 
groups. This indicates that even for quinhydrone 
where the solid structure is found to consist of 
infinite chains of quinol-quinone molecules (e), 
the hydrogen bridges are far from symmetrical, 


~_>- w-o— Kona on {=o aan Cl 
(é) 


as the 0—0 distance between adjacent molecules 
can hardly be less than 2.4A. This conclusion 
supports the revised results for the x-ray crystal 
study of quinhydrone.” 

Two long-chain isomeric compounds, the so- 
called cis- and trans-dihydroxybehenic acids™ 
were found by Marsden and Rideal'* to have 
very different properties as monomolecular films. 
As these differences appeared to be related to the 
different orientations of the hydroxyl groups, an 
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Fic. 11. Absorption of solid isomeric 9,10-dihydroxy- 
methyl stearates. (a) Trans-compound, d=0.060 mm. 
(b) Cis-compound. 


1 Badger, J. Chem. Phys. 3, 710 (1935). 

2 Palacios and Foz, Anal. Fis. Quim. 34, 779 (1936). 

13 Hilditch, J. Chem. Soc. (London) 1828 (1926). 
ass and Rideal, J. Chem. Soc. (London) 1163 





584 MANSEL 
infra-red study of their absorptions was of 
interest. In solution and in the solid state the 
isomers resembled one another markedly. The 
results, while clearly interpretable, are some- 
what involved owing to the presence of an 
additional free OH group in the carboxylic 
structure. Accordingly, it seems preferable to 
give the comparable results for the isomeric 
9.10-dihydroxymethyl stearates: CH3(CHs2);-- 
CHOH -CHOH -: (CH2);COOCHs. In carbon tet- 
rachloride the cis-compound shows an absorption 
with two peaks at 2.774 and 2.80yu (Fig. 10a). 
While the curves are not so clear as for the 
benzene derivative (in particular a peak near 
2.864 appears to be real) this presumably corre- 
sponds to a similar state of affairs to that in 
catechol. The fact that the trans-compound in 
solution also shows a doublet peak, perhaps 
best given as 2.76;u, 2.79;u, is confusing only so 
long as these compounds are considered to be 
really cis- and trans-isomers. Actually they are 
diastereoisomers and may be represented : 


(cm,),64, Ho lea), th, 


no 


(f) 


4 


(3) 


The point that emerges is that, insofar as the OH 
absorptions for both compounds are similar, the 
two hydroxyl groups occupy the same relative 
positions in both cases. By analogy with catechol 
where, according to Pauling," the structure (d) is 
necessarily confined to one plane, in this case 
the free rotation about the 9.10 C—C bond may 
bring these groups opposite one another in the 
two esters, as shown in (f) and (g). The mono- 
molecular films may well show differences of 
stability arising from the different ease of 
packing of the two structures (f) and (g). 

The general similarity in the absorptions of 
the solids is again not unexpected on the above 
formulation: Fig. 11a and b. The curves are, 
however, different from that of solid catechol, 
the trans-ester showing a considerable resem- 
blance to hydroquinone in the region 2.9u to 
3.24: in the cis-compound the band is somewhat 


H,), C0064 


EA), corCH, 


% Pauling, J. Am. Chem. Soc. 58, 94 (1936). 


M. 
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narrower and the peaks at 3.01y and 3.09 better 
defined than usual. The differences in the melting 
points and solubilities of these isomers suggest 
different stabilities of their packing in the solid 
and the minor differences in their infra-red 
absorptions can be related to the same change 
in molecular orientation. The fact that for both 
isomers in the solid the absorption is far more 
typical of complete ‘‘hydrogen bonding”’ than in 
the case of catechol may be due to the possi- 
bility of rotation about the C—C bond between 
the two carbon atoms carrying the hydroxy] 
groups. The intramolecular interaction will by no 
means entirely eliminate this possibility and if 
the intermolecular attraction becomes sufficient 
(as it reasonably might in the solid) the two OH 
groups could swing away from one another and 
take part individually in hydrogen bonding with 
other molecules. 


DISCUSSION 


There are two principal features in the fore- 
going results which need further discussion. The 
first is the general shift of the absorptions to 
longer wave-length in the condensed phases: the 
other is the doublet nature of the absorption of 
the ‘‘free’’ hydroxyl groups as observed in the 
solid for triphenylmethy] carbinol, ethyl mande- 
late, etc. 

The former of these effects has long been 
recognized as of very general occurrence. For 
simple molecules, more particularly the hydrogen 
halides, it has been carefully studied by Hettner, 
West and co-workers, and Sutherland, Lee and 
Wu." Thus, for HCI, the fundamental absorption 
centered at 2886 cm™ in the gas at room tem- 
peratures, has its center shifted in carbon tetra- 
chloride solutions to 2833 cm-', in the liquid 
(at 173°K) to 2785 cm, and in the solid 
(at 100°K) it is at 2768 cm—!. The largest change 
we have observed from solution to solid is 
Av=137 cm~ for triphenylmethyl carbinol, and 
the smallest, for the B-hydroxyester, of Av=48 
cm~!, It must be emphasized that these figures 
refer to the sharp absorptions in the solid state 
and not to the change between monomeric and 

16 Hettner, Zeits. f. Physik 89, 234 (1934); West and 
Salant, Phys. Rev. 37, 108 (1931); West and Edwards, 


J. Chem. Phys. 5, 14 (1937); West, J. Chem. Phys. 7, 795 
(1939); Sutherland, Lee and Wu, Nature 142, 669 (1938). 
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associated structures, for which Av can be ob- 
served in the gas, solution, liquid, or solid phases 
to be of the order of 350 cm~', depending on the 
type of association. Actually, of course, the 
association band itself. shifts (e.g., for cetyl 
alcohol or hydroquinone) by about 140 cm~! to 
lower frequencies on passing from solution to 
the solid. 

Of the three general sources of interaction 
between saturated molecules, that is, the dipole- 
dipole, dipole-induced dipole, and dispersion 
effects, calculation for molecules as close together 
and maintaining a fixed orientation as are those 
in the solid state, shows that the first of these 
items largely predominates for polar molecules.’ 
Accordingly, the purely dipole forces, which tend 
to assume a directed character in the solid phase, 
would seem to be mainly responsible for the 
perturbation in the OH bond vibration: that is, 
the molecular dipole field in the solid is sufficient 
to stretch the hydroxyl dipole the small amount 
necessary for the observed lowering of its fre- 
quency. On the basis of the Badger-Clark rule, 
the change in bond length involved for tri- 
phenylmethyl carbinol is 0.012A, in a bond of 
0.98A, i.e. a matter of 1.2 percent. Similarly, for 
HCI in passing from the gas to solid, the bond 
length increases 1.4 percent. (If this change were 
effected on isolated molecules in the gas it can 
be estimated, on the assumption of Hooke’s law 
applying to the extension, that it would necessi- 
tate an energy of 0.0047 ev.) That electro- 
static forces are largely responsible for these 
changes is borne out by the following com- 
parisons: 

RO-H HCl HBr C-H 


Avgas solid Cm™ 180 118° 93 19 
xin Debyes, 1.7 1.03 0.78 0.3. 


This conclusion is in agreement with that of 
West and Edwards, whose paper should be con- 
sulted for a discussion of the other possibilities 
and a calculation of the dipole distortion made 
by Kirkwood. It also provides an explanation of 
the fact that the greater the intramolecular 
binding of the hydrogen, the smaller the fre- 
quency change on passing from solution to the 
solid. As a measure of the former we may take 
the displacement of the main OH frequency from 
its expected position in solution (dv): Av is the 


subsequent change on passing to the solid: 
Triphenyl- 
methyl Ethyl B-hydroxy- 
carbinol Catechol mandelate ester 
bv 0 70 92 132 
Ap 137 115 108 48 


Clearly, the more firmly the hydrogen of the OH 
group is held internally by dipole interaction, the 
smaller will be its external field, with a corre- 
sponding decrease in intermolecular perturbation 
in the solid. 

The doublet nature of the relatively sharp 
bands of a number of simple vibrators in the 
condensed state is also of fairly general occur- 
rence. That the effect is not always confined to 
this phase is most clearly shown by hydroxylic 
compounds. Naherniac” in an examination of the 
second overtone of the valence OH vibration in 
a number of alcohols in the vapor and liquid 
states reports a doublet peak for twelve out of 
fifteen alcohols examined in the gas phase. 
Badger and Bauer® had earlier found the same 
doublet structure for three primary alcohols in 
the vapor. The frequency separation in these 
cases is usually about 50 cm~. For triphenyl- 
methyl carbinol, mandelic acid and a 8-hydroxy- 
ester, we have found doublet separations of 42, 
30, and 46 cm in the solid. Two possible 
interpretations of this splitting may be con- 
sidered. (i) The peaks may arise from the 
occurrence of the hydroxyl groups in two states 
owing to varying degrees of interaction with 
adjacent groups, either within the same molecule 
or in neighboring molecules. (ii) The doublet 
may correspond to a combination of the OH 
valence vibration with a low frequency vibration, 
or libration, of the hydrogen in this group, 
probably in a direction roughly perpendicular to 
the OH and C—O bonds. These two factors 
almost certainly operate separately in various 
hydroxylic structures. The former is equivalent 
to Pauling’s interpretation of the infra-red results 
for orthochlorphenol, etc.'® and is the explanation 
favored by Badger and Zumwalt"® for the slitting 
of the OH frequency in simple alcohols. Methyl 
alcohol is undoubtedly a special case in this 


17 Naherniac, Ann. d. Physik 7, 528 (1937). 
18 Badger and Bauer, J. Chem. Phys. 4, 711 (1936). 
19 Badger and Zumwalt, J. Chem. Phys. 7, 87 (1939). 
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respect; thus the two components reported by 
Naherniac with Av only 30 cm refer to promi- 
nent features in the fine structure of the band, 
resolved, owing to the smallness of the moments 
of inertia of this molecule, by Badger and 
Bauer.”° Borden and Barker*! have shown that 
the hydroxyl group is not free to rotate, at least 
in the ground state in methyl alcohol, and it 
seems likely that the restricting potential barrier, 
due to dipole and other interactions between the 
OH and adjacent CH groups, will increase for 
the higher homologues. 

The second possibility has been used in 
accounting for the great breadth of the associa- 
tion band in alcohols and carboxylic acids,! and 
one important confirmation of its operation 
seems worth mentioning. It is to be expected 
that in the carboxylic dimer combination of the 
valence OH vibration will be most pronounced 
with “‘deformation”’ vibrations involving the OH 
groups. The frequencies of the latter will be 
notably decreased by substitution of deuterium, 
which should lead to a narrowing of the associa- 
tion band.* Wall and Claussen”? have recently 
published results for the association band in both 
deutero- and normal acetic, propionic, and 
benzoic acids. They conclude that in the hydro- 
gen acids the association band lies between 
3.24 and 4.0u (Av=625 cm~!), in the deutero- 
acids between 4.34 and 4.94 (Av=289 cm~). 
The difference is most pronounced and in the 
expected direction. 

As the doublet structure both in the vapor 
and in the solid corresponds to quite narrow 
bands and the second factor results in a very 
broad absorption, it seems that the first interpre- 
tation will be the more satisfactory in the 
present circumstances. Some additional evidence 
in favor of this view may be mentioned. 
Naherniac” fails to comment on the fact that 
while the alcohols show this doublet structure, 
of six carboxylic acids examined in the vapor 
phase only one, and that formic acid, shows it 
for the COOH structure. In this grouping, 


20 Badger and Bauer, J. Chem. Phys. 4, 469 (1936). 
21 Borden and Barker, J. Chem. Phys. 6, 553 (1938). 
*This point was raised in a discussion with Dr. 
Sutherland. 
Wall and Claussen, J. Am. Chem. Soc. 61, 2812 


(1939) 
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however, the hydrogen is tied so firmly in the 
orientation 


O 
VA 


—C H 
a ed 
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(with an estimated energy of 7400 cal. per g mole 
relative to the trans-structure), 


O 


—-C 
\ 
1 
H 
that the restoring force for the OH valence 
vibration must be uniquely determined. The 
anomaly of the first member of the series again 
results from the extended structure of the band 
due to the low moments of inertia of the molecule, 
as is proved by a comparison of the high dis- 
persion work on formic and acetic acid vapors.” 
Closely similar to this point is the observation 
that for the two molecules where the intra- 
molecular field is sufficient to tie the OH groups 
all in one orientation in solution and which give 
a ‘free’ hydroxylic absorption in the solid, i.e., 
ethyl mandelate and the 8-hydroxyester, the 
sharp absorption in the solid was not clearly 
resolved into a doublet. 

While remembering that we are here dealing 
with the free hydroxyl in contrast to the associ- 
ated group, it is not immediately obvious why 
the doublet structure is not generally observed 
for the alcohols in solution or in the liquid phase. 
In carbon tetrachloride, where the majority of 
the reliable measurements have been made, the 
monomeric band is usually a single peak of 
half-width about 60 cm! shifted by a similar 
amount from its position in the vapor phase. 
Presumably owing to the constantly varying 
interaction with adjacent molecules in this phase, 
the potential field to which the OH group is 
subject takes on a continuous series of values 
over a range determining the position of the 
band. That such an interaction with solvent 
molecules occurs is highly probable from the ob- 


8 Badger and Bauer, J. Chem. Phys. 5, 605 and 852 
(1937). 
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servation of the apparent doublet nature of the 
alcoholic OH absorption in benzene solutions. 
This is shown by the corrected absorption curve 
for a dilute solution of benzyl alcohol in benzene, 
Fig. 12. The center of the absorption maximum 
is now at 2.80u (cf. 2.7644 in CCly, 2.770u in 
CHCl;) and the peak is no longer symmetrical: 
while sharp on its short wave side, a distinct 
shoulder is shown at 2.815, not resolvable as a 
separate peak with the quartz prism. In this 
solvent, owing to the ready polarizability of the 
double bonds of the benzene ring, the dipole- 
induced dipole interaction, added to the always- 
present dispersion effect, may become quite 
considerable. The result is that the distribution 
of the solute molecules between at least two 
energy states (corresponding to different degrees 
of interaction with the solvent) becomes detect- 


able. The splitting of the OH frequency in this 
way has also been reported by Naherniac." It is 
likewise very apparent for HCI solutions® where 
its magnitude varies with the dielectric constant 
of the medium. 

Obviously the interaction of the OH with 
adjacent groups will be favored by the fixed 
orientations and close packing in the solid: thus 
the fixed orientation of the dipoles with respect 
to one another means that their interaction ap- 
proaches the maximum possible (U = — 2yusu2/ R*) 
where the dependence is upon R= rather than 
R- as occurs for random orientation at large 
distances. Taking won =1.6D and wou =0.3D, at 
a separation of 3.5A, which is a reasonable esti- 
mate for groups in neighboring molecules in the 
solid, these dipoles give U= — 325 cal. per g mole. 
The frequency separation of Av~40 cm~ for the 
doublet in the solid may be roughly estimated on 
the empirical proportionality of energy difference 
and frequency separation’ to correspond to an 
energy of 660 cal. per g mole. Certainly other 
contributions to the interaction energy should be 
added to the simple dipole effect: especially in 
structures with phenyl groups the dipole-induced 
dipole interaction may be considerable. (This 
should lead to a larger doublet separation in 
phenolic vapors than in aliphatic alcohols. For 
resorcinol, where the groups can hardly interact 
being in the meta-position, Naherniac gives 
Av=104 cm, but reports only one peak for 
phenol.) The occurrence of a doublet peak for 
the alcohols is then throughout referred to the 
same cause, the only difference being that in 
the vapor the interaction is intramolecular 
whereas when it occurs in solution and the solid 
it is more likely intermolecular. 

The author wishes to thank Professor E. K. 
Rideal, Professor T. P. Hilditch, and Dr. P. 
Maitland for specimens of some of the com- 
pounds studied, and he gratefully acknowledges 
the helpful discussion of this paper with Dr. G. B. 
B. M. Sutherland. The work was carried out 
during the tenure of a Senior Award of the 
D.S.I.R. and was aided by a grant from the 
Chemical Society. 
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The behavior of isothermal packed fractionating columns is deduced from their differential 
equation. The relation between rates of production and over-all fractionation is studied, as is 
the rate of approach to equilibrium. The results are illustrated by problems of isotope separation. 





HE theory of plate columns has been exten- 

sively treated. The behavior of packed 
columns is generally related to that of plate 
columns by determining the “height equivalent 
to a theoretical plate.’ 

For certain problems. arising in the use of 
packed columns equivalent to hundreds of plates, 
the graphical and algebraic methods developed 
for small plate columns are inconvenient. Other 
problems, of importance for long columns only, 
such as time phenomena, have not yet been 
treated. 

The present paper discusses these problems, 
starting from the consideration of a differential 
element of the column. The two steps in the 
theory of plate or packed columns, i.e., the deter- 
mining of the relation between liquid velocities 
and surface and the number of theoretical plates, 
and then the finding of the properties of a column 
of n theoretical plates and given reflux ratio, are 
not separated. The solution compares favorably 
in simplicity with the second part of the usual 
procedure taken by itself. 

We restrict ourselves to isothermal columns 
with constant liquid and gaseous hold-ups. Even 
for isotope concentration columns, these condi- 
tions are not strictly true, since the gaseous hold- 
up and to a smaller degree the liquid hold-up 
are affected by the pressure-drop across the 
columns. In the fractional distillation of two 
closely-related liquids, the hold-ups will depend 
on the equations of heat transfer. The idealiza- 
tions of constant hold-ups and temperature 
could, of course, be obtained in practice: in the 


* Publication assisted by the Ernest Kempton Adams 
Fund for Physical Research of Columbia University. 

1C. S. Robinson and E. R. Gilliland, Elements of Frac- 
tional Distillation (McGraw-Hill, New York, 1939) (con- 
tains bibliography). 

2 J. R. Huffman and H. C. Urey, J. Ind. Eng. Chem. 29, 
531 (1937). 


former case by compressors, in the latter case, by 
appropriate heating or cooling. The reality of 
these conditions in a particular case is, however, 
inessential. We consider here only the kernel of 
the problem. But the modifications necessary for 
any actual process will, we believe, be fairly 
obvious. 


THE DIFFERENTIAL EQUATION 


Refer to Fig. 1. Consider a two-component 
two-phase* system. One component will be 
described as ‘“‘heavy,’’ the other, ‘‘light.’’ We will 
take N to be the mole fraction of heavy material 
in the liquid, [N] the corresponding concentra- 
tion (moles/liter) ; 1— N the mole fraction of light 
material in the liquid, [1/7] its concentration; 
H the hold-up per unit length (moles/cm) of 
combined heavy and light material in the liquid; 
L the flow (moles/sec.) of combined material in 
the liquid : and let lower-case symbols denote the 
corresponding quantities in the gaseous phase. Z 
is the length of the column. z is the distance along 
the column, measured from the “‘condenser.”’ 

If HéN/ét be the rate of transfer of heav 
material to the liquid across the liquid-gaseous 
interface, the conservation of matter gives us 


HdN/dt= —LON/d2+H6N/6t 
hdn/dt=ldn/dz—HbN/ ot. 


We will take 
HbN/st= —k{(N][m]—o[M][n]}. 


We will justify the use of (1.1) later. At this 
point we only remark that it is of the same form 
as 6[ NV ]/ét for a bimolecular reversible reaction. 
k is a constant depending on the permeability of 


(1.1 


* Throughout we take these two phases to be liquid anc 
gas. Obviously the two phases could equally be two im- 
miscible liquids (extraction), etc. 
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the interface, a is either the chemical equilibrium 
constant of a reaction such as 


HC?!8N(g) +C!?N~-(aq.) 
=HC"N(g)+C'!8N-(aq.) (1.2) 


(for the case of the concentration of isotopes), or, 
should we be fractionating two closely related 
liquids such as isomeric heptanes, the ratio of the 
partial pressures. In either case, at equilibrium 





sey ES 
i-nN/ 1-n [M][n] 


Substituting for 76N/ét from (1.1), and intro- 
ducing mole fractions, the differential equations 
governing the behavior of the column become 


ON ON 
H-—— = —L—-—kcC{ N(1—n) —an(1—N)}, (1.4) 
ot Os 


On On 
h—=Il—+kceC| N(1—n)—an(1—N)}, (1.5) 
ot az 


where 


c=[m]+[n]=a constant 
C=[M]+[N]=a constant. 


Equations (1.4) and (1.5) are nonlinear, and, 
except for the stationary state, cannot easily 
be solved. Considering N as always small, and 
thus writing 1 for 1— J, is not justified; the en- 


(1.6) 





2-0 rR 











a-Z ree 





Drain HCN 
Generator 


Fic. 1. Schematic diagram of fractionating column. 
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tire purpose of the use of columns is to obtain N’s 
which are not small. Let us introduce R= N/1—N 
and r=n/1—n. Then (1.4) and (1.5) transform 
into 





OR oR 1+R 
H—= —L——kcC- (R-ar), (1.7) 
at Oz 1+r 
or or 1+r7 
h—=Il—+kcC——(R-—ar). (1.8) 
ot as 1+R 


Since R/r varies only between 1 and a, we 
have, to terms of the order of (1—«a)* 


OR OR 
H—= —L—-—kcC(R-ar), (1.9) 
at Oz 
or or 
h—=l—+kcC(R—-ar). (1.10) 
Oot O02 


In the steady state (0R/dt=dr/dt=0), (1.9) 
and (1.10) give 
lar /dz— LOR/dz=0 
or 
lr—LR=constant, (1.11) 


whereas the exact Eqs. (1.4) and (1.5) possess as 
their integral of conservation 


ln — LN=constant. (1.12) 


This discrepancy is unimportant so far as the 
rate of approach to equilibrium is concerned. 
But since the question of conservation of matter 
is of the utmost importance in the separation of 
isotopes, in addition to solving the approximate 
Eqs. (1.9) and (1.10) for the steady state, which 
gives easily the general sweep of the relationship 
between over-all fractionation and production, 
we shall also solve the rigorous Eqs. (1.4) and 
(1.5) for the steady state. From this latter solu- 
tion we can find in any particular case the exact 
relationship between fractionation and _pro- 
duction. 


THE STEADY STATE: PRODUCTION AND FRACTION- 
ATION (APPROXIMATE SOLUTION) 


The set-up is essentially as shown in Fig.1. 
In the case of a fractional distillation, the NaOH 
solution which here absorbs the HCN gas and 
transforms it into NaCN solution is replaced by 
a total condenser. In the case of the concentration 
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of NH3, the roles of the liquid and gas are 
interchanged, and the NaOH solution serves to 
strip the ammonia from the liquid. In all cases 
the boundary conditions are 


(2.1) 
(2.2) 


At 2=Z :r=p=normal abundance ratio 
Atz=0:r=R. 


If P be the production in moles/sec. of mixed 
gases, 


1=P+L. (2.3) 


Introducing p= P//=fraction of flow drawn off, 
we find 


L/l=1-—p. 
We are now to solve (1.9) and (1.10) with 
dR/dt=dr/dt=0, subject to the boundary condi- 


tions (2.1) and (2.2). The solution is easily 
found to be 


(2.4) 


e~*?-+ a6 
~ L— pewZ 40 
(1—p)e*?+6 
(= pe##+0 
where 6=p/(1—a) and 
u=(kcC/L)(1—a)(1+26). 


When there is no production, (2.5) and (2.6) 
reduce to 





(2.5) 


(2.6) 





r 


R=r= pe™*(1-2/Z), (2.7) 


introducing the quantity o, which will be useful 
to us later. 


kcC 


o=—(1-—a)Z. (2.8) 
2L 


COHEN 


The over-all fractionation is given by 


1 kceCZ/L 
-) . 
a 


the approximation being valid when a is near 1. 

For the system defined by (1.2) and (1.3), 
a<i and the column therefore accumulates 
heavy carbon at the production end. Since R=r, 
(2.7) is an exact solution of the rigorous Eqs. 
(1.7) and (1.8) or (1.4) and (1.5) for no produc- 
tion. The expression (2.9) for the over-all 
fractionation agrees with the well-known formula 
for plate columns, with kcCZ/L the number of 
theoretical plates. We conclude that (1.4) and 
(1.5) are of the correct form, and therefore that 
H6N/ét is correctly given by (1.1). HéN/ét is 
generally given as* HéN/ét=k’(N—N*) (for a 
liquid-film-controlling transfer) or HéN/ét=k”’ 
(n—n*) (for a gaseous-film-controlling transfer), 
where N* and n* are the concentrations in equi- 
librium with the opposing phase. These expres- 
sions depend on the distance from equilibrium in 
the same way (to terms in (1—q@)?) as does (1.1). 
(1.1) is preferable since k (but neither 2’ nor k’’) 
is independent of concentration. 

Let us see whether (2.8) cannot be verified. 
The reaction represented schematically by (1.2) 
is composed of two parallel reactions. 


(HCN (g)+C?N-(aq.) 
=HC!N(g)+C!N-(aq.) 
| HC™N(g)-+HC#®N(aq.) 
L =HC"N(g)+HC!®N (aq.), 


(2.9) 


(2.10) 





whose equilibrium coefficients we will take, i 


TABLE II. 








EXPERI- 
MENTAL 


0.0752: ° 
ae.is” 
a 


(Pg)opt 
FROM (2.12) 


0.063 atmos. 


REACTION 


NHsNO;—NH; 19. 
NaCN—HCN 4. 
a 
Be 


C1 





NaHSO;—SOz 
NaHCO;— CO, 




















* At the foot of the columns. 

2H. C. Urey, J. R. Huffman, H. G. Thode and M. Fox, J. Chem. 
Phys. 5, 11 (1937). 

+H. G. Thode and H. C. Urey, J. Chem. Phys. 7, 1 (1939). - 

¢I. Roberts, H. G. Thode and H. C. Urey, J. Chem. Phys. 7, 13: 
(1939). 

4C, A. Hutchison, D. W. Stewart and H. C. Urey, reference 6. 

¢ Results to be published. 


3 For example, reference 1, p. 127. 
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round numbers, as 0.98 and 1.005, respectively. 
Supposing that the rate of transfer is controlled 
by a relatively slow diffusion across a gaseous 
interface, k will be proportional to the self- 
diffusion constant of HCN, and we may consider 
the reaction in solution to be sufficiently defined 
by 





0.02 Cen- — 0.005 Cren (aq.) 
Con-+ Cuen (aq. 


The C’s represent concentrations. 

Diffusion constants vary inversely as the pres- 
sure. Taking account of the presence of water 
vapor in the simplest possible way, 


k~1/(pucn+ pnz.c). 


Remembering that c= pucn/22.4 where pucn 
is measured in atmospheres, using the same 
value for Cucn(aq.) aS for pure water (namely 
Cucn(aq.) =10 Pucn), and taking Cceyn-=4.7 ac- 
cording to our usual practice, 


4.7Z Pucn 


x 
44.8. pucn +0.03 
0.02 X4.7 — 0.005 X 10pucn 
x 


4.7 + 10pucNn 


1-—a= 





o=const. X 





2s 6789) 
x/ 


Fic. 2. Fractionation vs. production. e®’ is fractionation at 
no production. 
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Fic. 3. Ratio of abundance ratios in liquid and in gas at 
foot of column, vs. production. Transports are proportional 
to deviation of curves from unity. 


has a maximum at pucn=9.1 atmosphere. AIl- 
though o has a flat maximum, the fractionation 
e** varies very considerably with the pressure. 
Table I shows the relationship between fractiona- 
tion and pressure according to (2.11) for a 
column whose fractionation factor is 30 at 0.1 
atmosphere partial pressure of HCN. The pres- 
sure in the system is therefore quite critical. 

In general, for any system similar to (2.10), if 
C, represents the concentration of the ions in 
solution, ay,o0 the equilibrium constant of the 
exchange reaction between the gas and its water 
solution, @ion the exchange equilibrium constant 
between the gas and the ions, m, the molarity of 
the gas in solution under partial pressure p,, 
w= —(1—a)n,0/(1—a@) ion and f=m,/p,, then 
the pressure at which fractionation will be a 
maximum is given very nearly by 


be=0.17 Ci/(fCi(1+w))}. (2.12) 


It will be observed that p, depends only slightly 
on w (which is of course less than 1). Taking 
w=}, taking for f its limiting value at zero pres- 
sure and zero concentration of salt‘ and for C, the 
concentrations used in experiments in these 
laboratories, we obtain Table II. In the last 
column are given the experimental optimum 
pressures, determined by trial-and-error. For 
those reactions which have been successfully used 
in the production of isotopes, the agreement is all 
that could be desired. 

Returning now to the general formulae (2.5) 


4 All values from Lewis and Randall’s Thermodynamics 
(1923). 
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and (2.6), p¥0 
0) (1p) +0 
r(Z) (1—p)(e*) +246 
R(Z) — (e“*)#+#9-4-08 
r(Z) (1—pyle™) 409 


r(0)/r(Z) is practically a function of e** and @ 
alone. Only for those values of 6 near 1 does the 
value of a begin to count. For concreteness, we 
take a=0.98 and we obtain for various values of 
e** the curves shown in Fig. 2. 

We see from these curves that even a very 
small production (if 6=10- the fraction of flow 
drawn off is 0.02 percent), reduces the fractiona- 
tion very markedly. The obverse of this phe- 


(2.13) 





(2.14) 








KARL COHEN 


nomenon is that the transport, the amount of 
heavy material saved from the drain, is practically 
zero for small 6 (large fractionation) and is only 
at its maximum for large 6 (fractionation nearly 
1). The transport is proportional to 


1—R(Z)/r(Z) 


and is readily followed from Fig. 3, which gives 
R(Z)/r(Z) against 6 from Eq. (2.14). 
From (2.13) 
_{1+0/—p) 8 
r(0)/r(Z) 1-p 


1/(1+ a8) | 


(2.15) 


e~ 20 





This gives the necessary value of o to obtain a 
fractionation r(0)/r(Z) and a production p. 


THE STEADY STATE: PRODUCTION AND FRACTIONATION (Exact SOLUTION) 


We return to (1.4) and (1.5), 0N/dt=0n/dt=0. The integral of conservation is 
—(1-—p)N+n=pno, 


where %p is the value of m at z=0. Solving for N, and inserting in (1.5) we find 


dn kcC 


Yigg | — pnotn[(1—a)(1+ pmo) +ap]+n*(a—1)} 


dz l(1—p) 
kcC 


dn 





_ dz= 
K(1—p) 


a+bn+cn? 


with obvious definitions for a, 6 and c. The integral of (3.3) is 


kcC 


— z+constant = tanh7! 


W(1—p) 


-) 2cn+b 


(3.4) 
(—D)! 


(—D)} 


where D=4ac—b?. Only values of 6<.05 are of interest. With this restriction, to terms of the order 


of (1—a)? 


b=1-—a+p 
= —(1l-—a) 
(—D)'=(1—a){1+0—2n,8} 


and (3.4) becomes 


1+0@—2n=(1+6—2n,0) tanh | 


(1—a@)(14+0—2n8)kcCz 


(3.5 





+const. , 
2l 


When 2=Z, n=£=mole fraction of heavy isotope in ordinary matter. This determines the constant. 


Therefore 


(1—a)(1+0—2n08)kcC(z—Z) 


+tanh-! (3.6 





1+6—2n=(1+6—2n,8) tanh 


1+6-—2é ! 
21 1+0—2n,0) 


The mo in (3.6), introduced by (3.1), is as yet unknown. However, taking z=0, (3.6) becomes an 





unt of 
ctically 
is only 

nearly 


(2.15) 


tain a 


p.- 
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implicit equation for mo: 


1+6—2n )=(1+6—2n,0) tanh } —o(1+6—2n,.0)+tanh-! 


1+0—2é 


————— }. (3.7) 
1+6—2n08 


(3.7) contains all the information we wish to know, and its solutions for various values of 6 give 
the relationship between fractionation and production. For 6=0 we find 


No/1 —Npo 


—__ = 20 


g/1—¢ 


which is precisely (2.9). The other values must be found by computation. The roots have been found 
for e*= 30, €=.01; r(0)/r(Z) calculated from them, and the results plotted as the heavy solid line 
on Fig. 2. Taking e?’=702, £=0.0038 gives a more striking difference between 7(0)/r(Z) calculated 
from (2.13) and r(0)/r(Z) calculated from (3.7) (Table III). And in fact it can be shown that (2.13) 
becomes a worse approximation as pe” increases: it ceases to be useful about pe?” =0.5. 


THE RATE OF APPROACH TO EQUILIBRIUM 


We take the production zero. (2.4) becomes L//=1. To the boundary conditions (2.1) and (2.2) 
we add 
Att=0: R=r=p. (4.1) 
(1.10) may be written as 


(4.2) 


ot hoz 


h (or lor 
R= art —| 
kcC 
and introducing in (1.9) we obtain 

. ae 0°*r ’r or or 
—(v—V) — Vs—= —A—+BV—. 
ot? Otdz 02" ot Oz 


v, V, A and B are defined by 
kcC 

A=—(h+Ha) V=L/H, 
hH 


kcC 
me v=l/h. 





Eliminating R from the boundary conditions, these become 


At s=0: dr/dt—vdr/dz=Br, 
(4.5) 
Ats=Z: r=p. 
Letsteaieail At t=0: r=p, or/dz=0, dr/dt=Bo. (4.6) 
ntroducing 


G(z, p= f e~?'r(z, t)dt, (4.7) 


(4.3) becomes® 





P p(p+D) (—-V)pt+BVdG Vv @G 
P+Ap p-+Ap dz pr+Ap dz?’ 


* This parametric p will not be confused with the p defined in (2.4). 

* A problem perfectly analogous to this one, the time dependent equation of the Clusius column, has been treated by 
J. Bardeen, Phys. Rev. 57, 35 (1940). The various steps of the operational method, which are sketched in here, are there 
exposed in more detail. 


(4.8) 
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calling 

A+B=D=(kcC/hH)(h+H). 
The boundary conditions transform into 


—vdG/dz=(B—p)G+p atz=0 
Gp=p at s=Z. 


The solution of (4.8) satisfying the boundary conditions (4.9) is 
ies p(p+D) — p(BD/v)e-*°*y—! sinh vs eS. Z) 
p+A (p+A)g(p, Z) 
(v—V)p+BV 
Ey= ’ (4.11) 
2Vv 
(v+V)p—BV 
E\= : (4.12) 
2Vu 
f(b) = — pBe*®*2+ p(BD/v)y“' sinh yZ, (4.13) 
g(p, 2) = Eyy~! sinh yz+ cosh yz, (4.14) 
(v—-V)p+BV)? p?+Ap7} 
=4| hace 
2Vv 





(4.10) 


where 











(4.15) 
Vv 


And therefore, by the Heaviside expansion theorem 


f(bi)g( Pi, 2) 


r= per7(t-2/Z) 4 ePit—Eolp;) 2 (4.16) 


j (pj +A)pig' (pi, Z) 





where the p; are the roots of 
g(p, Z)=0. (4.17) 


For any column, v>V and A>B. Introducing the quantities Y=1/(1—a@) and e=BV/v=2cL/HZ, 
and neglecting terms of the order of B/A or V/v 


Eo=(0/Z)(p/e+1), (4.18) 
E,=(¢/Z)(p/e—1), (4.19) 


ofp? p ‘ 
y=2(= +4741) (4.20) 
Z\é € 


and (4.17) becomes, discarding a root at y=0 


2 } 2/244 +1)3 
tanh o(—+40-+1) == (ets Weie+1) " (4.21) 
é? € 1—p/e 





(4.21) has no positive roots. The smallest (and most significant) roots of (4.21) may best be obtained 
by successive approximations. One solves first the equation 


p” p° 4 p” 0 4 
tanh o(—+4—-+1) = (= +441) , 
€ € 


€ € 
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then the equation 





bY BNE (eh /e + 4bpl/et 1)! 
tanh o(—+4v—+1) = , 
2 1— p°/e 


€ € 
etc. The convergence is excellent : p°— p! is in the third significant figure for y= 80. The first approxi- 
mation gives py/e as a function of o alone. 

For all ordinary purposes, only the smallest negative root is necessary, and the first approximation 
to this will be ample. Figure 4 gives a plot of p:y~/e against o from (4.22). 

No simple curve can be constructed to read off the coefficient f(p1)/(p1 +A) pig’ (pi, Z) since this 
depends on y as well as o. The coefficients of the immediately ensuing terms can, however, be ade- 
quately approximated to by p/ap;. 

Figure 4 gives us immediately the most important information about the approach to equilibrium. 
From (4.16), setting z=0 

r(0) = pe2® — Cye—!”i!t— Cae! melt... (4.23) 


where the C,, C2:-- are positive constants, and 


p(e’*’—1) >C\y>C.>C;:-- 
|Pil|<|p2| <|psl---. 


r(0) is given by a fixed term pe”’ and a series of transient terms, of which — C,e~'”:'‘ is the largest and 
the slowest to vanish. The time of coming to equilibrium is therefore determined by /,. 

Let us define the quantity r= —1/p,, the ‘‘period’”’ of the column: it is the time required for the 
principal transient term to diminish to 1/eth of its initial value. p:= —eF(c)/y where F(c) is the 
ordinate in Fig. 4. Therefore 

y H 


r= = (4.24) 
eF(c) (1—a)*vL F(o) 





with y=kcC/L=number of plates/unit length (2.9). Huffman and Urey? found, by other considera- 
tions, the same dependence on H, a, v, and L. 
For concreteness, we choose the figures 


L cm 2o 20V 
=50, —=V=0.238—, —=1.25X10-3/cm, «= =25.9deg.-'=1.08hour—!. (4.25) 
l-—a H sec. Z Z 


Holding y and ¢ (thus V and the number of plates per unit length) constant, we can vary o by 
varying the length of the column (2.8) (Table IV). For all practical purposes equilibrium will have 
been reached at 3r. 

To find the precise form of an r—¢ curve, and as an example of a detailed solution, we have com- 
puted the first five terms of (4.16), with the roots correct to the third significant figure for Y= 1/0.0126 
=79.4; «=25.9 day; ¢=0.90, e??=6.05. The result is 


r(0) 
—— = 6.05 —4.61e—% 9982 — 0), 18¢-2: 11t — 0. 05-8 06 — 0,0 3e-12-0t —. . (4.26) 


p 
Substituting (4.16) into (4.2), and evaluating, we find likewise 


R(Z) 
—— = 1—0.0139e—°- 9982¢ +. 9.0025e—2: 113! —0,0014e—* 964. .-, (4.27) 
p 
In both of these ¢ is measured in days. It will be seen that after the first day, only the first two terms 
contribute to 7(0) or R(Z). (4.26) is plotted on Fig. 5, (4.27) on Fig. 6. The transport, which is pro- 
portional to 1— R(Z)/p, is also represented on Fig. 6. 
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It should be noticed that with the aid of the 
relations between production and fractionation 
it is possible to determine the time necessary to 
attain a given rate of production. For purposes of 
illustration, suppose we want to find how long 
we must wait before the sixfold column con- 
sidered in Fig. 5 will produce at maximum frac- 
tionation with reflux ratio p=1X10~. In this 
case 0= 1X 107*/1.26X10-?=7.9 X10~, and from 
Fig. 2 we find r(0)/r(Z) =5. This is the maximum 
possible fractionation with the given reflux ratio. 
And, reading off the value from Fig. 5, we see that 
it requires 14 days to attain this fractionation. 


INTERPRETATION OF EXPERIMENTS 


In laboratory procedure a single long column 
(say 30X) is never used since (Table IV) its 
approach to equilibrium is too slow. Thus the 
best experimental material is on compound 
columns of two or three units. Although there is 
no conceptual difficulty in treating such com- 
pound columns as above, the amount of com- 
putation would be discouraging, while there is 
no reason to expect that we could gain much 
insight into the behavior of the columns. The 
following comparisons of properties of compound 
columns with the above results on simple columns 
are therefore made with due reserve. 


‘th pte F(o) 
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Fic. 4. —piy/e= F(c) vs. o. p; is the smallest negative root 
of (4.21). 
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The approach to equilibrium is found® (Fig. 7) 
to be adequately described by a straight-line 
relationship between r(0)/p and ¢ right up to 
the maximum value. Figure 5 indicates that 
equilibrium is approached very slowly after an 
initial steep rise. The discrepancy is undoubtedly 
caused by experimental fluctuations in L, which 
are of the order of 5 percent, and errors in deter- 
mination of 7(0)/p. Clearly, experimental equi- 
librium will be reached when fluctuations in 
r(0)/p from these causes is about double the 
theoretical daily increase in 7(0)/p. This point 
would be reached, on Fig. 5, after about 15 days. 
The column would then be at only 5/6 of its 
true equilibrium. Up to this point, the curve is 
sufficiently steep to give an experimental straight 
line. Taking 33 as the experimental equilibrium, 
the true equilibrium is near 33 X (6/5) = 40. 

It was found® that drawing off 0.046 mole of 
HCN per day resulted in a product for which 
r(0)/p=28. Taking 1—a=0.026, we find 6=1.9 
X10-*. From the exact curve on Fig. 2, we see 
that for 6=0.019, 7(0)/p is reduced from 30 to 
21.3. If this ratio were maintained, it would 
indicate that for @6=0 the HCN units would 
come to equilibrium at 7(0)/p=39. This is in 
good agreement with the estimate of true equi- 
librium made in the preceding paragraph. 

It has been the practice, in these laboratories, 
to take the transport as invariable, proportional 
to 1—a. Knowing the distribution in the columns, 
and having determined experimentally the time 
necessary to arrive at this distribution, the 
transport per day and thus 1—a was calculated. 
It is clear from Fig. 6 that the transport is 
always less than 1—a, and the average efficiency 
of the column there considered for the first 15 
days is 63 percent. The values of 1—a calculated 
from the transport must represent an effective 


TABLE IV. 








e2? = FRAC- 


TIONATION 7 = PERIOD 





is 3.5 hours 
2. 12 hours 
4, i 2.8 days 
5. F 4.1 days 
6. 5.7 days 
0. 


3 47 days 

















®C, A. Hutchison, D. W. Stewart and H. C. Urey, J. 
Chem. Phys. 8, 532 (1940). 
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Fic. 5. Approach to equilibrium of a typical column: frac- 
tionation vs. time. 
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Fic. 6. Ratio of abundance ratios in liquid and in gas at 
foot of typical column, vs. time. Transports are propor- 
tional to deviation of curves from unity. 


1—a, and we should expect them to be less than 
the theoretical values calculated from band 
spectra. In the system CN’—HCN(g) the com- 
puted value of 1—a is 0.026, and the experi- 
mental value from the transport is 0.0126: an 
efficiency of 50 percent. 

Presuming the theory, whose postulates cer- 
tainly seem eminently reasonable, to be suffi- 
ciently confirmed, we obtain interesting light on 
the operation of the columns. Because of the 
slowness of approach to equilibrium, and also 


IS 


30 








Sa aa See eS eS 
10 2 /4 16 18 20 22 24 26 2650 
Doys 





Fic. 7. Fractionation vs. time, experimental. Two units 
working on the system CN~(aq.)—HCN(g). Hutchison, 
Stewart and Urey, reference 6. 


because of the tremendous loss of efficiency with 
production, it is not profitable either to wait for 
equilibrium or to produce near maximum frac- 
tionation (for this would mean producing very 
little). It is preferable to start drawing off as 
soon as the r—¢ curve begins to flatten out. 
One thus reduces the nonproductive time of 
waiting for equilibrium, and, since the column is 
at low fractionation efficiency, obtains a re- 
spectable yield. Finally, the most economical 
way of producing a given concentration of 
isotope is to build a column whose equilibrium 
fractionation is at least 30 percent larger. 

One can obtain from Eqs. (4.24), and (2.8) 
the best conditions for the run. It can be at 
once seen that there is every advantage to 
increasing the concentration of the solution. 
(1—a) varies inversely as the absolute tempera- 
ture T, k varies as 7'*. There is a slight advan- 
tage in fractionation, if it is not counterindicated 
by other considerations, such as polymerization, 
to increasing the temperature. However, the 
period which varies as (1—a)*k, will be also 
increased. 
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Two hypothetical types of structure, derived by folding from basic patterns containing 
extended zigzag polypeptide chains, joined together by hydrogen bridges, are described and 
critically discussed. The (idealized, undistorted) 288-residue open-ended octahedron model has 
suitable horizontal dimensions but is somewhat too large vertically to {it satisfactorily into the 


unit cell of air-dried insulin. 





RINCH! has proposed that, except for 

differences in the R groups, globular pro- 
teins in general have structures derivable from 
the same type of ‘‘fabric,”’ folded, in such a way 
as to maintain stable bond angles and distances, 
to give closed polyhedra. Starting with a given 
type of fabric, she has shown that the possible 
sizes (i.e., numbers of amino acid residues) of 
closed polyhedra of a given class (e.g., truncated 
tetrahedra) are strictly limited. For the ‘‘basic 
cyclol fabric,” for example, she deduces struc- 
tures for molecules containing 72 n? residues, 1 
being an integer. This ingenious postulate thus 
accounts for the definiteness of molecular weight 
for a given globular protein and for the apparent 
“quantumicity” of molecular weights (or of 
numbers of amino acid residues) in globular pro- 
teins as a class. 

The “basic cyclol pattern” and molecular 
models deduced from it have met with a large 
number of objections,” some of which seem to the 
writer to be insurmountable. The idea of folding 
to form a closed system of definite size is, how- 
ever, also applicable to other types of pattern. 
Other types of pattern have, in fact, been sug- 
gested** but models of globular proteins derived 


* Communication No. 762 from the Kodak Research 
Laboratories. Presented before the Division of Colloid 
Chemistry at the Cincinnati Meeting of the American 
Chemical Society, April, 1940. The models discussed in 
this paper were also briefly described in a paper on ‘‘The 
structure of insulin,’ presented before the Division of 
Physical and Inorganic Chemistry at the Baltimore Meet- 
ing, April, 1939. 

1D. M. Wrinch, Proc. Roy. Soc. A161, 505 (1937); 
Trans. Faraday Soc. 33, 1368 (1937). 

2 See, for example, M. L. Huggins, J. Am. Chem. Soc. 
61, 755 (1939) and L. Pauling and C. Niemann, ibid. 61, 
1860 (1939). 

(1936) M. Wrinch and D. Jordan Lloyd, Nature 138, 758 
(1939) T. Astbury and D. M. Wrinch, Nature 139, 798 


from these have not been discussed in print. 

One of these proposed patterns’ consists of 
polypeptide rings held together by means of 
hydrogen bridges (“hydrogen bonds’). From 
chemical evidence and by analogy with other 
known, comparable structures, hydrogen bridg- 
ing would seem to be much more probable than 
cyclol condensation as a means of joining together 
polypeptide chains or rings. 

If globular proteins were composed of small 
rings, containing only a few amino acid residues, 
however, one would expect to find such small 
rings among the decomposition products. So far 
as the writer is aware, they have never been 
found. 

In fibrous proteins, the best available evidence 
indicates structures consisting of polypeptide 
chains (extended or coiled) held to each other in 
sheets or layers by means of hydrogen bridges. 
One is tempted to try structures of this sort as 
basic patterns, to see if, on folding, they might 
produce satisfactory models for globular protein 
molecules. 

In this connection it may be pointed out that 
Astbury and co-workers,* after an x-ray study of 
several globular proteins, state, “The x-ray in- 
terpretation of protein denaturation suggests 
that it always involves the liberation or genera- 
tion of peptide chains which aggregate on coagu- 
lation into parallel bundles like those found in 
the structure of §-keratin and similar fibers.” 
With regard to excelsin, for instance, ‘‘At the 
moment it appears that the original crystal 
structure is symmetrical about a triad axis 
perpendicular to three dyad axes and that it 
degenerates into six sets of fibrous crystallites 


5 W. T. Astbury, S. Dickinson and K. Bailey, Biochem. 
J. 29, 2351 (1935). 
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Fic. 1. Hypothetical protein pattern. 


lying approximately parallel with the basal plane, 
with the side-chain spacing also roughly parallel 
and the backbone spacing roughly perpendicular 
to this plane.”’ 

In this paper will be considered only fabrics 
which are composed of zigzag polypeptide 
chains, fully extended or nearly so. Analogy with 
other known structures (silk, keratin, glycine, 
diketopiperazine, isatin) suggests that these 
chains are held together by means of NHO 
bridges. Tentatively postulating this, we arrive 
at the two patterns represented in Figs. 1 and 2. 
In the pattern of Fig. 1, the sequence —CHR 
—NH—CO- runs in opposite directions in ad- 
jacent chains; in the pattern of Fig. 2, the se- 
sequence is the same in all the chains. 

The first type (Fig. 1) can be folded, without 
breaking any bonds or doing violence to any of 
the accepted ideas regarding bond distances and 
angles, about lines at right angles to the chain 
axes. The angle between adjacent “‘planes’’ of the 
fabric at each fold is 120°, if the bond angles are 
the usual ‘‘tetrakedral” angles. By repeated 
folding, an open-ended hexagonal prism struc- 
ture can be produced (Fig. 3), in which the 
polypeptide chains are now large rings. (Folding 
parallel to the chain axes, to give rodlike prism 
molecules of indefinite length, is also hypo- 
thetically possible. For the present, we shall 
neglect such models.) If we make the reasonable 
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Fic. 2. Hypothetical protein pattern. 


assumptions that all the faces of the basic frame- 
work are alike and that all corners are alike, the 
number of polypeptide rings must be an even 
number and the number of amino acid residues 
in each ring must be a multiple of 12. These re- 
quirements thus limit the total number of 
residues per prism molecule to 24 ”, where is an 
integer. There would be too much crowding in the 
interior of prisms composed of rings containing 
only 12 residues, unless the average interior R 
group is very small, hence one would expect a 
minimum ring size of at least 24 residues, corre- 
sponding to a molecule size of at least 48 residues. 
There is nothing apparent in the nature of this 
model to limit the size further, although one 
might expect a 2-ring structure to be less stable 
than one containing 4 or more rings. A molecule 
containing 288 residues, if built according to this 
scheme, would be expected to consist of 6 rings 
of 48 residues each, or 12 rings of 24 residues each. 

From the bond distances and angles in diketo- 
piperazine and glycine, one would predict an 
average distance per residue, along an extended 
zigzag chain, of about 3.6A. The chain-to-chain 
distance would be expected to be in the neighbor- 
hood of 4A, with this type of structure. The prism 
faces would, therefore, be expected to be about 
29A wide and 20A high (between top and bottom 
chain axes), for one 288-residue model or 14.5A 
wide and 44A high for the other. Neither of 
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Fics. 3 ANpD 4. Illustrating the distribution of poly- 
peptide chain axes in hypothetical molecules derived by 
folding from structures such as those represented in 
Figs. 1 and 2. 


these alternatives seems to be in accord with 
the x-ray data for insulin.® 

The pattern of Fig. 2 folds without strain to 
give open-ended octahedra (Fig. 4). Making the 
requirement that all corners are alike, the num- 
ber of residues per ring must be 12 m and the 
number of rings 6 n. The total number of residues 
per molecule is thus 72 n?. Except for the 72-resi- 


6D. Crowfoot, Proc. Roy. Soc. London, A164, 580 
(1938). 
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due structure, in which there would seem to be 
too much crowding in the interior, the smallest 
possible molecule would have 288 residues in the 
framework. The top and bottom triangular rings 
would be about 29A on a side (the same as the 
octahedron edge in the Wrinch cyclol model). 
Assuming about 4A between chain axes, the 
molecule would have a vertical height (bottom- 
chain axis to top-chain axis) of, roughly, 43A. 
Minor alterations in the structure might decrease 
this figure somewhat, but it appears unlikely 
that it could be small enough to enable molecules 
of this sort to fit into the structure of (air-dried) 
insulin, for which the unit distance in the direc- 
tion of the hexagonal axis is 31A. 

A discussion of structures similar to the fore- 
going but containing kinked or folded chains— 
rather than extended zigzags—will be postponed. 
For the present, the writer merely wishes to 
point out that structures of this general nature, 
which are reasonable from many points of view, . 
have many of the advantages of the cyclol models 
previously proposed, without some, at least, of 
their disadvantages. 
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The Radial Distribution Method in Electron Diffraction 


JoHN WALTER AND J. Y. BEACH 
Frick Chemical Laboratory, Princeton University, Princeton, New Jersey 


(Received April 8, 1940) 


A new radial distribution function for use in interpreting the electron diffraction photographs 
of gas molecules has been obtained. Instead of replacing the inverted Fourier integral by a sum 
the integration is performed by assuming that the shapes of the maxima and minima are those 
of cosine functions. The amount of scattering matter 7*D(r) is calculated rather than the 
density D(r). Comparison between radial distribution functions calculated by this method and by 
other methods is made for a diatomic molecule, a nonlinear triatomic molecule (to test resolv- 
ing power), a linear triatomic molecule and carbon tetrachloride. 





N the study of the structure of gas molecules 
by the electron diffraction method the radial 
distribution method is frequently used to assist 
in the interpretation of the photographs. The 
method consists of taking the theoretical in- 
tensity formula which gives intensity of scatter- 
ing as a function of the distribution of scattering 


ing matter in the molecule as a function of the 
scattering pattern. This method was first de- 
scribed by Pauling and Brockway.' The Pauling 
and Brockway formula has been modified by 
Degard? and Schomaker* to give more reliable 
results. The formulae differ because of the 
different approximations which are used in their 


matter in the molecule and inverting it to obtain derivations. We wish to discuss these approxi- 
a formula giving the radial distribution of scatter- mations and to propose a new formula. 


METHOD 


The complete expression for the coherent scattering by gas molecules is given by the formula‘ 
pate ae (1) 
s4 ST ij 
where s=((47/X) sin 6/2) and — 
A ij =6r,;7/2. 
This expression can be simplified in accordance with the usual visual method by omitting the s* 
in the denominator, the Z— F and temperature factor, however, being retained. We have calculated 
the product (Z;— F;)(Z;— F;)e~4«i* for several interatomic distances in several molecules and have 
found that when s is between 5 and 20 it remains reasonably constant and roughly proportional to 
Z,Z;. As this product, (Z;— F;)(Z;— Fj)e~4:i*, is dependent only on r it can be replaced by r*D(r). 
rD(r) is equal to the product whenever r is equal to 7;; and is equal to zero otherwise. As D(r) is a 
density it must be multiplied by 7? to be proportional to the amount of matter at the distance r. 
The double sum (1) can now be expressed as an integral 


sin sr “ 
dr or slI(s) =K'f rD(r) sin srdr. (2) 
0 


I(s)= Kf r?D(r) 
0 


ST 


This is a Fourier integral and can be inverted to give 
2 


rD(r) =K" f sI(s) sin srds. 
0 


1L. Pauling and L. O. Brockway, J. Am. Chem. Soc. 57, 2684 (1935). 

2 C. Degard, Bull. Soc. Roy. Sci. Liége, 115 (1938). 

’V. F. H. Schomaker, paper presented to the American Chemical Society, Baltimore, April, 1939. 
*L. Pauling and L. O. Brockway, J. Chem. Phys. 2, 867(1934). 


601 








602 J. WALTER AND J. Y. BEACH 
I(s) consists of a series of maxima and minima, the positions and intensities of which are measured 
on the photographs. The shape of a maximum or a minimum can be represented approximately by 
(sp—S) 
I(s) =I(s,) cos ————, 
2a; 
where s; is the measured position of the maximum or minimum, J(s,) is the intensity® of the maximum 
or minimum, and @;, is the distance, in s, from s; to the zero of intensity (minima are given negative 
intensities). a, is roughly the half-width of the maximum or minimum. The integral (3) can now be 








































































































































broken up into a sum of integrals, one for each maximum and minimum, giving 





8kt+ak (s.—S)a 
rD(r)=K" > (sx) Ss cos sin rsds, 
k 8k—@k 2a; 
which is nearly equal to 
sktak (Sx —s)r 
rD(r)=K" > 5,1 (sx) cos — sin rsds, (4) 
k 8k—ak 2a: 


since s is nearly equal to s; in the range s;—a, to s,+a,. This integral can be evaluated in a fairly 


straightforward manner and is equal to 





47a}. 
rD(r)=K" > s,I(s,.)——— cos a;7 sin Sy7 
k a? —4a,?r? 
or 
COS a,r sin S,r 
rD(r)=K'"r? > ay] (sx) Sx (5) 
k w*—4a,?r? Ser 


In applying this expression, the product a;J(s,), which is proportional to the area of the maximum, 


should be estimated directly from the photographs. If the peaks are all of the same width, a, is the 
same for all maxima and minima, and Eg. (5) becomes 


vr? COS Ayr 


r’?D(r)=K"" 


APPLICATIONS AND DISCUSSION 


If the radial distribution function is to repre- 
sent the amount of scattering matter at a dis- 
tance r, the density D(r) must be multiplied by 
4rr® (the 42 has been dropped in Eq. (6)). The 
integral (3) has in previous derivations been set 
equal to a sum with one term for each maximum 
and minimum. Introducing a cosine function for 
the shape of a maximum and performing the 
integration leads to the factor 


a;, cos ayr/(m?—4a,2r*) 


in Eq. (6). In place of this factor Degard? used 
e~°"* where a is equal to our a. This factor, how- 


5 In this treatment J(s) is assumed to be represented by 
Eq. (2). In this equation the intensities of the maxima and 
minima decreases approximately as 1/s. Therefore the 
visually estimated intensities should on the average de- 
crease as 1/s, 


a? —4a;2r? 


The factor cos a,r/(m?—4a,2r*) is relatively insensitive to a; except at large values of r. 


Sin S;P 
Dx Sx2L (Sx) ’ (6) 
5,7 








ever, decreases much too rapidly with increasing 
r. When a=1 it decreases by a factor of 500 
when 7 changes from 1 to 3A°, whereas the factor 
a, cos a,r/(m?—4a,r) decreases by about 40 per- 
cent. e~*” is nearly proportional to 


a, cos a,r/ (2? —4a;2r") 


if a is about one-tenth of ax. 

The factor s,2 in Eq. (6) was originally omitted 
by Pauling and Brockway. Including it increases 
the resolving power considerably, makes the 
peaks sharper and improves their positions 
slightly. Degard? and Schomaker* ® included s,’. 

6 Schomaker’s function, as it has been used in published 
electron diffraction work, is: 

D(r) =D p9spe-2*e? sin ser /sxr, 


where J;° is essentially Js, and e~*% =0.1 for the largest 
ring. 
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Fic. 1. Radial distribution function for a diatomic 
molecule (7 =2.00A). Curve A, ordinary radial distribution 
function, D(r), using nine maxima. Curve B, ordinary 
radial distribution function, D(r), using nine maxima and 
minima. Curve C, modified radial distribution function, 
D(r), using nine maxima and minima. Curve D, radial dis- 
tribution function, r?D(r), present research, using nine 
maxima and eight minima. 


Degard introduced a factor e~*** into Eq. (6). 
This was to take into account the (Z;—F;) 
(Z;— F;) in Eq. (1). However, if the temperature 
factor is also included in Eq. (1), we have shown 
that nothing need be done about the (Z;— F;) 
(2;—F;)." 














Fic. 2. Radial distribution function for a molecule having 
interatomic distances of 2.00A and 2.30A. In each case 
six maxima and five minima were used. Curve A, ordinary 
radial distribution function, D(r) (maxima only). Curve B, 
modified radial distribution function, D(r). Curve C, 
radial distribution function, r?D(r), present research. 





"If it were desired to use the radial distribution method 
to estimate the amplitude of oscillation of an oscillating 
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Fic. 3. Radial distribution function for a linear triatomic 
molecule having interatomic distances 1. 50A, 2.50A and 
4.00A, using 6 maxima and 6 minima. Curve ‘A, ordinary 
radial distribution function, D(r) (maxima only). Curve B, 
modified radial distribution function, D(r). Curve C, 
radial distribution function, r2D(r), present research. 


The radial distribution function (Eq. (5) or 
(6)) should possess peaks whose heights are 
proportional to the average value (with respect 
to s) of n(Z;— F,;)(Z;— F,)e~4ii®* (nis the number 
of interatomic distances at r=r;;). We shall now 
test the method and compare it with the other 
methods. 

We have applied the various radial distribution 
methods to the theoretical intensity curve for a 
diatomic molecule (y= 2.00A). Curve A in Fig. 1 
is calculated by the original method of Pauling 
and Brockway, using maxima only. A bad har- 
monic appears at about 3.8A. If the minima are 
included in the same method, curve B is obtained. 
Curve C is calculated by the method of Scho- 
maker. The peak is much sharper and the base 
line much flatter. Curve D, calculated by Eq. 
(6), is even sharper but the base line is not so flat, 
although it stays closer to the line of zero density 
(the horizontal line). 

To test the resolving powers of the methods 
we have calculated the theoretical curve for a 
X Y2 molecule in which r,, = 2.00A and r,,=2.30A. 
Using these maxima and minima we have calcu- 
lated the radial distribution functions. In Fig. 2 
curve A is calculated by the method of Pauling 
and Brockway (maxima only), curve B by the 
method of Schomaker, and curve C by Eq. (6). 


group the temperature factor should be included to D(r), 


and some such function as Degard’s e~*** introduced to 
F;)(Z;— Fj). 
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Fic. 4. Radial distribution function for CCl4, using ten 
maxima and nine minima. Curve A, ordinary radial dis- 
tribution function, D(r), (maxima only). Curve B, modified 
radial distribution function, D(r). Curve C, radial dis- 
tribution function, r?D(r), present research. 
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Fic. 5. Radial distribution functions, r2D(r), present 
research, for CCly. Curve A, three maxima and two 
minima. Curve B, five maxima and four minima. Curve C, 
seven maxima and six minima. Curve D, ten maxima and 
nine minima. 


The resolution is seen to improve. In curve C the 
resolution is a little too good, the two peaks 
being at 1.98A and 2.32A instead of 2.00A and 
2.30A. This is probably due to the fact that only 
six maxima and five minima were used in the 
calculation. 
Figure 3 shows radial distribution functions 
. calculated from the theoretical intensity curve 
for an unsymmetrical linear triatomic molecule 
having distances of equal importance at 1.50A, 
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Fic. 6. Radial distribution function for isobutane, r?D(r), 
present research. 


2.50A and 4.00A. The peaks on curve A are at 
1.52A, 2.57A and 3.93A; those on curve B are 
1.50A, 2.52A and 2.98A; and those on curve C 
are at 1.51A, 2.49A and 4.00A. 

The theoretical scattering curve was calculated 
for carbon tetrachloride and radial distribution 
curves were calculated from it. These curves are 
shown in Fig. 4. Curve A is slightly different from 
the one calculated by Pauling and Brockway be- 
cause their curve was calculated from actual data 
instead of from the theoretical curve. Curve B 
calculated by the method of Schomaker, is much 
smoother than the curve A. The peaks are at 
1.75A and 2.87A (assumed values were 1.76A and 
2.87A). The peaks on curve C are somewhat 
sharper and are located at 1.76A and 2.87A. The 
ratio of their heights above the zero line, 3.1, is 
nearly equal to the ratio of the average values of 
(Z;—F;)(Z;—F,)e-4*i” for the two distances, 
this ratio being 3.2. In Fig. 5 we have reproduced 
the radial distribution curves for carbon tetra- 
chloride using a smaller number of maxima to 
show roughly how many maxima are needed to 
obtain a certain accuracy. 

In Fig. 6 is shown a radial distribution for 
isobutane. The fairly important C—H peak at 
2.16A which was not present on previously pub- 
lished curves? is now present in the right propor- 
tion to the other peaks. 

We wish to express our appreciation to Profes- 
sor Linus Pauling for valuable criticism and 
advice regarding this paper. 


8 J. Y. Beach and John Walter, J. Chem. Phys., 8, 303 
(1940). 
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Conductivity in Insulators and Its Interpretation* 


A. von HIPPeEL 
Electrical Engineering Department, Massachusetts Institute of Technology, Cambridge, Massachusetts 


(Received April 17, 1940) 


The migration of charge carriers in insulators is discussed in this paper on the basis of a large 
amount of experimental material on hand. Beginning with deviations from Ohm’s law, the 
trapping of charge carriers is considered and the formation of the ‘‘F” band in alkali-halide 
crystals. Then the laws of electronic conductivity in solid dielectrics are formulated. 





HE old conception of insulators as “‘non- 

conductors” of electricity is dead. It 
became undermined by accumulating experi- 
mental evidence demonstrating that conductivity 
can be created in insulators if only the proper 
conditions are chosen. These conditions are 
manifold and changing from case to case; time 
and temperature, field strength and illumination, 
structure and previous history, even the sur- 
rounding atmosphere and the contacting elec- 
trodes prove to be of importance. The definition 
“insulator” has become vague; it can mean any- 
thing but a metal. Consequently theoretical phys- 
ics defined insulators directly as ‘‘nonmetals” 
with the help of the zone structure of the elec- 
tronic levels.! But this is not sufficient ; detailed 
knowledge is needed about the elementary proc- 
esses providing mobile charge carriers and about 
the laws regulating their motion. These laws of 
motion are discussed in the following pages. 


DEVIATIONS FROM OuHM’s LAW 


The concept of conductivity is normally based 
upon Ohm’s law. It implies constant density of 
the charge carriers and assumes that they move 
in the field direction with an average velocity 


v=bE (6= mobility) (1) 
(E= field strength). 


This assumption of a friction factor 1/5 inde- 
pendent of the driving force (1) has proved in- 
valid lately in several cases. 

M. Wien,? measuring the conductivity \ of 


*Invited paper presented before the Division of In- 
dustrial and Engineering Chemistry at the Spring Meeting 
of the American Chemical Society, Cincinnati, 1940. 

‘See, for instance, F. Seitz and R. P. Johnson, J. App. 
Phys. 8, 84, 186 (1937). 

(1928) Wien, Ann. d. Physik 83, 327 (1927); 85, 795 


electrolytes with high impulse voltages (10° 
volt/cm), found an increase of X with field 
strength to a final value Ax, which corresponds 
to the equivalent conductivity for infinite dilu- 
tion. This fact was explained’ with Debye’s con- 
cept of the ionic atmosphere: Each ion in solution 
surrounds itself preferentially with ions of the 
opposite sign, with a space-charge cloud char- 
acterized by a radius and a relaxation time. At 
the high voltages used, the mobility increases be- 
cause the ions begin to move faster than their 
atmosphere can form. 

Another deviation from Ohm’s law was found 
by the author‘ while studying the laws of electric 
breakdown in ionic single crystals. In metals the 
number of free electrons about equals that of the 
atoms, local space-charge effects cannot develop, 
and Fermi statistics result. Consequently the 
electrons have a long free path between inelastic 
collisions because the coupling with the lattice 
ions is weak, the time of interaction is short, and 
the number of free final states into which the 
electrons can go is limited. In nonmetals such as 
the alkali-halides the situation is quite different. 
If surplus electrons exist, their concentration is 
low ; they are slow, because Boltzmann statistics 
apply, and act as space charges distorting the 
surroundings. The electrostatic coupling between 
electron and structure is strong and far-reaching ; 
the electrons are in a state of permanent interac- 
tion with the lattice. The probability of exciting 
lattice vibrations dominates the friction factor 
1/b of Ohm’s law. This probability may be 
roughly described by an excitation function co- 
ordinated to the spectrum of vibrations of the 


3G. Joos and M. Blumentritt, Physik. Zeits. 28, 836 
(1927); G. Joos, ibid. 29, 755 (1928). 

4 A. von Hippel, Zeits. f. Physik 75, 145 (1932); J. App. 
Phys. 8, 815 (1937), 
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material. The interaction time between electron 
and surroundings decreases with increasing veloc- 
ity of the electron; therefore the function passes 
over a maximum above a critical field strength, 
the friction decreases, and acceleration, impact 
ionization, and breakdown result.® 


TRAPPING OF CHARGE CARRIERS 


The two deviations from Ohm’s law discussed 
are interrelated, as Fig. 1 illustrates. Debye’s 
ionic atmosphere can be pictured as a potential 
minimum floating along with the ion under con- 
sideration. It represents a polarization of the 
medium, which is normally slight because the 
number and mobility of the charge carriers of 
opposite sign are about equal. In extreme cases, 
on the other hand, it can become very large, as 
the movement of protons in palladium indicates ;® 
the swarm of electrons around the Ht in the 
metal is so dense that the ion responds only with 
zs of its charge to an outside field. Wien’s effect 
indicates that the potential cup flattens out in 
electrolytes when the ions reach velocities of 
meters per second. 

This final state of ionic migration is, thanks to 
the small electron mass, more or less the initial 
state of electronic motion. The charge carriers 
move along, distorting temporarily the electron 
clouds of the molecules; the nuclei tending to 
swing into new equilibrium positions set up sound 
waves.’ While the heavy ions stay practically at 
the lower end of this friction barrier of vibrations, 


5 The quantitative theory is being developed in the last 
years by H. Froéhlich, R. J. Seeger, and E. Teller, W. 
Franz and F. Seitz, but has not yet reached a final form. 

6 J. Franck, Nach. Ges. der Wiss., Gottingen 44, 293 
(1933); B. Duhm, Zeits. f. Physik 95, 801 (1935). 

7A calculation of the longitudinal vibrations excited is 
given by E. Teller and his co-workers, Phys. Rev. 57. 
i084A (1940). 


the electrons in high fields may pass it (break- 
down effect). 

The space-charge potential around ions in 
electrolytes is the first indication of ionic binding; 
it comes to a full development in lattice struc- 
tures like rocksalt, where each charge carrier is 
surrounded by six opposite ions and kept locally 
trapped in a potential cup several electron volts 
deep. Figure 1 suggests that also around electrons 
a polarization of the medium should take place if 
they are slow enough. In crystals with strong 
binding forces they do not even need to be very 
slow, because the ions respond quickly. Hence 
also electrons should be trapped if they are 
strongly coupled to locally bound particles. 

Such a capture of a surplus electron in a NaCl 
lattice, for instance, may be described quite 
graphically: An electron traveling with the 
kinetic energy of 1/10 e volt has a wave-length 
\=h/mv=3.8 10-7 cm, and acts therefore like a 
space charge of one elementary charge spread 
over a volume }* containing about 2500 lattice 
ions. The distortion produced is small but far- 
reaching, thanks to the slow decay of the Cou- 
lomb potential with distance. The ions respond 
with a top frequency higher than 5.7 10” vibra- 
tions per second (‘‘Reststrahl’’-frequency) ; there- 
fore the distant ions have time enough to be 
somewhat displaced by the electron speeding 
with a velocity of 1.9-107 cm/sec. through the 
structure. A small potential cup begins to form 
around the electron ; it slows down by dissipating 
energy into lattice vibrations. But the Coulomb 
field springing up around the electron, thanks to 
the lattice distortion, increases the total kinetic 
energy of the electron because it must oscillate or 
rotate in this field with a kinetic energy equaling 
half the potential energy.* Therefore the wave- 
length shortens, space-charge density and inter- 
action time increase, the potential cup deepens 
into a local trap, and the electron becomes 
localized around a positive ion of the lattice. 


THE ‘‘F”’ BAND OF THE ALKALI-HALIDES 


If this picture of electron capture is correct, 
an optical absorption band should appear, indi- 
cating by its location the position of the surplus 
electrons in the material and by its color the 


8 N. Bohr, Phil. Mag. 26, 24 (1913). 
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energy needed for their photoelectric liberation. 
Now Pohl and his co-workers? have found an 
absorption band in the alkali-halides, the ‘‘F”’ 
band, which can be created by photo-effect, heat 
treatment in alkali vapor, or electric current at 
high temperature. This band moves in the electric 
field like an electron cloud, as Stasiw”® first found, 
and has been identified by the author with the 
trapping band he expected." But there is a 
strong controversy at the moment about the cor- 
rectness of this model; another possibility exists, 
which we shall consider next, because the discus- 
sion of it introduces concepts of importance. 
Landau” first proposed a quantum-mechanical 
trapping mechanism for electrons, taking the 
standpoint that slow electrons cannot be trapped 
on account of their long wave-length, but that 
electrons might be captured if speeded up suffi- 
ciently by an activation energy. The author 
came, without knowing of this paper, to the op- 
posite conclusion—that slow electrons will be 
trapped by the mechanism outlined above. By 
comparing both processes it can be seen that the 
Coulomb field of the distorted surroundings per- 
forms the service of the activation energy to 
shrink the wave-length of the electron. Gurney 
and Mott" took first our standpoint—that an 
electron can be trapped in the ideal lattice—but 
later they discarded this model" in favor of one 


Distance from missing ion———> 








Fic. 2. Potential energy of an electron in the neighbor- 
hood of a missing negative ion according to de Boer and 
Mott. The dotted lines give the curve —e?/xv. 


ass for instance, R. W. Pohl, Proc. Phys. Soc. 49, 3 
10 QO. Stasiw, Nach. Ges. Wiss., Géttingen, No. 50 (1933). 
A. von Hippel, Ergebnisse d. exakt. Naturwiss. 14, 

113 (1935); Zeits. f. Physik 101, 680 (1936). 

2 L. Landau, Physik. Zeits. Sowjetunion 3, 664 (1933). 
1 R. W. Gurney and N. F. Mott, Proc. Phys. Soc. 49, 

32 (1937). 

“That the reason given at the Bristol Conference 

(Proc. Phys. Soc. 49, 36 (1937)) is invalid, has been pointed 

out already (A. von Hippel, J. App. Phys. 8, 832 (1937)). 
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given by de Boer.'® De Boer assumes that an 
electron is captured at a lattice point where a 
negative ion is missing (Fig. 2). That such vacant 
lattice points exist can be concluded from the fact 
of ionic conduction. The ions in an ideal lattice, 
like rocksalt, fill out the crystal structure by 
very tight packing; their activation energy is 
much too high to explain the large conductivity 
observed. There must be defects in the structure, 
such as ions in interstitial positions or transferred 
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Fic. 3. Lattice disorder. 


to the surface, leaving vacant places (Fig. 3) for 
ionic migration.'* In NaCl, in thermal equi- 
librium at 1000°K, about 10-° of the lattice 
points should be disordered to account for the 
ionic mobility measured. Frenkel, Wagner, Jost, 
and Schottky mainly have developed this 
theory, and Schottky and others!’ have calcu- 
lated that in the alkali-halides the dominant kind 
of disorder is his “Type 4” of vacant places 
(see Fig. 3). 

It is of importance to establish which of the 
two models given above represents the “‘F’’ band, 
because a deeper insight into the properties of 
ideal and defective crystal structures is involved. 
From the standpoint of our model, a decision 
might be reached by the following cycle process 
(Fig. 4): Given the ideal lattice of NaCl and a 
migrating surplus electron in its conduction 
band. Now (1) remove one Na? ion into vacuum 
by spending the lattice energy Up=183 kcal. 
minus a polarization energy U,, because the lat- 
tice around the hole becomes distorted and settles 
down into a new equilibrium position. U,, ac- 


9st) H. de Boer, Rec. Trav. Chim. Pays-Bas 56, 301 

16 J. Frenkel, Zeits. f. Physik 35, 652 (1926); C. Wagner 
and W. Schottky, Zeits. f. physik. Chemie 11, 163 (1930); 
W. Jost, J. Chem. Phys. 1, 466 (1933); W. Schottky, 
Zeits. f. physik. Chemie 29, 335 (1935). 

17 W. Jost and G. Nehlep, Zeits. f. physik. Chemie 32, 1 
(1936); N. F. Mott and M. J. Littleton, Trans. Faraday 
Soc. 34, 485 (1938). 
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Fic. 4. Cycle for calculation of release energy. (F band.) 


cording to Mott and Littleton,!’ amounts to 
about 77 kcal. (2) Remove the electron from the 
crystal by spending the energy of the work func- 
tion x=11.5 kcal.!® (3) Combine electron and 
ion, gaining the ionization potential J=118 kcal. 
(4) Place the sodium atom back into the hole; a 
small energy gain Q in the order of 37 kcal. re- 
sults.!§ (5) Remove the electron by light absorp- 
tion back into the conduction band ; the quantum 
hy, necessary should represent the energy of the 
“F” band. Result: 


hye= — Unt U,-—xt+14+2=37.5 kcal. 
=1.63ev (2) 


as compared with the experimental value 2.73 
ev.!® This difference can be in the limits of error, 
and the value calculated indicates that an elec- 
tron will be trapped in the ideal lattice, but a 
decision between the two models can only be 
made by new experimental evidence. We are 
engaged in this study. 


ELECTRONIC CONDUCTIVITY IN SOLIDS 


Conductivity in solid insulators, as the preced- 
ing discussion shows, is normally not produced 
by a free flow of charge carriers through a resist- 
ant medium but by a progression of the charges 
in jerks. The carriers are trapped after an average 
displacement distance #; they stick until the 
heat vibration of the structure supply with 
statistical probability the activation energy for 
their release. 


18N.F. Mott, Trans. Faraday Soc. 34, 500 (1938). 
19 E. Mollwo, Zeits. f. Physik 85, 56 (1933). 


HIPPEL 


For electrons in crystals this conduction 
process can be*followed up step by step in direct 
experiments. The liberation of the charges can be 
produced photoelectrically by a narrow light 
beam; if the insulator is kept cool, no second re- 
lease of the electrons after capture will take place. 
The current measured as function of distance 
between light beam and anode (Fig. 5) reflects 
directly the decay of the number of charge car- 
riers with the distance of migration.?° An expo- 


Photo Current 


—Width of light beam 
in 10° amp. 


Light Beam 


‘ibe 





LULLULLL 











ilF-|I] 


X, distance from anode to light beam, in mm. 


Fic. 5. Measurement of the displacement distance of 
electrons in AgCl at —186°C. 


nential decay law has been found, 
N= Noe=*!”; (3) 


the displacement distance wW increases propor- 
tionally to the field strength applied and has a 
value Wp, of 2.5-10-5 cm in AgCl, and of about 
10-7 cm in NaCl at a field of 1 volt/cm.”! The 
large difference between the two crystals of 
identical lattice structure is to be expected from 
the picture of electronic migration given in the 
preceding paragraphs. While the alkali halides 
are strongly ionic crystals coupling the electron 
intensely to an extended volume of the material, 
the silver halides have a large portion of atomic 
binding, and a much weaker and more limited 
interaction results. Values for the free path of 
surplus electrons found by ‘‘Hall-effect’’ meas- 
urements give additional evidence. No effect 
could be found in NaCl ” indicating that the free 
path was below 10-7 cm, while in diamond 
crystals a large free path of 10-* cm has been 
observed.” It is quite possible that in atomic 
structures like diamond the coupling is so weak 


20 W. Flechsig, Physik. Zeits. 32, 843 (1931); K. Hecht, 
Zeits. f. Physik 77, 235 (1932). 

21 G, Glaser, Nach. Ges. der Wiss. Gottingen, 3, No. 2 
(1937). 

2]. Evans, Phys. Rev. 54, 47 (1940). 

23H. Lenz, Ann. d. Physik 77, 449 (1925); 82, 
(1927). 
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CONDUCTIVITY 


that electrons are not trapped in the free volume 
but only on boundaries or impurities; such 
materials should also have a comparatively low 
breakdown strength if they do not possess elec- 
tronic excitation levels as a second barrier 
against impact ionization. 

The velocity of the surplus electorns in alkali- 
halide crystals can be measured by observing the 
motion of the color centers forming the ‘‘F’’ band 
or by calculating the density of the centers and 
the electronic current transferred.** An exponen- 
tial temperature dependence 


V=V9-e-U/kT (4) 


has been found in this way which can be easily 
understood. The outside observer, integrating 
over the intermittent motion of the charge car- 
riers records the average velocity 


v=W/ (titty), (5) 


where ¢; marks the mean time of capture at one 
point, and ¢, the mean time of free travel between 
two capture acts. The time ?¢; is ended statisti- 
cally by the temperature motion of the crystal 
supplying the activation energy U for release of 
the electron 


tp=a-eU!/*?, (6) 
t; is normally <t,; therefore the form (4) results 
v=(w/a)-e-U/*T (7) 

and the conductivity can be written 


h= N-e-(Wo/a)-e-U!/*T 
(N=number of surplus electrons per cm*). (8) 


The factor a measures the degree of coupling 
between the trapped electron and its surround- 
ings; 1/a can be interpreted as the frequency of 
oscillation produced by this binding. It should be 
appreciably lower than the ‘Reststrahl’’ fre- 
quency of the lattice points because the size of 
an alkali atom is large; therefore the electron will 
be shared in the trapped state by one central ion 
and a group of neighbors. With Smakula’s 
values for NaCl:?* U=0.94 ev, v9 =20 cm/sec./ 
volt/em, T7=973°K; and Glaser’s”! w29=1.10-% 


(1934) Smakula, Nach. Ges. der Wiss. Géttingen, 1, No. 4 
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cm*5) for a concentration of color centers 10!7/cm? 
we find: a=5-10- sec., ¢;=4.9-10-° sec.; 
vp=1/a=2-10° sec.—! as expected Kv,=5.7-10" 
sec.—'. With increasing temperature the lattice 
structure widens, the oscillation frequency goes 


- down, and the maximum of the ‘‘F”’ band shifts 


towards the red, as observed.”* The ratio between 
thermal activation energy U and optical absorp- 
tion energy hy; is 1 : 2 in the limits of error; this 
it should be, according to the Franck-Condon 
principle, if the electron is trapped by lattice 
distortion. 

Equation (8) takes care of the conductivity 
produced by surplus electrons, but this is only 
one side of the problem. If these movable 
charge carriers are created by transferring elec- 
trons froui a lower filled zone into a conduction 
band, normally two possibilities for conduction 
are created in each case—one excess electron 
and one hole (Fig. 6). By treating the hole as 
positron of electronic mobility, the anomalous 
Hall effect has been explained,?’ and since that 
time the theoretical discussion has normally dis- 
tinguished between the two modes of conduction 
only by the sign. But the processes may be very 
different, as the case of the alkali-halides demon- 
strates. If an electron is transferred from the halo- 
gen ion into the conducting state, it represents the 
excess electron moving along distances w between 
trapping, as discussed; the observer sees a 
sodium atom moving relatively fast towards the 
anode. The missing electron at the halogen ion 
represents the hole; it may also move and trans- 
fer the chlorine atom towards the cathode, but 
it can do so only in jerks ®, of atomic distances 
by electron exchange from neighbor to neighbor. 
Keeping this difference in mind, the hole conduc- 
tivity can be represented by an equation like (8) 


bar) Atom 
BAND PICTURE OF 
ELECTRON AND HOLE 








ELECTRON AND HOLE 
MIGRATION IN NoCI 


Fic. 6. 


It is hard to see how such small distances can be 
explained with the de Boer-Mott model of trapping. 

26 U and v in (4) are therefore dependent on T. 

27 W. Heisenberg, Ann. d. Physik 10, 888 (1931). 
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and semiconductors with excess- and defect- 
conductance can now be treated as individuals 
characterized by N, w, a, and U. 


CONCLUDING REMARKS 


The discussion of conductivity in the preceding 
paragraphs has limited itself to the laws of mo- 
tion. But charges have not only to be transported, 
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they have to be generated; they have to enter 
the material and to leave it; they have to be 
balanced by counter-charges providing the neces- 
sary electroneutrality. Only this more general 
picture, which includes the space-charge effects 
and field distortion, can explain the peculiar 
reactions of insulators observed in rectifiers and 
barrier-layer photo-cells. The author will come 
back to these questions elsewhere. 
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The present paper describes a continuation of the work on the heat capacities of lower hydro- 
carbons by the Lummer-Pringsheim adiabatic expansion method. The apparatus and the pro- 
cedure described in the first paper of this series, henceforth to be denoted as Part I, was used in 
the present research without important modifications, 


INCE an absolute calibration of the resistance 

thermometer forms the most important 
feature of the present work, Fig. 1 shows the 
resistance temperature coefficients of the Wol- 
laston wire which was used in most of the meas- 
urements here described. To appreciate the re- 
liability of these coefficients it may be pointed 
out that in a year’s time the resistance of the 
wire at the ice point has changed by only 0.02 
ohm out of a total of 280 ohms. Determinations 
of the thermal resistance coefficients which were 
made with every gas studied gave points which 
fell perfectly on the curve given in Fig. 1. 

As discussed in Part I the temperature of the 
wire does not remain constant after expansion 
but changes linearly with time. This was tenta- 
tively explained as due to absorption of thermal 
radiation from the warmer walls of the expansion 
vessel by the gas and it was proposed to mini- 
mize the effect by using a vessel with surfaces of 
very low emissivity. With this idea in mind the 
present research was conducted with the same 
expansion vessel and thermometer wires as in 
Part I except that the inside of the expansion 


1 Part I, G. B. Kistiakowsky and W. W. Rice, J. Chem. 
Phys. 7, 281 (1939). 


vessel was gold-plated and polished to a mirror- 
like surface. Unfortunately it was found that the 
constancy of the wire temperature after expan- 
sion was thereby only slightly improved. 

The failure of the gold plating to improve the 
appearance of the experimental records and to 
reduce experimental errors is very disappointing. 
However, this failure does not eliminate radiation 
as the possible source of the temperature in- 
constancy because at the low temperatures of 
the experiments rather long wave-lengths may 
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Fic. 1. The resistance-temperature coefficient of the 
Wollaston wire thermometer plotted against the tem- 
perature in °C, 
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be chiefly responsible for the heating of the gas 
and for these wave-lengths the ordinary un- 
polished brass surface is hardly a poorer reflector 
than polished gold plate. Dr. Bernard Lewis has 
called to our attention the possibility that the 
changing temperature after expansion may be 
the evidence that only incomplete energy ex- 
change occurs between the translational and some 
internal degrees of freedom of the molecules. In 
that case extrapolation to zero time would give 
erroneous values of the heat capacities and hence 
the whole method would be in doubt. However, 
the perfect agreement of thus measured values 
with the theory in the case of carbon dioxide and 
that of ethane (cf. Part I) practically eliminates 
this possibility and we consider therefore the dis- 
cussion of errors as given in Part I to apply to 
the present data also. 

An additional error which was insignificant in 
experiments on ethane makes its presence felt in 
the present research. It is the correction of the 
experimental data to the ideal gas state? which 
is so uncertain with some of the gases here dis- 
cussed that the resultant error is larger than all 
the other errors. 

For some of the hydrocarbons, fortunately, 
accurate equation of state data have recently 


TABLE I. Correction to ideal state. All values in 
cal./mole degree. 








Cp®—C,’ Cp—Cy’ 


ComMPOUND Catc.1 Cartc.2 Catrc.3/Catc.1 Catc. 2 





Acetylene : 07 
Ethylene d .06 .06 
Ethane ; AZ 
Propylene ‘ 30 
Allene , 33 
Trans-butene-2 | . 78 
Ethylene i 05 05 
Ethane d .08 
Propylene ; 18 
Allene BF .20 
370 ~Trans-butene-2 | .: A7 














Calc. 1 from the known equation of state data. 

Calc. 2 from the known critical data and the Keyes equation of state. 

Cale. 3 from the known critical temperature, graphically estimated 
critical pressure and the Keyes equation of state. 


? The recalculations of the second virial coefficient of 
ethane by E. E. Roper (J. Chem. Phys. 8, 290 (1940)) from 
the data of Eucken which was used in Part I, changes the 
heat capacity data there given by less than 0.01 cal./mole 
degree. The value of the gas constant (R) used in all these 
calculations was taken as 1.9864 cal./mole degree. The 
i en temperature of the ice point was taken as 
4/5. ‘ 
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Fic. 2. A graph of the critical pressure vs. the critical 
temperature for the C2, C3 and C, hydrocarbons. Circles 
represent the saturated hydrocarbons; circles with crosses 
represent the unsaturated hydrocarbons. 


been obtained by Lamb and Roper’ and they 
form, so to speak, landmarks by which others 
can be estimated. The general mathematical pro- 
cedure for the correction to the ideal state was 
the same as described in Part I. The gas correc- 
tion for the hydrocarbons allene, propylene and 
both butene-2 isomers reported here were ob- 
tained from Roper’s measurements. For propane 
critical data are available and these were substi- 
tuted into the Keyes equation of state* thus 
giving the second virial coefficient. The gas cor- 
rection for propane was also calculated using the 
Beattie equation of state® and the correction was 
found to agree within 0.02 cal./mole degree with 
that obtained by the use of the critical data. For 
propane the critical data used were T, = 368.7°K 
and P,=43 atmos. For methyl acetylene only the 
critical temperature has been determined and is 
given as 394.8°K.® We have collected therefore, 
the known critical data for lower hydrocarbons 
which show that a simple relation exists between 
the critical temperature and pressure of these 
compounds (with the exception of acetylene) as 
seen in Fig. 2. From this graph the critical pres- 
sure of methyl acetylene is obtained as 40.5 at- 
mos. if it is assumed that this compound obeys 


3 A.B. Lamb and E. E. Roper, J. Phys. Chem. (to appear 
shortly). 

4F, G. Keyes, J. Am. Chem. Soc. 60, 1761 (1938). 

5 Beattie, Kay and Kaminsky, J. Am. Chem. Soc. 59, 
1589 (1939). 

6 F. R. Morehouse, Can. J. Research 5, 306 (1931). 
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TABLE II. Heat capacity of allene. 








dR OHM 


Pinal FT sea 
DEG. 


AV. MM 


RANGE OF 
AP MM 


Tw 
°K 


C,° 
CAL. 


Rap. Corr. a 
DEG. MOLE“! 


PERCENT Cy’ 





272.16 
300.00 
334.00 
366.45 


40.230-29.048 
33.614-16.081 
33.731-18.100 
34.752-16.649 


759.65 
754.70 
752.20 
754.80 


1.2245 
1.1110 
1.1004 
1.0887 

















13.20 
14.18 
15.25 
16.22 


12.93 +.03 
13.96+.02 
15.06+.03 
16.08+.02 




















* The wire used in these measurements was annealed less extensively. 


TABLE III. Heat capacity of methyl acetylene. 








Pfinal 
AV. MM 


RANGE OF 
AP mM 


dR OHM 
dT DEG. 


Cp® 
CAL. 


Rap. Corr. 
DEG. MOLE~! 


PERCENT C,’ 








37.094—28.454 
37.080-28.274 
34.923-31.447 
35.555-24.721 








768.10 
772.10 
770.80 
768.30 





1.1182 
1.1110 
1.1010 
1.0876 





13.42+.01 
14.29+.02 
15.31+.01 
16.31+.03 





13.76 
14.58 
15.56 
16.52 











TABLE IV. Heat capacity of cyclopropane. 








Tay 
°K 


RANGE OF 
AP MM 


Pinal 
AV. MM 


dR OHM 
dT DEG. 


PERCENT 


Rap. Corr. 


Cy’ 


Cp® 
CAL. 
DEG. MOLE™! 





272.15 
300.48 
333.70 
368.46 





36.754-28.622 
37.785-28.231 
37.710-26.456 
39.244-25.083 








762.80 
766.00 
767.10 
769.20 





1.1182 
1.1108 
1.1006 
1.0878 


0.13 
12 
14 
.20 





11.83+.02 
13.28+.02 
14.91+.03 
16.62 +.03 








12.10 
13.50 
15.10 
16.77 











the normal relation. These critical values for 
methyl acetylene were used in the Keyes equa- 
tion to obtain the gas correction for this com- 
pound. For cyclopropane, finally, no critical data 
are available; since the boiling point of this sub- 
stance is very close to that of allene, it was as- 
sumed that the second virial coefficients are 
identical. 

Table I gives the results of calculations of the 
correction to the ideal state by the different 
methods here outlined, insofar as enough data 
were available to perform more than one calcula- 
tion on the same substance. The data of the first 
column were obtained from reliable equation of 
state data and should be accurate to 0.02 to 0.04 
cal./mole degree depending on the magnitude of 
the correction. For acetylene, measurements of 
Schaefer’ were available: for ethylene, measure- 
ments of Roper and for ethane, measurements 
of Eucken.® 

The agreement of the first column in Table I 
~ 1K, Schaefer, Zeits. f. physik. Chemie B36, 85 (1937). . 


8 Eucken and Parts, Zeits. f. physik. Chemie B20, 184 
(1933). 


with the other two shows that the errors due to 
extrapolation to the ideal state are about 0.05 
to 0.08 cal./mole degree for gases for which no 
accurate equations of state are available. For 
cyclopropane they may rise to as much as 0.1 
cal./mole degree because of the complete lack of 
experimental data. These errors are estimated to 
apply only to the correction from C,’ to C,°. 
When the Keyes equation is employed to deter- 
mine C, (the heat capacity at a pressure of ap- 
proximately one atmosphere) very much larger 
errors are committed. This is demonstrated in 
Table I for the case of propylene and it is found 
that the Keyes equation gives too low values for 
both dB/dT and d?B/dT?, so that the calculation 
of C, results in too low values. This is especially 
true at temperatures near the boiling point of the 
compound. Fortunately, however, in the further 
correction to C,° the first error is almost com- 
pletely compensated. 

In the following paragraphs we shall take up, 
one by one, the five hydrocarbons with which 
this paper is concerned. 
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ALLENE (C3H,) 


The sample of allene used in these heat capac- 
ity determinations was part of the material used 
in the measurements of the heat of hydrogenation 
by Kistiakowsky and co-workers.’ It was con- 
taminated by air and moisture and was therefore 
subjected to repeated bulb to bulb distillations 
in vacuum, pumping after each distillation, and 
was finally distilled into the expansion vessel and 
the storage flasks. 

Table II presents the summary of the experi- 
ments on allene. None of the experiments have 
been disregarded because of too large deviations 
from the mean in forming the averages of this 
table. The same is true of other heat capacity 
tables given in this paper. 

In this table and those following, 7 repre- 
sents the average of the temperature of the gas 
before and after expansion, AP gives the pressure 
change of the gas upon expansion, Pyinai repre- 
sents the average barometric pressure for the 
series of runs, dR/dT is the temperature coeffi- 
cient of resistance for the thermometer wire. 
The column labeled ‘‘Radiation Correction” is 
discussed in Part I. The next column gives C,’, 
the uncorrected heat capacity as given by Eq. 
(3b) in Part I and the mean deviation from the 
average. C, gives the heat capacity of the gas at 
atmospheric pressure, see Eq. (3a) in Part I, 
and C,° represents the heat capacity of the gas 
at infinitely low pressure. Finally the last column 
of the table gives the theoretical heat capacities 
calculated using the harmonic oscillator rigid 
rotator model with the frequency assignment 
which has been given by Linnett and Avery :!° 
2992(1), 1430(1), 1069(1), 820(1), 3000(1), 
1968(1), 1389(1), 3062(2), 1031(2), 852(2), 
353(2), all frequencies in cm~! and the numbers 
in parenthesis indicating the statistical weights 
attached to each frequency. 

The agreement between the experimental and 
theoretical values is excellent. Only at the lowest 
temperature is the difference slightly larger than 
the mean deviation of the individual experi- 
mental values from the average. The experiments 
at this temperature however were made much 


_* Kistiakowsky, Ruhoff, Smith and Vaughan, J. Am. 
Chem. Soc. 58, 146 (1936). 

10]. W. Linnett and W. H. Avery, J. Chem. Phys. 6, 
686 (1938). 
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earlier than the rest of the experiments, using a 
wire which was not so thoroughly annealed. This 
same wire gave less satisfactory results also in 
the work on air (see Part I). The agreement be- 
tween the experimental and the theoretical heat 
capacities may be regarded as an indication of 
the correctness of the above vibrational assign- 
ment. In fact, the agreement of experimental and 
theoretical values in Table I is somewhat better 
than we should expect allowing for the accidental 
and systematic experimental errors and the un- 
certainty of extrapolation to the ideal gas state. 
Together these errors may amount to as much as 
0.4 percent. On the other hand the molecular 
model used in the theoretical calculations is 
probably too simple but the effect of this on the 
heat capacity we cannot estimate. It may be only 
suggested that all these errors cancel more or less 
accidently. 


METHYL ACETYLENE (C3H,) 


A sample of methyl acetylene was available 
from the work on its heat of hydrogenation" and 
it was used without purification except for an 
ordinary bulb to bulb distillation in vacuum. 

Table III presents a summary of the experi- 
mental data for methyl acetylene. The theoretical 
heat capacities have been calculated using the 
harmonic oscillator, rigid rotator model and the 
vibrational frequencies proposed by Crawford” 
who suggests the following: 926(1), 1380(1), 
2125(1), 2910(1), 3300(1), 333(2), 642(2), 1035(2), 
1444(2), 2975(2). The frequencies were ob- 
TABLE V. Corrected values for the free energy, entropy and 

heat capacity of cyclopropane. 








T°K _ —(F -Eo)/T Ss Cy? 


100 38.71 46.66 7.96 
150 41.93 49.92 8.25 
200 44.25 52.43 9.30 
250 46.11 54.69 11.14 
300 47.73 56.92 13.44 
350 49.21 59.18 15.88 
400 50.60 61.46 18.33 
500 53.22 66.03 22.64 
600 55.72 70.47 26.14 
700 58.15 74.73 29.04 
800 60.47 78.78 31.45 
900 62.71 82.61 33.57 
1000 64.90 86.24 35.37 




















1 Conn, Kistiakowsky and Smith, J. Am. Chem. Soc. 
61, 1868 (1939). 
2B. L. Crawford, J. Chem. Phys. 7, 140 (1939). 
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TABLE VI. Heat capacity of propylene. 








Pinal 
AV. MM 


RANGE OF 
AP mM 


dR OHM 


dT DEG. 


Cp® 
CAL. 


Rap. Corr. _— 
DEG. MOLE! 


PERCENT C,’ 








765.40 
765.30 
769.75 
772.30 


37.937-32.379 
37.798-23.211 
36.182-22.441 
35.656-19.363 














1.1182 
1.1110 
1.1006 
1.0884 


14.06+.03 
15.22+.03 
16.54+.04 
17.76+.02 


14.36 
15.47 
16.74 
17.93 




















TABLE VII. Calculated heat capacity of propylene. 








Cvib Cir. 
ASSIGNMENT | ASSIGNMENT 
C,?® USED BY USED BY 
Exp. me An Re K. L. R. BY P 


Cyib 


ASSIGNMENT 


Cir, 
ASSIGNMENT 
, r 


AND W. 


Cvib 
Cir. ASSIGNMENT 

ASSIGNMENT BY W 
BY P. AND W. 


Cir. 
V =2100 
CAL./MOLE 





14.36+.03| 5.268 1.15 5.021 
15.47+.03} 6.399 1.12 6.115 
16.74+.04| 7.882 91 7.566 
17.93+4.02| 9.279 .70 8.924 

















4.331 
5.337 
6.679 
7.995 


2.08 
2.19 
2.11 
1.99 


2.13 
2.11 
1.98 
1.89 




















tained by Crawford from the Raman spectrum 
of liquid methyl acetylene and hence it is not 
excluded that the true frequencies of the gaseous 
molecule are slightly higher than the above. How- 
ever the resultant change in the theoretical heat 
capacity should be quite small and will not alter 
materially the excellent agreement of the experi- 
mental and theoretical values given in Table ITI. 
The present work therefore provides strong sup- 
port for Crawford’s analysis of the vibrational 
spectrum of methyl acetylene. 


CYCLOPROPANE (C3Hg) 


The sample of cyclopropane was obtained 
through the courtesy of Dr. V. H. Wallingford 
of the Mallinckrodt Chemical Company to whom 
we wish to express our indebtedness. It was part 
of the material which was analyzed by Dr. F. D. 
Rossini of the National Bureau of Standards who 
found 0.25 mole percent of some impurity." The 
sample was freed from air and moisture by bulb 
to bulb distillations in vacuum so that the re- 
maining impurity should be composed of sub- 
stances with heat capacities not greatly different 
from that of cyclopropane. The uncertainty in 
the experimental heat capacity values resulting 
from this effect is certainly less than 0.25 percent. 
Table IV gives the experimental heat capacities 
obtained for cyclopropane. 

Linnett" has studied the Raman and infra-red 


18 Private communication from Dr. Wallingford. 
4 J. W. Linnett, J. Chem. Phys. 6, 692 (1938). 


spectrum of cyclopropane and was able to give a 
frequency assignment to all normal vibrations 
but one. Using the heat capacity data of the hot 
wire method, then available in this laboratory, 
he concluded that the unknown frequency is 
1250 cm. This value gave satisfactory theoret- 
ical heat capacities and was used to explain the 
observed 385 cm@! frequency as a difference 
(1250-865). 

Theoretical heat capacities calculated on Lin- 
nett’s original assignment for cyclopropane 
deviate from the present experimental data by 
approximately 0.1 to 0.2 cal./mole degree, at the 
lowest and the highest temperatures, respectively. 
The theoretical heat capacities presented in 
Table IV have been calculated therefore with the 
following assignment: 3000(1), 1485(1), 1187(1), 
1070(1), 3000(2), 1435(2), 1022(2), 860(2) 
1000(1), 3050(1), 1041(1), 3080(2), 1125(2), 
740(2), which differs from that of Linnett in that 
1125 cm is taken as the frequency of the un- 
known vibration. They are seen to be in excellent 
agreement with the experimental data. The new 
choice of the unknown frequency has the further 
advantage that the observed 385 cm frequency 
can be represented now as a transition from the 
lowest vibrating state of the molecule, or (1125 
—740= 385). This transition is permitted by 
the selection rules as an inspection of Linnett’s 
paper will readily show. 

The values of the thermodynamic functions of 
cyclopropane as calculated by Linnett are 
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changed only very slightly by the suggested 
modification of the fundamental frequencies, cf. 
Table V. 


PROPYLENE (C3Hg) 


The propylene used in the present measure- 
ments was part of the material obtained for the 
determination of its heat of bromination,'® and 
was stated to be of 99.9 percent purity. It was 
contained in the original cylinder but like the 
other hydrocarbons discussed in this paper was 
subjected to repeated bulb to bulb distillations 
in vacuum, to remove traces of air and water be- 
fore using. Table VI presents the summary of the 
experimental heat capacities obtained for 
propylene. 

In a preliminary note Kistiakowsky, Lacher 
and Ransom!'® have discussed the heat capacity 
data of propylene at low temperatures obtained 
by the hot wire method. They concluded that 
the experimental data was in agreement with the 
suggestion of Pitzer!’ that the internal rotation of 
the methyl group in propylene is hindered by a 
potential barrier of 600 to 800 cal./mole, there 
being three positions of minimum potential 
energy. 

As a means of comparison we have compiled 
Table VII which gives in the second column a 
summary of the new experimental data on the 
heat capacity of propylene. Column three records 
the vibrational contribution to the heat capacity 
at the various temperatures according to the 
assignment used by Kistiakowsky, Lacher and 
Ransom. Column four gives the residues obtained 
on deducting the vibrational heat capacity plus 
4R from the experimental values for the heat 
capacity. This should be the heat capacity due to 
internal rotation, but since the values decrease 
considerably below the classical value of (43)R 
for free rotation at the higher temperatures it is 
evident that the vibrational assignment is er- 
roneous. The next two columns contain the vi- 
brational heat capacity and the contribution due 
to internal rotation calculated from an approxi- 
mate vibrational assignment given by Pitzer,!’ 
and the following two columns similar calcula- 


® Conn, Kistiakowsky and Smith, J. Am. Chem. Soc. 
60, 2766 (1938). 

16 Kistiakowsky, Lacher and Ransom, J. Chem. Phys. 6, 
900 (1938). 

“K. S. Pitzer, J. Chem. Phys. 5, 473 (1937). 
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tions using the recently developed vibrational 
assignment of Wilson and Wells,'* based on an 
assignment which is still only provisional, al- 
though based upon an experimental study of the 
infra-red and Raman spectra of propylene and a 
normal coordinate treatment of the molecule. 
This assignment is given as: 432(1), 580(1), 
921(1), 940(1), 1000(1), 1022(1), 1100(1), 1111(1), 
1297(1), 1399(1), 1445(1), 1468(2), 1647(1), 
3000(6), cm-'. The last two columns of the 
Table VII contain internal rotation contributions 
calculated from Pitzer’s tables!® for the potential 
barrier of 800 cal./mole and Crawford’s method*° 
for the higher barrier of 2100 cal./mole. 

On the basis of Table VII alone a final decision 
between Pitzer’s assignment and a low potential 
barrier or the Wilson-Wells assignment and a 
high barrier would be difficult, although the 
differences between experimental and calculated 
data can be better smoothed out by a small 
change in the height of the barrier in the latter 
case. The better agreement between the Wilson- 
Wells calculation and experiment is more pro- 
nounced at lower temperatures, where the hot 
wire method has been used. Finally, the third 
law entropy data on propylene, obtained re- 
cently by Powell and Giauque,”! favors strongly 
the higher barrier. 

These authors themselves draw the conclusion 
that the third law is inapplicable to propylene 
because their experimental entropy is lower by 
about 1 E. U. than that calculated using the low 
potential barrier and either the Pitzer or the K.., 
L., and R. vibrational assignment. However, with 
the Wilson-Wells assignment and the 2100 cal./ 


TABLE VIII. Theoretical and experimental entropy 
of propylene. 








225.35°K 298.1°K 
S(t.4r.) 57.00 59.23 


S (vib) .89 2.00 
Sc. r.) (V=2100 cal./mole) 1.98 


S (total calc.) 59.87 
S (third Law) 59.93 











* S2s.1 —S225.3 obtained by graphical integration of the experimental 
heat capacity curve. 


18 Crawford, Kistiakowsky, Rice, Wells and Wilson, 
J. Am. Chem. Soc. 61, 2980 (1939). 

19K. S. Pitzer, J. Chem. Phys. 5, 469 (1937). 

20 B. L. Crawford, J. Chem. Phys. 8, 273 (1940). 

*1T. M. Powell and W. F. Giauque, J. Am. Chem. Soc. 
61, 2366 (1939). 
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TABLE IX. Heat capacity of propane. 








Pfnal 
AV. MM 


RANGE OF 
AP MM 


Cp° 
Rap. Corr. CAL. 
PERCENT C,’ - DEG. MOLE™! 





765.36 
761.80 
753.85 
763.05 


37.954-30.245 
41.380-30.669 
41.878-32.450 
40.367-31.583 














15.84+.01 16.19 
17.28+.02 17.58 
18.98+.02 19.25 
20.67 +.01 20.90 




















mole potential barrier, the experimental entropy 
is perfectly reproduced, as shown in Table VIII 
where the experimental entropy increase from 
225.35°K to 298.1°K was calculated from a 
graphical integration of the now available experi- 
mental heat capacity data (including hot wire 
data at the lower temperatures). The correction 
for the hindrance of internal rotation was made 
with Crawford’s treatment.!* It becomes there- 
fore unnecessary to suppose that in crystals of 
propylene end-over-end randomness exists, a 
randomness which is difficult to understand in 
view of very different properties of the CH: and 
CH; terminal groups. Pitzer has_ recently 
criticized our selection of the higher potential 
barrier for propylene because he is unable to ob- 
tain agreement between calculated equilibrium 
constants in the propane, propylene and hydro- 
gen reaction and those measured by Frey and 
Huppke.”’ However he is using an unquestionably 
wrong vibrational assignment for propane (see 
the following section) and if the assignment here 
proposed and tested on the heat capacity data for 
propane is taken, the statistical equilibrium con- 
stant is found to be 0.0014 at 673°K while the 
experimental value is 0.0022. The discrepancy 
between these figures is not as serious as Pitzer 
leads one to believe. Even the determination of 
the ethane, ethylene equilibrium by Frey and 
Huppke gave too high equilibrium constants” 
and the authors themselves point out that less 
satisfactory results were obtained for higher 
hydrocarbons. 


PROPANE (C3Hs) 


The sample of propane was obtained by the 
hydrogenation of 99.9 percent propylene of the 


* Petroleum Symposium, Am. Chem. Soc. Meeting in 
Cincinnati, 1940. 

%F, E. Frey and W. F. Huppke, Ind. Eng. Chem. 25, 
54 (1933). 

* Smith and Vaughan, J. Chem. Phys. 3, 341 (1935). 


same material described above and according to 
the method of Kistiakowsky and co-workers.”5 In 
order to secure as pure a sample of propane as 
possible the first portion of the hydrogenation 
was entirely discarded, a sample was then col- 
lected and tested for olefins. The test showed 
about 0.002 mole percent of unsaturated com- 
pounds present. In obtaining the sample extreme 
care was taken to allow no moisture to contami- 
nate it. The sample of propane was transferred 
to the vacuum system and subjected to repeated 
bulb to bulb distillations. The middle portion of 
the sample only was used. Table IX presents the 
summary of the experimental data on propane. 

Sage, Webster and Lacev“* have measured the 
heat capacity of propane over a similar tempera- 
ture range as given in Table IX. Table X gives a 
comparison of the heat capacities determined by 
Sage, Webster and Lacey and our measurements 
interpolated to their temperatures by graphical 
means and calculated to 1 atmos. pressure of real 
gas. The molecular weight of propane was taken 
as 44.06 in computing Sage, Webster and Lacey’s 
data to our scale. We are at a loss to explain the 
large discrepancy between the two sets of 
measurements which is many times greater than 
our experimental errors including the gas cor- 
rection. 

Kemp and Egan” have measured the third 


TABLE X. Heat capacity of propane. 








Cp (K. AND R.) 


17.41 
18.25 
19.04 
19.81 
20.61 


T°K 


294.31 
311.00 
327.60 
344.30 
361.00 


Cp (S. W. AND L.) 


17.83 
18.29 
18.77 
19.27 
19.78 











% Kistiakowsky, Romeyn, Ruhoff, Smith and Vaughan, 
J. Am. Chem. Soc. 57, 65 (1935); ibid. 57, 876 (1935). 

—_ Webster and Lacey, Ind. Eng. Chem. 29, 1309 
(1937). 

% Kemp and Egan, J. Am. Chem. Soc. 60, 1521 (1938). 
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law entropy of propane obtaining for the ideal 
gas at one atmosphere pressure and 231.04°K the 
value of 60.45+0.10 E. U. Their statistical calcu- 
lations for the entropy showed that the experi- 
mental value for the entropy and the statistical 
calculations could be brought into agreement 
with the assumption of a potential barrier of 3300 
cal./mole restricting rotation of the methyl 
groups; this confirmed the conclusions drawn by 
Pitzer.!* Kistiakowsky, Lacher and Ransom!* 
found from low temperature heat capacity 
measurements that a barrier of 3200 to 3600 
cal./mole gave the best fit for the heat capacity 
data, but that the vibrational assignment used 
by Kemp and Egan and that given by Pitzer 
both gave systematic trends of the barrier with 
temperature, indicating that the vibrational 
assignments were in error. As Table XI shows, 
neither of these vibrational assignments is con- 
sistent with the present heat capacity data. In 
the third and fourth columns of this table have 
been given vibrational heat capacities and those 
due to internal rotation on the basis of Pitzer’s 
assignment, while columns five and stx reproduce 
similar calculations taking the Kemp and Egan 
choice of frequencies. Variable internal rotation 
restricting potentials reaching 10,000 cal./mole 
are necessary to reconcile these data with experi- 
ment,?® as will be seen on inspecting the last 
column in which are given heat capacity contri- 
butions from the restricted rotation of two 
methyl groups with a potential barrier of 3300 
cal./mole. 

Kohlrausch and Koppl®® have studied the 
Raman spectrum and Bartholome*®® has studied 
both the Raman and the infra-red spectrum of 


TABLE XI. Calculated heat capacity of propane. 








Cyib . 
ASSIGN- Cir. cale 
MENT K V =3300 
BY K. CAL./ 
Cir.) ann E. | Cir.| R. | Cir] Move 


3.42] 4.938 . . 3.94 4.04 
kd A ; 3.52 . e . 4.18 4.16 
334.05 . é 3.50 A . x 4.36 4.26 
368.55 s B 3.40 J ° -571| 4.38 4.32 



































** The use of the experimental data of Sage, Webster 
and Lacey would make the situation still less satisfactory. 

* Kohlrausch and Koppl, Zeits. f. physik. Chemie B26, 
209 (1934). 

® Bartholome, Zeits. f. physik. Chemie B23, 152 (1933). 
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TABLE XII. Theoretical and experimental entropy of propane. 








231.04°K 


Strans +rot assuming free rotation (as 


given by Kemp and Egan) 62.80 65.33 
Svib 1.03 2.05 


298.1°K 





63.83 67.38 
60.45 64.62 
Scale —Sexp 3.38 2.76 
Internal barrier from Scajc —Sexp and 


Stotal calc 
Sexp 


Pitzer's tables 3200 cal./mole} 3100 cal./mole 














propane. With these observed frequencies and by 
a comparison with other molecules we have ar- 
rived at the following assignment for the vibra- 
tional frequencies for propane: 373(1), 870(1), 
940(1), 3000(8), 1450(7), 735(1), 1150(1), 1278(1), 
1050(2), 1250(2). Insofar as was possible, only 
the observed frequencies have been used. How- 
ever the last two frequencies have been chosen to 
fit the internal rotational heat capacity curve 
calculated on the basis of a barrier of 3300 cal./ 
mole, using Pitzer’s tables. The results are shown 
in columns seven and eight of Table XI. It is seen 
that the agreement is quite good, but of course 
cannot be regarded as conclusive until the assign- 
ment has been checked by a normal coordinate 
treatment of the molecule. 

The theoretical entropy calculations given by 
Kemp and Egan” for propane were corrected at 
231°K and 298°K using our frequency assign- 
ment. We have also integrated our heat capacity 
curve (using hot wire data at lower tempera- 
tures) in order to get the experimental entropy at 
298.1°K. Table XII gives the summary of these 
calculations. 

In concluding this article we should like to 
point out that surprisingly good agreement be- 
tween experimental heat capacities and statis- 
tical calculations using the harmonic oscillator, 
rigid rotator model are obtained for molecules 
having no internal rotations (cf. data on allene, 
cyclopropane and methyl acetylene). The dis- 
agreements between similar calculations for 
molecules having internal rotation, assuming free 
rotation, and experimental data, are to be taken 
therefore as positive evidence for the existence 
of potential barriers restricting internal rotation. 
An exact calculation of such barriers, however, is 
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impossible unless a complete vibrational analysis 
is available and the Third Law Entropy data are 
preferable for this purpose in case of more com- 
plex molecules since the vibrational contribution 
to the entropies is relatively small. But only by 
combining such entropy data with experimental 
heat capacities can one be assured of a reasonably 
accurate vibrational assignment and hence of a 
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reliable as well as a precise value for the internal 
potential. 

It is a pleasant duty to thank Professor E. B. 
Wilson, Jr., for his interest and many valuable 
contributions to this work. We would further like 
to express our appreciation to Mr. E. E. Roper 
of this laboratory for the use of his equation of 
state data before publication. 
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The present paper presents a continuation of measurements on the gaseous heat capacities by 
the adiabatic expansion method. The apparatus and the’ experimental procedure are exactly 
the same as described in the previous papers denoted herein as Part I and Part II. The com- 
pounds with which the present paper deals are dimethyl ether, ethylene oxide, dimethyl acetv- 
lene, cis-butane-2 and trans-butene-2. For the correction of the experimental data to the ideal 
gas state several procedures had to be used, as discussed in Part II. 


OR dimethyl ether critical data are available* 
and they were put into the Keyes equation 
of state* to obtain the second virial coefficient 
(B) and its derivatives. As pointed out in Part IT, 
the resultant correction of the heat capacity data 
to atmospheric pressure (C,) comes out too low, 
by as much as 10 to 25 percent, and therefore 
we do not list this heat capacity in the following 
tables. In the further correction to the ideal state 
(C,°) the error is fairly well compensated. 

For ethylene oxide only the critical tempera- 
ture’ is known and the following somewhat 
dubious, procedure was resorted to. Maass and 
Boomer® have determined the gas density of 
ethylene oxide at several temperatures and 
pressures. Their data are not of sufficient 
accuracy to evaluate the second virial coefficient 
(B) as a series expansion in temperature but a 
satisfactory average value of (V.—RT/P) could 


1 Part I, G. B. Kistiakowsky and W. W. Rice, J. Chem. 
Phys. 7, 281 (1939). 

? Part II, G. B. Kistiakowsky and W. W. Rice, J. Chem. 
Phys. 8, 610 (1940), preceding article in this issue. 

3 Int. Crit. Tab., Vol. 3. 

4F. G. Keyes, J. Am. Chem. Soc. 60, 1761 (1938). 
asm and Boomer, J. Am. Chem. Soc. 44, 


1726 


be obtained. This, together with the critical 
temperature, was substituted into. the Keyes 
equation? and the correction to the ideal state 
thus calculated. The critical pressure calculated 
on this basis is found to be 49.1 atmos. as a 
weighted average. 

For dimethly acetylene also only the critical 
temperature 7 =489°K is known.* We have 
used the plot of critical temperatures vs. critical 
pressures given in Part II to estimate the latter 
as 31.3 atmos.’ 

Accurate equations of state were available for 
both butene-2 isomers® and were used directly. 
For these latter gases, therefore, the errors due 
to the correction are quite small, but in the other 
cases also, the final data should be accurate to 
better than 0.1 cal./mole degree, except perhaps 
in the neighborhood of the boiling points. 

In the following sections the gases will be 
considered one by one. 


6 Maass and Morehouse, Can. J. Research 11, 640 (1934). 

7Osborne, Garner and Yost (J. Chem. Phys. 8, 131 
(1940)), estimate this by a different procedure and find 
41 atmos. 

8A. B. Lamb and E. E, Roper, J. Phys. Chem. (to 
appear shortly) 1940. 
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DIMETHYL ETHER (C;2H,O) 


The sample of dimethyl ether used in these 
measurements was taken from a cylinder having 
a stated purity of 99.95 percent and obtained 
from the Ohio Chemical Company. It was sub- 
jected to repeated vacuum distillations to elimi- 
nate small traces of air and moisture which 
might be present. 

The following Table I gives a summary of the 
experimental results obtained. As emphasized in 
Part II, none of the completed experiments have 
been disregarded in forming the averages of this 
table because of too large deviations from the 
mean. The same is true of other heat capacity 
tables given in this paper. 

In Table I, 74,°K is the average temperature 
of the gas before and after expansion for each 
series of runs. The column headed (range of AP), 
gives the maximum and minimum values of the 
pressure drop of the gas upon expansion. PrinaiAv. 
is the average barometric pressure for the series 
of runs. The temperature coefficient of resistance 
of the Wollaston wire thermometer is given in 
column five. Column six records the radiation 
correction (see Part I). Column seven shows the 
apparent heat capacity of the gas (C,’) calcu- 
lated as in Part I, and finally the last column 
gives the heat capacity of the gas at zero pressure 
(C,°) which is to be used for comparison with 
the theoretical calculations. 

Crawford and Joyce® have studied the infra- 
red spectrum of dimethyl ether and combining 
this with the available Raman data have given a 
tentative assignment of the vibrational fre- 
quencies: 2900(6), 1466(6), 1180(4), 1122(1), 
940(1), 412(1); all but the lowest frequency 
412 cm— were taken from the infra-red spectrum. 
The numbers in the parenthesis indicate the 
statistical weights attached to each frequency. 
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Two degrees of freedom are allowed for the 
rotation of the methyl groups. ” 
We have calculated the theoretical heat capaci- 
ties on the above assignment, at the temperatures 
of the experimentally determined heat capacities, 
assuming the harmonic oscillator and classical 
rigid rotator as the molecular model. The results 
are shown in Table II. 
In Table II, the third column gives the theo- 
retical heat capacities calculated from Craw- 
ford’s assignment exclusive of the internal rota- 
tion of the methyl groups, (C;j.r.). Column four 
gives the contribution to the heat capacity due 
to the internal rotation of the methyl groups, 
OF Cn een) — (CO veate) ~ Ci-r.). The last column 
gives the contribution to the heat capacity for 
two methyl groups calculated from Pitzer’s 
tables” for a barrier of 2500 cal./mole. The 
moments of inertia used were: J,°=21.666, 
I= 82.474, IP=93.604 and 5.18, the latter being 
for the methyl group. These moments are in 
units of gem? 10” and were calculated on basis 
of the dimensions given by Pauling and Brock- 
way." J,° is the moment of inertia about the axis 
through center of gravity and perpendicular to 
the carbon plane. J,, the moment about axis 
bisecting the C—O—C angle in carbon plane. 
22=moment about axis through center of 
gravity and perpendicular to Y+-Z axis. The cor- 
relation with experimental data is seen to be 
excellent if a potential of 2500 cal./mole is taken 
to restrict the rotation of the methyl groups. 


TABLE II. Internal rotational heat capacity of dimethyl ether. 








Cp® Cir. calc 
T°K (Exp.) (Cp°—Ci.reale Cir. |(V =2500 CAL./MOLE) 


272.20} 14.82 10.79 4.20 
300.76} 15.75 11.59 4.24 
333.25 | 16.81 12.63 4.20 
370.42 | 17.96 13.89 4.12 























TABLE I. Heat capacity of dimethyl ether. (C2H,O). 








RANGE OF Pfnal 


Tay °K AP MM AVERAGE MM 


dR OHM 


i. Rap. Corr. Cp® 
dT DEGREE 


PERCENT C,’ CAL./MOLE DEG. 





37.609—33.882 757.10 
38.385-33.611 759.73 
38.168-32.864 760.80 
37.163-30.563 746.60 


272.20 
300.76 
333.25 
370.42 














14.52+.02 14.82 
15.49+.02 15.75 
16.59 +.02 16.81 
17.77+.02 17.96 


1.1182 0.107 
1.1108 113 
1.1008 141 
1.0871 174 

















*B. L. Crawford, Jr., and L. Joyce, J. Chem. Phys. 7, 
307 (1939). 


10K. S. Pitzer, J. Chem. Phys. 5, 469 (1937). 
1 Pauling and Brockway, J. Am. Chem. Soc. 57, 2684 
(1935). 
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TABLE III. Heat capacity of ethylene oxide (C2H4O). 








P final 
AVERAGE MM 


RANGE OF 
AP mM 


dR Onn 


ox Rap. Corr. Cp 
dT DEGREE 


PERCENT Cy’ CAL./MOLE DEG. 





36.712-30.090 
37.778-31.342 
37.010-33.173 


762.05 
763.70 
764.30 














0.157 
183 
219 


1.1090 
1.0994 
1.0868 


11.57+.01 
12.57+.01 
13.76+.03 


11.80 
12.79 
13.96 

















The only significant deviation of calculated from 
experimental values is at the lowest temperature 
(272°K). Here, however, the correction to the 
ideal state is least certain and the discrepancy 
may be due to this cause. It has been shown” 
that the more complete calculation involving the 
consideration of the coupling of the angular 
momenta of the top and of the whole molecule 
may give results differing appreciably from the 
values given in Pitzer’s tables for the heat 
capacity contribution due to restricted rotation ; 
therefore the value of the restricting potential 
used in Table II should be regarded as pro- 
visional. 


ETHYLENE OXIDE (C2H,O) 


The sample of ethylene oxide used in these 
measurements was a commercial sample ob- 
tained from the Matheson Chemical Company 
and had a stated purity of 99.9 percent. The 
sample was subjected to repeated vacuum dis- 
tillations and pumping to remove traces of air 
and moisture. Only the middle portion of the 
sample was used for actual measurements. Since 
ethylene oxide boils at 10.7°C, the heat capacity 
was determined at only three temperatures. 
Table III gives the experimental results. 

The infra-red and Raman spectra of ethylene 
oxide have been investigated by Linnett!* who 
has given the following tentative vibrational 
assignment: 3000(2), 1494(1), 1122(1), 1270(1), 
811(1), 1453(1), 1163(1), 868(1), 3062(2), 1172(1), 
673(1), 807(1), (1), cm~, with the 811(1), 
1172(1), and a frequency which may be corre- 
lated with the methylene rocking motion being 
uncertain. 

We have calculated the theoretical heat capaci- 
ties for ethylene oxide assuming, as suggested by 
Linnett, a value of 1350 cm~ for the missing 
methylene rocking frequency. The following 
Table IV gives the results. It is seen that the 


2 B. L. Crawford, J. Chem. Phys. 8, 273 (1940). 
13 J. W. Linnett, J. Chem. Phys. 6, 692 (1938). 


agreement between the experimental and the 
theoretical values is not good. The discrepancy 
cannot be attributed to an error in the correction 
of the experimental heat capacities to the ideal 
state, as it is around 0.4 cal./mole degree, 
whereas the total gas correction as shown in 
Table III is only around 0.25 cal./mole degree 
and its uncertainty is at most about 0.1 cal./mole 
degree. 

We have changed therefore the 811 cm"! 
frequency to 1100 cm™ because it is only a 
matter of conjecture as to whether there are 
two frequencies at 807 and 811 cm™ or only one. 
The other two uncertain frequencies namely 1172 
and 1350 cm remain unchanged. Column four 
in Table IV gives the corrected calculations. 
The agreement with the experimental data is 
seen to be quite good but until heat capacity 
data over a larger temperature range are avail- 
able, one cannot be sure that the new assignment 
is substantially correct. 


DIMETHYL ACETYLENE (C4He) 


The dimethyl acetylene used in these measure- 
ments was part of the same material which was: 
prepared and used in the heats of hydrogenation 
by Kistiakowsky and co-workers.!* The material 
was vacuum distilled before introducing it into 
the expansion chamber. Because of the high 
boiling point of dimethyl acetylene, the heat 
capacity was determined at only two tempera- 
tures. The experimental results are given in 
Table V. 


TABLE IV. Statistical heat capacity of ethylene oxide. 








Cp® Cp® 
Cp° LINNETT’S NEw 
(EXP.) ASSIGNMENT ASSIGNMEN1 


11.80 11.80 
12.79 12.78 
13.96 13.94 


— 


4 J. B. Conn, G. B. Kistiakowsky and E. A. Smith, J. 
Am. Chem. Soc. 61, 1868 (1939). 


T °K 


307.18 
337.04 
371.23 





12.16 
13.18 
14.34 


























rae 
OLE DEG. 


80 
1.79 
}.96 








nd the 
epancy 
rection 
e ideal 
degree, 
own in 
degree 
|. /mole 


| em! 
only a 
ere are 
ly one. 
ly 1172 
an four 
lations. 
data is 
apacity 
2 avail- 
znment 


easure- 


ich was’ 


enation 
naterial 
it into 
ie high 
1e heat 
-mpera- 
ven in 








1.80 
2.78 
3.94 


——————— 
———_— 


Smith, J. 








GASEOUS HEAT 





CAPACITIES. 


TABLE V. Heat capacity of dimethyl acetylene (C4Hs). 










Ill 























No. oF | | RANGE OF Pfnal dR _OuM Rap. Corr. Cp® 
RUNS Tw °K AP mM AVERAGE MM dT DEGREE PERCENT Cy’ CAL./MOLE DEG. 
5 336.07 40.796-33.643 771.55 1.0998 0.162 19.41+.03 20.21 
369.46 40.102—34.384 770.90 1.0874 .200 20.76+.04 21.43 











Crawford'® has studied the infra-red and 
Raman spectra of dimethyl acetylene and has 
given an assignment for the vibrational fre- 
quencies of this molecule: 725(1), 1380(1), 
2270(1), 2916(1), 1126(1), 1380(1), 2976(1), 
213(2), 1050(2), 1468(2), 2976(2), 371(2), 1029(2), 
1448(2), 2966(2). 

Crawford and Rice'® using this assignment and 
the heat capacities given in Table V have calcu- 
lated the contribution to the heat capacity due 
to the internal rotation of the methyl groups 
relative to each other and find the restricting 
potential to be at most 500 cal./mole. Table VI 
reproduces these calculations. The experimental 
results for (C;j.2.), the contribution from the 
internal rotation, are in good agreement with 
the value 0.99 cal./mole degree expected for a 
completely free rotator, or a restricting potential 
of zero; 500 cal./mole is undoubtedly the upper 
limit. 

Recently Yost et al.? have published measure- 
ments on the third law entropy of dimethyl 
acetylene, from which they also conclude that 
the restricting potential in this molecule is very 
low, possibly zero. Their conclusion is less de- 
pendent on a correct vibrational assignment than 
the evidence given here and together the two 
studies offer indisputable proof of a practically 
complete absence of hindrance to the rotation of 
the methyl groups in this molecule. It seems to 
us that the theory of restricting potentials 
offered some time ago by Gorin, Walter and 


TABLE VI. Internal rotational heat capacity of 
dimethyl acetylene. 


























Cir. calc 
— Cp° (V =500 
TavK (ExP.) (Cp® —Ci.r eale Cir. CAL./MOLE) 
336.07 | 20.21+.1 19.14 1.07+.1 1.16 
369.46] 21.43+.1 20.47 .96+.1 1.14 























* B. L. Crawford, Jr., J. Chem. Phys. 7, 555 (1939). 
*° B. L. Crawford, Jr. and W. W. Rice, J. Chem. Phys. 
7, 437 (1939). 






Eyring!’ will find it rather difficult to account 
simultaneously for the existence of a rather high 
potential in dimethyl ether and the absence of 
same in dimethyl acetylene. 


CIS-BUTENE-2 (C4Hs) 


The sample of cis-butene-2 was part of the 
original material prepared by Kistiakowsky and 
co-workers for the determination of the heats 
of hydrogenation. The sample was vacuum dis- 
tilled several times before being introduced into 
the expansion vessel but after the experiments 
were carried out it was found that a small trace 
of water vapor remained in the sample. Because 
of this the heat capacity values given in Table 
VII may be 0.1 to 0.2 percent too low. The 
experimental results are shown in Table VII. 


TRANS-BUTENE-2 (C4Hs) 


The sample used was also part of the material 
prepared by Kistiakowsky and co-workers. It 
was subjected to the usual vacuum distillations 
and was definitely freed of water vapor. The 
experimental results are shown in Table VIII. 

No vibrational analysis of the two isomeric 
butenes has been suggested as yet, although 
Gershinowitz and Wilson'® have studied the 
infra-red and Raman spectra of these molecules. 
The important conclusion which can be drawn 
from their work is that the spectra of the two 
molecules have very different appearance and 
the present measurements, which show a differ- 
ence in vibrational plus internal rotational heat 
capacities of as much as 12 percent, bring 
convincing evidence that the fundamental fre- 
quencies of the two isomers are significantly 


1 Gorin, Walter and Eyring, J. Am. Chem. Soc. 61, 
1876 (1939). 

18 Kistiakowsky, Ruhoff, Smith and Vaughan, J. Am. 
Chem. Soc. 57, 876 (1935). 

19 Gershinowitz and Wilson, J. Chem. Phys. 6, 247 
(1938). 


G. B. 


KISTIAKOWSKY AND W. W. 


RICE 


TABLE VII. Heat capacity of cis-butene-2. (C4Hs). 








RANGE OF 


AP MM AVERAGE MM 


Pfinal dR Oum 
dT DEGREE 


Rap. Corr. 
PERCENT Cy’ 


p 
CAL./MOLE DEG. 





38.424-26.831 
38.088-32.972 
38.554-32.275 


759.20 
769.00 
769.40 











1.1114 
1.1009 
1.0868 


19.39 
21.09 
23.01 


0.135 18.47+.05 
.164 20.37 +.03 
157 22.46+.05 




















TABLE VIII. Heat capacity of trans-butene-2. (C4Hs). 








Pfnal 
AVERAGE MM 


RANGE OF 
AP MM 


dR OHM 
dT DEGREE 


Rap. Corr. Cp® 
PERCENT C,’ CAL./MOLE DEG. 
P 





38.472-35.244 
38.350-33.436 
38.322-36.302 


766.50 
767.75 
765.95 














1.1114 
1.1009 14 
1.0866 aa 


0.13 20.04+.03 20.98 
21.96+.03 22.69 


23.97 +.04 24.53 




















different. The neglect to allow for such differ- 
ences in the schematized frequency assignments 
occasionally proposed”® leads evidently to sig- 
nificant errors in the comparative values of the 
calculated thermodynamic functions. 

A different approach to this problem of the 
thermodynamic functions of the more complex 
organic molecules is through a frankly empirical 
set of recursion type formulas derived from ex- 


20K. S. Pitzer, J. Chem. Phys. 5, 473 (1937); M. 
Huggins, ibid. 8, 181 (1940). 


perimental measurements on the gaseous heat 
capacities. The data collected in the first three 
articles of this series show indeed the possibility 
of setting up such equations but we prefer to 
await the accumulation of a little more experi- 
mental material on heavier hydrocarbons, the 
work now in progress in this laboratory, before 
claiming their general validity. 

Professor E. B. Wilson, Jr. has advised us on 
several phases of the present work, for which 
we wish to express to him our appreciation. 
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Statistical Mechanics of Cooperative Phenomena 


Joun G. KirKwoop 
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A direct method of successive approximation to the configurational partition function, based 
on the use of local free energies, is described. The method serves to clarify the assumptions 
underlying the Bethe method of local configurations and the quasi-chemical method of Guggen- 


heim and Fowler. 





HE theoretical treatment of cooperative 
phenomena initiated by Bragg and Wil- 
liams! has received extension by statistical me- 
chanical methods at the hands of several investi- 
gators. Three methods of approximation to the 
configurational partition function have been 
proposed, the method of local configurations by 
Bethe,” the method of moments by the writer,’ 
and the quasi-chemical method by Guggenheim 
and Fowler.* ® The method of local configurations 
and the quasi-chemical method led to closed ex- 
pressions for the partition function, but both are 
based upon approximations, the nature of which is 
not entirely clear. The method of moments leads 
to an exact expansion of the logarithm of the 
partition function in powers of the reciprocal 
temperature, at the sacrifice of closed expressions 
for the thermodynamic functions of the system. 
Bethe and the writer® have recently shown that 
the method of local configurations is valid to the 
third moment of the configurational energy. 
Moreover, Fowler and Guggenheim® have de- 
monstrated the equivalence of the method of local 
configurations and the quasi-chemical method. 
In the present article a method of successive 
approximation to the partition function based 
on the use of local free energies is described. In 
the light of this method, the nature of the as- 
sumptions underlying both the Bethe method of 
local configurations and the quasi-chemical 
method is clarified. In addition, closed expres- 
sions for the configurational entropy and free 
energy of solid solutions are developed. These 
' Bragg and Williams, Proc. Roy. Soc. A145, 699 (1934) ; 
A151, 540 (1935); A152, 231 (1935). 
* Bethe, Proc. Roy. Soc. A150, 552 (1935). 
* Kirkwood, J. Chem. Phys. 6, 70 (1939). 
* Guggenheim, Proc. Roy. Soc. A169, 304 (1938). 


FE and Guggenheim, Proc. Roy. Soc. A174, 189 
40). 


° Bethe and Kirkwood, J. Chem. Phys. 7, 578 (1939). 


expressions were first presented by Bethe and 
the writer® without proof and were later proposed 
by Fowler and Guggenheim! on the basis of the 
quasi-chemical method. 

Although the general features of the method to 
be described are applicable to the calculation of 
the partition functions associated with other 
types of cooperative phenomena, we shall restrict 
ourselves here to the problem of order and dis- 
order in solid solutions. We therefore consider a 
binary solid solution of arbitrary composition 
built up of two interpenetrating lattices A and B 
of N sites each. An arbitrary site of either lattice 
is assumed to have as neighbors z sites of the 
opposite lattice. The lattice sites are supposed to 
be occupied by atoms of two types, denoting by 
N;4, Ni, N24, N2® the numbers of atoms of the 
respective types on the two lattices. We denote 
the corresponding atom fractions by x:“, «1°, 
x24, and x28. To facilitate the calculation of the 
partition function we introduce the long range 
order parameter of Bragg and Williams, 


s=x,4-—x,'; x14=N;4, N=xi+5s 2 
xyB=x,—s/2; xeA=x2.—s/2 (1) 
x22 =x2+5/2, 


where x; and x2 are the atom fractions 
(NiA+N,1®)/2N and (N2A+N28)/2N. of the 
crystal as a whole. If f(s) is the partition function 
of the crystal in a specified state of long range 
order s, the configurational partition function is 


f= ¥ f(s), (2) 


s=—]1 


where the sum extends over all accessible values 
of s between zero and unity. Since f(s) possesses 
a sharp maximum, log f may be set equal to the 
logarithm of the maximum term in the sum (2) 
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with an error of negligible order of magnitude 
0(N-*). In order to calculate the partition func- 
tion f(s) corresponding to a fixed value of s, a 
unique specification of the set of configurations 
corresponding to this value of s is required. Such 
a specification is provided by a set of variables 
£,---&y, m1, ***, nv, each limited to a domain of 
two values, zero and unity. The variable £, states 
the number of atoms of type 1 in site a of lattice 
A and the variable 7, the number of atoms of 
type 1 in site } on lattice B. As a consequence of 
their definition, the configuration variables 
satisfy the following relations, 


N N 
ms &,=N,4, a m= N14. (3) 
a=1 b=1 


If interactions between atoms which are not 
nearest neighbors are neglected, the configura- 
tional energy of the crystal may be expressed in 
the form 


E=Eo—akTo(&, «++ ny), 
Eo= Ny2V4,/2+ No V 22, 2, 
a=((Virt Vo2)/2— Vie l/RT, 
N 
Lal Ea( 1 = no) +(1 = £.)n |, 


a, b=1 


where \q» is unity if sites a and b are nearest 
neighbors and zero otherwise. 


N N 
pas ‘= 2» ab = ( 5 ) 


We may therefore write the configurational par- 
tition function in the following form, 


fsya= Lo LY extrem, (6) 


Er EN m°"7N 
0 0 
Léa = N14 Sm = NF 


From Eq. (6), we obtain by differentiation 
d log f/da=(e)m, 
(o\n=( ¥ oee*)/( X en), 
E1---N f1-++N 


(7) 


where the summation variables in Eq. (7) as in 
Eq. (6) are subject to the restraints, Eq. (3). 
When a is zero, the partition function reduces to 
the total number of configurations consistent 
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with the specified value of s. Using Sterling’s 
formula for factorials, we may write 


log f(s,0)=—N > [xn4 log xn4 +2, log x, ]. 
n=1 (8) 


Integration of Eq. (7) with the use of Eq. (8), 
leads to the following expression for the free 
energy, A(s,a). 


A/NkT=¥ [xp4 log x,4-+x,2 log x," ] 
n=l 


-(1/N) J (o)nda. (9) 
0 


Equation (9) provides an exact expression for the 
configurational free energy, insofar as coupling 
between configuration and lattice vibrations 
can be neglected. Since all pairs of nearest neigh- 
bor sites are equivalent, Eqs. (4) and (7) yield 


(o)w=Nz{l&1—9JwtLl(l—&)nda}, (10) 


where [&(1—7n) ]4 and [(1—£)]w are mean 
values of the specified functions of & and 7 for 
an arbitrary pair of neighbor sites. The mean 
value of any function, y(£, 7), of two of the con- 
figuration variables only may be written as fol- 
lows, 

1 


(v)w= Dd v(&, ne“ VE w/kT, 


§, n=0 


WE, n) =A(é, 9)—A, 


1 1 


> > ero(kt, ++ -mN 


&rEN-1 m"7N —1 
0 0 


eA. n)/kT — 


N-1 N-1 
2 ta=NiA-E © m=NiF 

a=l1 b=l 
where W(é, 7) is the local free energy of a single 
pair of neighbor sites of specified configuration 
[é,], equal to the work of formation of the 
specified local configuration from the average 
crystal. The derivative of W(é, 7) with respect to 
a is easily calculated from Eqs. (7) and (11). 


1 OW(E, n, a) 


ET oo -” Lo(E, n) Iw Lo Jw. 


When a is equal to zero, — W(é, n, a) /RT reduces 
to the logarithm of the ratio of the number ol 
configurations, consistent with the specified 
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values of ~ and 7, to the total number of con- 
figurations of the crystal, 


— WE, n, 0)/kT = log x14+7n log x1? 
+(1—£) log x24+(1—7) log x28. (13) 
Integration of Eq. (12) with the use of (13) gives 
EWE DIRT = (4, A)E(grg4) IF (— 8) 


X (xB) new, DIET 
(14) 


—w(é, 9) /kT= J {Lo(E, 0) w—[o Iw ida. 


Eqs. (10), (11), and (12) lead to the following 
expression for [o ]a 
[ow = Nef xpAxgBe-e il. O/kT 

+xo4x Fev Det} (15) 
Since &, 7 each have a domain of two values only, 


zero and unity, we may express w(é, 7) as a func- 
tion of these variables in the following form 


w(t, n) = €00— €10€ — Eon + e1é7. (16) 


The four values of the function w(é, 7) are not 
independent, but satisfy three equations which 
may be obtained by substitution from Eq. (14) 
in the following generally valid relations 


1 
dy e VE mikT = 1, 


&, 7=0 
1 
Celw= YE eV EMT I yA, (17) 
§&, n=0 
1 
Unda = Zz ne WE /kT = xB, 
§, n=0 . 


The first of Eqs. (13) follows directly from the 
definition of W(é, 7), Eq. (11). The second and 
third of Eqs. (13) follow from Eqs. (3) and the 
fact that all sites of lattice A or lattice B are 
equivalent in the average crystal. Thus [£a Jw 
and [mo] are independent of the site indices a 
and 6 and equal to x;4 and x, respectively, for 
any value of a. By means of Eqs. (14) and (17), 
the coefficients €10, €01, and €o9 in the expansion 
of w(, n), Eq. (16), may be expressed in terms of 
11 the coefficient of &. By elimination of €,0, 
€o1, and €9 in this manner and substitution of the 
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resulting expressions for w(é, 7) in Eq. (15), we 
obtain after some tedious algebraic transforma- 
tions, 


[o lav 2s?+-(1—s?—q*)et/*7+-(1+5?—g?) 9 
Nz (14+s?—g?)+(1—s*+q2e#?42¢ 
g= {s?+(1—s?@—g2)e!*T + g2e2/ kT}, (18) 





g=X1—X2, s=x,4—x,5, 
e=w(1, 1)+w/(0, 0) —w(1, 0) —w(0, 1), 


where, for simplicity in notation, the symbol 1; 
has been replaced by e. Equation (1) has been 
employed to express x;4, x;", x24, x2? in terms 
of s and g, the latter denoting the difference in 
atom fraction of the two components in the 
crystal as a whole. By Eqs. (14) and (15), we may 
write 
€ ? d*Lo(E, n) Iw 
-—= — ——da. (19) 
kT 0 0£0n 





In order to calculate ¢ as a function of a, it is 
necessary to determine [o(£, 7) Jw. This could be 
done by a continuation of the method employed 
in the calculation of [¢],, with local free ener- 
gies of sets of four sites in specified configurations. 
An alternative procedure consists in the expan- 
sion of ¢ in powers of a by the method of mo- 
ments. Such an expansion leads to the following 
expression for w(£, 7). 


w(é, 7) 
kT 





=| 2e1— (144-98 (1+045)n 
(1+g)?—s? 
$F at 010". (20) 

The corresponding expansion of ¢ is evidently 
e/kT =2a+O(a’). (21) 


When the expansion (21) is substituted in Eq. 
(18), the following approximation to [o]« is 
obtained, 


2s?+ (1—s?—q?)e2#+(1+5?—-@2)¢ 
(1-+s?—q?)+(1—s?+q*)e**+2¢ - 
( 
g(a) = {s?+(1—s?—g?)e2*#+q7e4*} 3, 





[o |a/Ns= 


As we shall presently show, the approximation 
(22) leads directly to the Bethe theory of super- 
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lattices. The Bragg-Williams theory is obtained 
when the total free energy of the crystal is ex- 
panded in powers of a by the method of moments 
with retention of the first moment only. The 
Bethe theory, on the other hand, is obtained by 
introducing the same type of approximation for 
the local free energy «. Thus we may say that the 
Bethe theory is based upon the application of 
the Bragg and Williams approximation to the 
local free energy ¢, rather than to the total free 
energy of the crystal. One may ask why the first 
moment approximation to w(é,7), Eq. (20), 
rather than the corresponding approximation to 
e, Eq. (21), should not be used for the calculation 
of [o Jw in Eq. (15). If this is done, an approxima- 
tion agreeing with the exact expansion of the total 
free energy of the crystal to the second moment in 
energy is obtained. It may be verified that the 
Bethe approximation is valid to the third moment 
in energy, and therefore somewhat better. The 
reason for this is not difficult to understand. 
The w(é,7) of Eq. (20) do not yield W(é, 7) 
satisfying Eqs. (17), whereas these relations are 
satisfied automatically if Eqs. (18) and the first 
moment approximation to ¢e is employed in the 
calculation of [o ]w. 

Although it is possible to express the integral, 
JSv'[o nda, in terms of elementary functions with 
the integrand Eq. (22), the result is very cumber- 
some. We shall therefore give expressions for the 
free energy only for the special cases g=0 and 
s=0. When g=0 the atom fraction of the two 
components are equal. Equations (9) and (22) 
lead to the result; 


A i+s 1+s 
2NkRT 2 2 


1-—s i-s az 
+(+) log (—) +——2G(a, s*), 
2 2 
i+s 2(s+¢) 
G(a, s?) = —{ —— } log 
4 (1+ ¢)(1+s) 


1-s 2(¢—s) 
<i igest 
4 (g+1)(1—s) 
o(a, s) =[s?-+ (1 —s?)e?*]}, 


The equilibrium value of s is to be determined by 
the condition dA/ds=0. A discussion of the 
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order-disorder transformation and the accom- 
panying heat capacity anomaly, based upon Eq. 
(23), has been given by Bethe and the writer. 
The results are identical with those of the in- 
direct method originally formulated by Bethe. 
When G(a, s*) is expanded in powers of a and 
the result is compared with the exact expansion 
of A/2NkT by the method of moments, it is 
found that Eq. (23) is valid to terms in a’, that 
is to the third moment in the energy. 

When a is negative, the long range order s is 
constrained to the value zero. In this case, Eqs. 
(9) and (22) yield 


A (= (= 
= ) log —) 
2NkT 2 2 
1 


i—¢g 








) oe | 
(g+1)(1+q) 


2(¢—q) 





(g+1)(1—q) 
o(—a, g*)=[¢?+(1—g*)e*]}. 


When a is negative, there is instead of a phase 
change of the second kind, one of the first kind 
with a critical solution temperature below which 
the solid solution separates into two phases. The 
critical solution phenomenon has been investi- 
gated by the writer’ using the method of mo- 
ments, and by Rushbrooke,* using the method of 
local configurations. We shall therefore not dis- 
cuss in detail application of the closed expression, 
Eq. (24), to the problem of critical solution phe- 
nomena, but remark that, as in the superlattice 
case, the results agree with the method of mo- 
ments to the third moment in the energy, that is 
to terms in a’ in the free energy. 

The foregoing analysis serves to illuminate the 
nature of the approximations underlying the 
method of local configurations and the quasi- 
chemical method and demonstrates that they are 
both valid to the third moment of the configura- 
tional energy. All three methods, the method of 
moments, the method of local configurations, and 


7 Kirkwood, J. Phys. Chem. 43, 97 (1939). 
8 Rushbrooke, Proc. Roy. Soc. A166, 296 (1938). 
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JOULE-THOMSON EFFECT 


the quasi-chemical method reduce in first ap- 
proximation to the Bragg-Williams form and the 
higher approximations do not give rise to very 
large deviations from the Bragg-Williams theory. 
It is the writer’s opinion that improved agree- 
ment with experiment is not to be sought in re- 
finements in the statistical calculation of the 
partition function, but rather in improvements of 
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the physical model upon which this calculation is 
based. For example, the dependence of the inter- 
action energies Vi, V22 and Vi2 on the lattice 
parameter of the crystal is of considerable im- 
portance. ® 

In conclusion, the writer wishes to express his 
thanks to Professor Bethe for his comments on 
the method proposed in the present article. 
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The Joule-Thomson Effect in Mixtures of Helium and Argon* 


J. R. Roesuck AND H. OsTERBERG 
Department of Physics, University of Wisconsin, Madison, Wisconsin 
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The Joule-Thomson effect has been measured in four mixtures of helium and argon. The 
apparatus and methods are those used repeatedly in previous work. The measurements of the 
isenthalpic curves cover the range — 100° to 250°C and up to 200 atmos. From these measure- 
ments uw has been computed and plotted, and the trends with temperature, pressure, and com- 
position are shown. u is read from these curves and tabulated. In part of the region 1 =0 and 
these inversion curves are plotted as functions of pressure, temperature, and composition. 
Mixtures of helium and nitrogen have been similarly measured. 


HE measurements on the Joule-Thomson 
effect in air,’? in helium,** in argon® in 
nitrogen,® and in four mixtures of helium and 
nitrogen’ are here extended to four mixtures of 
helium and argon. Reference is made to the above 
articles for details of apparatus and methods. 
Our usual pressure range of 1 to 200 atmos. is 
maintained here, but our temperature is reduced 
* The senior author’s attention has been called lately to 
the reprinting by E. Justi in “Spezifische Warme, Enthal- 
pie, Entropie, etc.” of Table XIV from my second article, 
reference 2, on the Joule-Thomson effect in air. As indi- 
cated in that article the source of this data was a large 
chart from H. Hausen. Unfortunately in reading off the 
values from the chart serious errors were made in the 


100-atmos. group. The data for » at 100 atmos. should have 
been as follows: 








—150°C| —140] —120] —100} —75 | —50} —25) 0 | +25°C 


Roebuck | +0.019 | 0.044] 0.165] 0.285] 0.287|0.247/0.215/0.178| 0.152 
Hausen | —0.006 | 0.016] 0.164/ 0.283] 0.308|0.268/0.227|0.186) 0.158 






































The agreement is markedly improved. 

! Roebuck, Proc. Am. Acad. 60, 537 (1925). 

? Roebuck, Proc. Am. Acad. 64, 287 (1930). 

® Roebuck and Osterberg, Phys. Rev. 43, 60 (1933). 

‘ Roebuck and Osterberg, Phys. Rev. 45, 332 (1934). 

’ Roebuck and Osterberg, Phys. Rev. 46, 785 (1934). 

® Roebuck and Osterberg, Phys. Rev. 48, 450 (1935). 

* Roebuckand Osterberg, J. Am. Chem. Soc. 60, 341 (1938). 


to —100° to 250°C. The experimental work was 
done between May and November of 1935. It has 
seemed advisable to proceed with the collection 
of other data even though this led to consider- 
able delay in publication. Complete data on 
carbon dioxide has now been collected. 

The work on the Joule-Thomson effect in 
carbon dioxide has yielded as well data on the 
vapor pressure of the liquid. This data was 
just enough different from that in the literature 
to arouse suspicion of trouble. Hence the vapor 
pressure was measured statically using the piston 
gauge and the resulting data compared with the 
readings made on the open tube mercury 
manometer® at the ice point, and with the 
measurements of Meyer and Van Dusen.® The 
observed disagreement led to a careful search of 
the records on the piston gauge. In 1912 the 
piston gauge was checked against an open-tube 
mercury manometer very successfully, but in the 
computation of a factor to translate gauge-load 


8 Roebuck and Cram, Rev. Sci. Inst. 8, 215 (1937). 
® Meyer and Van Dusen, Bur. Stand. J. Research, 10, 
381 (1933). 
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TABLE I. Helium and argon mixtures. 








HELIUM % ARGON % 





100 0 
75.8 24.2 
50.6 49.4 
33.5 66.5 
16.6 83.4 

0 100. 








readings into atmospheres a numerical error was 
made. The factor 1.457 should have been 1.404. 
The factor computed from the work on the vapor 
pressure of carbon dioxide at the ice point by the 
open tube mercury manometer and by the piston 
gauge, is 1.410. This value (1.410) was used to 
calculate the vapor pressure in atmospheres read 
at a set of temperatures between 0°C and 30°C. 
All these values checked to less than 0.1 percent 
with readings of Meyer and Van Dusen. The 
factor has thus shifted from 1.404 to 1.410 or 
0.4 percent in the 27-year interval. 

In this interval a series'~’ of papers have been 
published on work employing this piston gauge 
in which the factor 1.457 was always used. The 
numbers representing the pressure in all these 
papers have to be multiplied by 0.9677 to correct 
this error. Suitable corrections must also be 
made in the derived quantities like u. As oppor- 
tunity offers this data will be reworked to make 
these corrections. They have been applied to the 
data of this article. 


LIQUEFACTION 


The presence of liquefaction anywhere in the 
system led to the same difficulties with erratic 
data, as described in the preceding article on 
He-Ne: mixtures.’ For similar reasons no attempt 
was made to take data with liquefying conditions 
present anywhere in the system. The bath 
temperature was in consequence never set below 
about —100°C and of the curves measured at 
this temperature only the helium-rich mixture 
(No. 1) curve was carried to the lowest pressures. 
Until data has been gathered in the simpler 
parts of this field, the more difficult may well 
be avoided. 


Gas SUPPLIES 


The helium was part of that supplied by the 
Bureau of Mines for the previous work on pure 


helium.* It was purified before use by passage 
over activated cocoanut charcoal cooled by 
liquid air as in the earlier work.’ 

The argon was a gift from the Air Reduction 
Sales Company through Mr. F. P. Gross, Jr. 
It had been prepared for use in incandescent 
lamps and was stated to contain 99.5 percent 
argon, 0.5 percent nitrogen, with negligible 
quantities of oxygen and water vapor. 

The mixtures were made up by admitting the 
desired measured quantities of the constituent 
gases to the flow system and mixing thoroughly. 
The mixtures were carefully analyzed by ab- 
sorbing the argon in liquid-air-cooled charcoal.’ 
The two methods of knowing the composition 
agree. The analytical data are given in Table I. 

Since both the He+air and the He+N, 
mixture work had shown that the Joule-Thomson 
effect varied less with composition for the helium- 
rich mixtures, these He+A mixtures were 
crowded somewhat toward the pure argon side. 
Sufficient of each mixture for the projected work 
was made up at once so that the compositions 
should be very constant. To avoid excess loss 
of argon the lubricating water was discharged 
only near atmospheric pressure where the dis- 
solved helium and argon are both small in 
quantity. 


EXPERIMENTAL RESULTS 


The data for mixtures 1-4 are given in Tables 
II-V. The individual runs are referred to by 
their approximate bath temperatures. The last 
temperature and pressure readings in a run are, 
respectively, the bath temperature and the inlet 
high pressure. The temperatures are in degrees 
centigrade and the pressures in atmospheres 
absolute. 

For reasons similar to those given in the 
He+N, mixture article,’ it was not considered 
necessary to duplicate the individual runs. The 
high precision with which the points fall on a 
curve and the consistency of the family of curves 
and of the family of their derivatives are suff- 
cient assurance of the reliability of the data. 


ISENTHALPIC CURVES 


The data of Tables II-V are plotted, respec- 
tively, in Figs. 1-4, with temperature as ordinate 
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and pressure as abscissa. The coincidence of the 
points with the curves in the accurate plots from 
which these figures were traced, leaves very little 
to be desired. 

The corresponding curves for pure helium* 
and for pure argon’ should be considered to- 


TABLE II. Mixture No. 1. 








Z 
9 
Z 
° 


pe —le 


—50° Curve 


1.9 48.38 
19.6 48.43 
42.4 48.64 
64.0 48.88 
84.0 49.19 

107.8 49.60 
130.8 50.06 
153.5 50.55 
194.8 51.61 


—100° Curve 


1.9 101.00 
100.39 

99.79 

99.34 

99.06 

98.86 

98.78 

98.84 

99.00 


pe lo 


250° Curve 
2.5 262.78 
19.6 261.81 
260.46 


pe le 


100° Curve 


3.8 109.61 
19.4 108.91 
41.2 107.82 
63.4 106.71 
85.3 105.56 

107.1 104.42 
130.3 103.18 
153.4 101.89 
194.8 99.73 





CSrmnrIaAUswn— 
Caen auvehlwne 
CeBNAUELWN 


50° Curve 
57.81 
57.33 
56.44 
55.55 
54.63 
53.62 
52.63 
51.60 


COM IAUESWNH— 
CeRNIANEWN— 
CeOIAWEWN— 


1 
2 
3 
4 
5 
6 
7 
8 
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SSRDNIAUEWNS 


— 








TABLE III. Mixture No. 2. 








Z 
° 
Z 
9 


pe le 


250° Curve 
3.8 258.39 
19.5 258.11 
43.3 257.28 
66.9 256.40 
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CONAN EWNe 
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46.22 
47.09 
47.84 
48.46 
48.89 
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gether with these four groups of curves for the 
four mixtures. A comparison of the curves of 
these six figures with the corresponding six 
figures for the He+N:2 mixtures’ shows them to 
be remarkably similar. In the six figures, a 


TABLE IV. Mixture No. 3. 








Z 
9 
Z 
° 


p2 lo 


250° Curve 


4.4 252.62 
19.5 252.69 
42.4 252.87 
64.4 252.73 
87.4 252.57 

108.5 252.35 
132.6 252.03 
155.6 251.66 
194.8 251.01 





—50° Curve 


2.5 92.98 
86.35 
78.57 
72.21 
66.69 
62.07 
58.23 
55.18 
51.25 


CBNAWVSWN— 
SOPMIAUNAEWNHE 
CeONAUEWN— 


— 


200° Curve 


199.36 
199.83 
200.31 
200.54 
200.75 
200.86 
200.85 
200.82 
200.66 


150° Curve 
4.5 145.74 
19.6 146.61 
42.4 147.54 
65.3 148.36 
87.1 148.92 
108.5 149.43 
132.6 149.90 
154.1 150.22 
170.5 150.35 
194.8 150.67 


OCeAULONS 
CSeONAUEWNH— 


SOeMNIAUNE WH 
SOMAIAUNEWNS 


= 
— 








TABLE V. Mixture No. 4. 








Z 
° 
Z 
° 


2 le pe le pe —te 





100° Curve —650° Curve 
75.41 
78.01 
81.51 
84.72 
87.71 
90.26 
92.63 
94.61 
96.50 
98.32 


250° Curve 


4.6 245.57 
19.4 246.44 
42.4 247.55 
65.2 248.42 
86.8 249.13 

109.0 249.74 
131.3 250.30 
154.9 250.76 
194.8 250.77 


CHeIAUSWNH— 
SORNAUEP WN 
SSOmMNAUEWN 


_ 
—_ 


50° Curve 


4.1 16.51 
19.8 20.74 
42.3 25.96 

30.61 
34.78 
38.31 
41.90 
44.93 
49.64 


200° Curve 


4.8 191.11 
22.1 192.33 
42.0 193.94 
66.2 195.47 
87.4 196.67 

109.5 197.80 
131.5 198.76 
154.2 199.59 
177.5 200.41 
194.8 200.62 


SOMA S wre 
CHOUIANUAWN— 


150° Curve 0° Curve 
5.4 134.58 2.5 —49.26 
19.9 136.34 —42.09 
42.0 138.86 —33.88 
66.2 141.22 — 26.56 
87.4 143.07 —20.49 
109.7 144.85 —15.45 
131.3 146.29 —10.38 
154.4 147.74 —6.37 
177.5 148.99 —2.85 
194.8 150.04 —0.18 
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100 ATM. 150 200 


Fic. 1. Mix No. 1. Temperature versus pressure at 
constant enthalpy. 


comparison of the isenthalpic curves for any one 
bath temperature shows how the straight lines of 
negative slope for pure helium shift regularly to 
the concave-downward curves of large positive 
slope for argon. The situation is the same as 
first observed in our mixture work with air in 
helium.’ 

One observes again the marked similarity 
between the family of isenthalpic curves for one 
pure gas (e.g. for argon) at a series of bath 
temperatures, and the family of curves corre- 
sponding to a series of compositions all at the 
same temperature. 

To observe the effect of the composition and 
temperature on the total drop, the At for a 
pressure drop from 200 to 15 atmos. was read 
off from these isenthalpic curves of Figs. 1-4. 
In Fig. 5 these values of Aé are plotted against 
composition giving a set of curves, one for each 
of four bath temperatures. The four curves have 
the same general character as Fig. 5 in the 
helium’ article, and Fig. 5 in He+Ne mixture’ 
article. Here, in Fig. 5, the rise at the left-hand 
end of the curves is not present, the 250° curve 
starting practically horizontally. This is probably 
due to the higher critical temperature of the 
argon over both air and nitrogen. It is probable 
that at bath temperatures higher than 250°C the 
rise would appear in this case also. 
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100 ATM. 150 


Fic. 2. Mix No. 2. Temperature versus pressure at 
constant enthalpy. 














































































































































































































50 100 ATM. 150 200 


Fic. 3. Mix No. 3. Temperature versus pressure at 
constant enthalpy. 


The slope of all these curves increases toward 
the pure argon side of the figure, and also in- 
creases with falling temperature. The data was 
not carried to low enough temperatures to enter 
the region where with pure argon this increase 
reverses. 


JouLE-THOMSON COEFFICIENT, = (dt/dp); 


As in previous work, the numerical values of 
are obtained by taking the ratio of the successive 
and corresponding differences of temperature and 
of pressure for each experimental run. These 
values of u were plotted against the correspond- 
ing average pressure to give a set of isenthalpic 
curves for each mixture. In the large scale plot 
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(25X40 cm) the experimental points drop on 
smooth curves at least as well as with previous 
data, leaving only small irregularities to be 
smoothed out. It has not seemed necessary to 
reproduce these curves. 

From these unpublished curves values of u for 
a series of selected pressures were read off and 
the corresponding temperatures read off from the 
isenthalpic curves of Figs. 1-4. These values of yu 
and ¢ were then plotted as isopiestics. The points 
obtained from successive differences show neces- 
sarily considerably more variation than the 
original isenthalpic curves. This cross plot was 
used to help adjust the yp vs. p curves through the 
experimental points to obtain as consistent 
families as possible, without in any case leaving 
the experimental points farther than the spread 
of these points justifies. 

The four families of curves, one family for each 
mixture, are reproduced in Figs. 6-9. 

The abscissa scale is the same in all four 
figures, but to offset the increase in the average 

















































































































































































































Ee 





50 100 ATM. 150 200 


Fic. 4. Mix No. 4. Temperature versus pressure at 
constant enthalpy. 


Fic. 5. Temperature change for a pressure drop of 200 
to 15 atmos. as a function of composition at several initial 
temperatures. 


value of yw, the ordinate scale is shifted by a 
factor of 2 between each Fig. 6, 7, 8, and 9. 
These families resemble each other less than the 
four similar families of the He+N. mixtures, 
due here to a closer approach to the critical 
temperature of the higher critical-temperature 
component. Some of the curves here enter the 
region where uw decreases with decreasing tem- 
perature. Again, however, each family looks like 
an appropriate section from the family for pure 
air, argon, or nitrogen. The theorem of corre- 
sponding states is again suggested. 

The previous observation for the isenthalpic 
curves of Figs. 1-4, that the effect of a change in 
composition resembles greatly the effect of a 
change in temperature, holds again in these 
families, Figs. 6-9. 

The numerical data for u have been picked off 
the large scale plots from which Figs. 6-9 were 
traced and are collected in Tables VI-IX. With 
mixture No. 1 the negative values of yw at the 
highest temperature and lowest pressure are the 
same as for pure helium; but these values fall 
numerically with rising pressure, while with pure 
helium they remain constant except at our 
lowest attainable temperature. With all these 
mixtures yw decreases with rising temperature 
and falling pressure except for 3 values in Table 
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Fic. 9. Mix No. 4. « versus temperature at 


Fic. 6. Mix No. 1. u versus temperature at constant pressure. 
constant pressure. 





TABLE VI. »X10?; Mix No. 1; He 75.8 percent, 
A 24.2 percent. 





















































\? 
tN. 1 20 60 100 | 140 | 180 | 200 
250 —5.83 | —5.95| —6.15| —6.37| —6.56| —6.77 | —6.85 
200 5.55] 5.66] 5.90] 6.13] 6.34] 655] 6.63 
150 5.11] 5.24| 5.52] 5.77] 5.99] 6.21| 6.34 
100 4.47| 4.61| 4.91] 5.18] 5.45] 5.72] 5.88 
50 3.61| 3.76] 4.08] 4.40| 4.68] 5.01] 5.19 
0 2.40] 2.57] 2.92] 3.30] 3.65| 4.03] 4.22 
—50 —0.69 | —0.92 | —1.32| -1.75] 2.21| 2.66] 2.82 
—100 +3.37 | +2.82 | +1.87| +0.79| -0.14| —0.65 | —0.78 
TABLE VII. »X 10°; Mix No. 2; He 50.6 percent, 
A 49.4 percent. 
\? 
t°c\, 1 20 60 100 | 140 | 180 | 200 
G. i. i . 2. uw versus tem r 250 —2.84| —3.19| —3.65| —4.04) —4.21| —4.33 | —4.34 
Fea. 7. Mix No ae aati 200 1.67} 2.07| 2.71] 3.15| 3.40| 3.55| 3.37 
f constant pressure. 150 —0.13 | —0.67| —1.50} 2.01] 2.32] 2.56] 2.62 
100 +1.84] 41.15] +011] —0.59} —1.01] -1.32] 1.48 
50 4.50} 3.66] 2.37) +1.39] +0.70] +0.14| —0.07 
0 8.19] 7.20] 5.51] 4.12) 2.96] 1.99] +1.57 
—50 13.84] 12.61] 10.27} 8.14] 6.28} 4.53| 3.63 
—100 — — | +17.79| +14.17| +10.36| +6.90 | +5.40 





























TABLE VIII. «10°; Mix No. 3; He 33.5 percent 
A 66.5 percent. 
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Fic. 8. Mix No. 3. u versus temperature at 
constant pressure. 








. TABLE IX. 4X 10?; Mix No. 4; He 16.6 percent, 
IX (Mix. 4) where the influence of the critical A 83.4 percent. 
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77 | —6.85 
55 6.63 
21 6.34 
72 5.88 
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03 4.22 
66 2.82 
65 | —0.78 
cent, 
) 200 
33 | —4.34 
55 3.57 
56 2.62 
32 1.48 
14] —0.07 
99 | +1.57 
53 3.63 
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B | —1.68 
B | —0.38 
2 | +0.83 
" 2.54 
| 5.01 
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Fic. 10. y« versus composition at a series of constant 
temperatures but all at one atmos. pressure. 


pressures of 1, 100, and 200 atmospheres in the 
three figures, respectively. These include also the 
necessary data for pure helium and pure argon. 
The points fall systematically on smooth curves 
as in the corresponding situation for the He+Ne 
mixtures. With the He+A mixtures however the 
curves exhibit less pronounced minima, but from 
the trend of the curves in the three figures it is 
apparent that they will become more pronounced 
at higher temperatures. The critical temperature 
of argon being higher than that of nitrogen is 
apparently a reason for this higher temperature. 
Several of the low temperature curves enter the 
temperature field in which yu for argon shows the 
rapid decrease characteristic of the approach to 
the critical temperature from above. Here evi- 
dently the properties of the mixture are largely 
controlled by those of the argon constituent. 


INVERSION CURVES, n=0 


In three of the four families of curves Figs. 1-4, 
the isenthalpic curves show maxima (u=0). In 
the mixture richest in helium, No. 1, the maxima 
show in the two lowest temperature curves; in 
No. 2 they lie about the middle of the group; 
in No. 3 they lie near the upper part of the tem- 
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Fic. 11. mw versus composition at a series of constant 
temperatures but all at 100 atmos. pressure. 





Fic. 12. « versus composition at a series of constant 
temperatures but all at 200 atmos. pressure. 
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Fic. 13, Pressure as a function of temperature for «=0 
for three mixtures. 


perature range; and in No. 4 the highest tem- 
perature (250°C) is not high enough for them 
to appear. 

The pressures and temperatures at which ».=0 
may be taken off the isenthalpic curves of Figs. 
1, 2and 3. They may also be read from the curves 
of Figs. 6, 7 and 8. These data are plotted in 
Fig. 13. Only the upper branch of each curve 
appears. The curves are much less regular than 
the corresponding group’ for He+Na, particu- 
larly the curve for mixture No. 3. Within this 
uncertainty they confirm the earlier conclusion 
that the curves approach zero pressure without 
any material change of their small curvature. 
The curves of Fig. 14 exhibit the effect of 
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Fic. 14. Temperature as a function of composition for 
u=0 and three constant pressures. 


composition on the temperature at which »=0 
for three pressures, 1, 100, and 200 atmos. The 
curves are not quite so regular as in the He+Ne 
case, but are very similar to them. The situation 
where »=0, shifts to lower temperature with 
increase in pressure or in He proportion. This 
applies of course only to these upper branches of 
the inversion curves. The lower branches fall in 
the general region where condensation occurs and 
were not measured. 

As in the case’ of He+Nez a mixture of He 
and A may be made, the reciprocal of whose 
average coefficient of expansion between 0°C 
and 100°C (ag!) is itself the temperature of the 
ice point on the Kelvin scale. It contains about 
35 percent by volume of argon, whereas it 
required about 46 percent Ne. This smaller pro- 
portion can be attributed to the effect of the 
higher critical temperature of the argon. 

We have estimated‘ the temperature as 23.6°K 
at which pure helium has »=0 at 1 atmos. For 
pure argon® this temperature is above the tem- 
perature range of our experiments. It can be 
estimated however from the nitrogen® data using 
the theorem of corresponding states. In the 
nitrogen® article this temperature for N» is given 
as 334°C (607°K). This gives 373°C (645°K) for 
argon. This data (# against composition for 
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Fic. 15. Temperature as a function of composition for 
u#=0 and p=1 atmos. for several mixtures. 


p=1, u=0) for the range of composition of one 
pure material to the other, and for the mixtures 
He+Air, He+N2, and He+A is plotted in 
Fig. 15. 

The single point for the air mixture coming 
from the preliminary work on He’, falls almost 
exactly on a straight line between the points for 
pure air and pure He. The He+N: points fall 
well on a smooth curve which is straight for the 
first half from pure He, but is concave downward 
for the rest. The temperature is not far from a 
linear function of the composition. 

The He+A points do not fall well on a smooth 
curve and the whole curve is concave downward. 
Like the He+N2 curve the greatest curvature is 
at the higher temperature. While the experi- 
mental data for individual runs is fully as 
consistent as in any of the previous work, there 
appear to be somewhat greater inconsistencies in 
the comparisons of data from the various parts 
of the field than is to be observed in the He+N2 
case. 

In the work with helium? tests for air impurity 
and repurification were relatively simple matters. 
With nitrogen the addition of air disturbs the 
purity more slowly, air being # nitrogen, and 
also the nitrogen being cheap could be discarded. 
But with the argon, air impurity has a relatively 
large effect on its properties, we had no adequate 
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TABLE X. Temperature of inversion for p =0 in reduced units. 











MIx No. 4PC Ti&K T&K Ti/Te 
He+A 1 — 54.1 209 404 5.17 
He+A 2 135 420 77.2 | 5.44 5.44 
He+A 3 290 583 102 3.71 














tests for the presence of air, repurification was 
beyond our available experimental means, and 
the argon was too expensive to discard. All of 
the data in this article were taken on a single 
cylinder of argon by mixing it with helium, 
measuring the argon-rich mixtures first. But also 
the measurements were carried out rapidly and 
with a minimum of difficulty with an apparatus 
in which leaks were small. It is, nevertheless, to 
be expected that the data contain a larger 
element of uncertainty on account of air im- 
purity, than in any of the previous work. 

The hypothetical critical temperature of these 
mixtures was calculated as described in the 
He+Ne article.’ Using these 7. values the 
observed inversion temperatures of the mixture 
were shifted to a reduced temperature scale. The 
data are given in Table X. Table X of the 
article’ on He+Ne mixtures gives the corre- 
sponding data for other mixtures and some pure 
gases. The average value of 5.44 is larger than in 
this previous work where, for example, He+Ne 
mixtures’ give 5.18. 


By chance, reference to the above use of a 
linear variation of 7. with composition, was 
found in a copy of Nernst’s Theoretical Chemistry 
where he cites two references.'® "! 

Our thanks are due to the Wisconsin Alumni 
Research Foundation whose generous assistance 
enabled full time to be devoted to the experi- 
mental part of this work. 

The helium used here and in other work was 
supplied by the Bureau of Mines of the United 
States Government. We thank them for their 
help which made this extensive work with helium 
possible. 

We also wish to thank the Air Reduction Sales 
Company and Mr. F. P. Gross, Jr., for their 
gift of the argon which made this work possible. 

In this as in some other problems, grants from 
the Rumford Fund of the American Academy of 
Arts and Sciences have been very helpful. 

Much of the reduction of the experimental 
data and the preparation of the curves were 
carried out by Mr. H. R. Kelly, whose assistance 
is gratefully acknowledged. Our thanks are also 
due to Mr. T. A. Murrell who carried out the 
major part of the recalculation and redrafting 
necessary for the correction of the pressure 
scale. 


10 Pawlewski, Berichte (Berlin) 16, 33 (1883). 
1G, C. Schmidt, Lieb. Ann. 266, 266 (1891). 
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HE structures of methyl cyanide (CH;CN) 
and methyl isocyanide (CH3NC) have been 
determined by Brockway! by means of electron 
diffraction, and their vibration frequencies have 
been determined by Badger and Bauer? from a 
study of the infra-red and Raman spectra. These 
studies provide all the necessary data for com- 
puting the thermodynamic properties of these 
substances in the ideal gaseous state on the basis 
1L. O. Brockway, J. Am. Chem. Soc. 58, 2516 (1936). 


2 R. M. Badger and S. H. Bauer, J. Am. Chem. Soc. 59, 
303 (1937). 


TABLE I. Molecular constants of methyl cyanide and 
methyl tsocyanide. 














METHYL METHYL 
CYANIDE ISOCYANIDE 
C—H distance 1.09A 1.09A 
C—C or C—N distance | 1.54A 1.48A 
C=N or NSC distance | 1.16A 1.17A 
ZHCH 1093° 1094° 
ZHCC or ZHCN 1094° 1094° 
qi 5.3 x 107 5.3 10-* 
I,=T1s3 97.3 X 10-* 86.8 x 10-4 
I,Iols3 50,000 x 10-2 | 39,900 « 10-120 
Molecular weight 41.05 41.05 
Symmetry number 3 3 
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TABLE II. Vibration frequencies of methyl cyanide and 
methyl isocyanide in cm, 


EWELL AND J. F. 


BOURLAND 


TABLE III. Calculated thermodynamic properties of methyl 
cyanide in the ideal gaseous state. 








METHYL 
ISOCYANIDE 


3000° 

928! 
2161! 
1414! 
1456? 
1041? 

290? 


METHYL 
CYANIDE 


3000* 

9171 
2250! 
13701 
1417? 
1040? 

376? 





C—H stretching? 

C—C or C—N stretching! 
C=N or N=C stretching! 
HCH bending? 


HCC or HCN bending? 
CCN or CNC bending? 








of the rigid rotator-harmonic oscillator approxi- 
mation. Tables I and II give the assembled 
molecular data. 

The calculations were made using the well- 
known relations 


F°— F,° 
a In M+4R In T7+1/2R In Glels 


12 
~R In 6 +:257.40—R ¥ In (1-e-*), 
1 


H°—E,° 12 et 
———=4R+R)> 


’ 
1 e*—1 


— H®—E,° F°—Ep® 
. — T T ’ 





12 x2er 


C,=4R+R , 
1 (e7— 1)? 


where 
x=hei/kT =1.43270/T 
R=1.9869. 


Inserting the data from Table I, gives 


CH;CN : 
Fo— FE, F°— F° 
———— =2.53+4R In r-(——) ’ 
z T vib 


CH;NC: 


F°—E,° Fo—E,° 
——-=2.30+4R In T— (+) : 
T T vib 


The Einstein vibrational functions were evalu- 
ated by using a table calculated by Sherman and 
Ewell. The calculated values of the thermo- 
dynamic functions are given in Tables III and 


3 J. Sherman and R. B. Ewell, unpublished. 


1°C —(F°—Eo°)/T 


0 47.72 
25 48.57 
50 49.36 

100 50.82 
150 52.14 
200 53.39 
300 55.69 
400 57.90 
500 59.69 


(H° —Eo°)/T So Cp® 


9.42 57.14 11.93 
9.65 58.22 12.48 
9.89 59.25 13.12 
10.38 61.20 14.12 
10.90 63.04 15.29 
11.40 64.79 16.19 
12.39 68.08 17.93 
13.34 71.24 19.62 
14.28 73.97 21.04 























TABLE IV. Calculated thermodynamic properties of methyl 
tsocyanide in the ideal gaseous state. 








rc —(F°—Eo°)/T | (H° —Eo®)/T So Cp® 


0 47.87 9.83 57.70 
25 48.75 10.05 58.80 
50 49.68 10.28 59.96 

100 51.10 10.74 61.84 
150 52.47 11.22 63.69 
200 53.76 11.70 65.46 
300 56.09 12.66 68.75 
400 58.21 13.56 71.77 
500 60.21 14.47 74.68 























IV. The heat capacity data were fitted to three 
constant equations as follows: 


CH3CN: C,=4.77+29.1XK10°T 
+10.45 x 10-*7? 
(0.5 percent, 273-773°K), 


CH3NC: C,=6.14+24.2k10°7+6.16X 10-7" 
(0.5 percent, 273-773°K). 


It is unfortunate that the heat of combustion 
data for methyl cyanide and methyl isocyanide 
are so uncertain,‘ otherwise the equilibrium 
constants for some interesting and possibly tech- 
nically useful reactions could be calculated. 
Some calculations have been made with the 
combustion data tabulated by Bichowsky and 
Rossini.4 While the results are too unreliable to 
warrant giving them in detail, the general con- 
clusions are briefly as follows: (1) Methyl cyanide 
could be obtained in good yields by dehydrogena- 
tion of ethylamine above 350°C; (2) Ethylene 
and ammonia could combine to give methyl 
cyanide and hydrogen in good yields above 
500°C. 


4F. R. Bichowsky and F. D. Rossini, Thermochemistry 


of Chemical Substances (Reinhold Publishing Corporation, 
New York, 1936). 
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ANGSETH and his co-workers! have re- 

cently applied the results of essentially 
incomplete spectroscopic studies of liquid cyclo- 
hexane, symmetrical tetrachloroethane, and 
ethylene deuterobromide to a discussion of the 
intramolecular forces restricting internal rotation 
about the C—C bond. We believe that none of 
their structural conclusions is correct. Their 
discussion is based on their conclusion that in 
these molecules the opposed or eclipsed configu- 
rations are the stable ones. Insofar as liquid 
cyclohexane and symmetrical tetrachloroethane 
are concerned this conclusion is most probably 
incorrect since it directly contradicts the results 
of a great number of essentially more straight- 
forward studies of these and similar molecules in 
the gaseous and solid phases, and since it is 
probable that molecules in the liquid generally 
have structures closely related to those in the 
gas and the solid, at least when these are not 
too different. But whatever may be the structures 
of these molecules in the liquid, it is certainly 
true that the structures of the free molecules are 
more pertinent to a discussion of intramolecular 
forces. These structures indicate that there is a 
general tendency for the staggered configurations 
to be the stable ones, and Langseth’s general 
discussion!* must accordingly be without signifi- 
cance. Some of the existing evidence on cyclo- 
hexane, symmetrical tetrachloroethane, and re- 
lated molecules is briefly given below together 
with a preliminary report on an electron diffrac- 
tion investigation of the structure of gaseous 
symmetrical tetrachloroethane. We know of no 
direct evidence regarding the stable configuration 
of the ethyl bromide molecule. We believe, 
however, that in ethyl bromide and in ethane it 
is the staggered configuration which is the more 
stable, as it is in the hydrocarbons and substi- 
tuted hydrocarbons for which unambiguous 


‘a. A. Langseth and B. Bak, J. Chem. Phys. 8, 403 
(1940). b. A. Langseth and H. J. Bernstein, ibid. 8, 410 
(1940). c. A. Langseth, H. J. Bernstein and B. Bak, ibid. 
8, 415 (1940). d. A. Langseth, H. J. Bernstein and B. Bak, 
ibid. 8, 430 (1940). 


structural information exists. It seems likely 
therefore that Langseth’s assignments of the two 
isomeric components!*? of the C—Br stretching 
frequency of ethylene deuterobromide should be 
reversed. 

X-ray studies of crystalline cyclohexane,** 
B-benzene hexabromide and hexachloride,’” and 
1,4-di-iodocyclohexane** uniformly eliminate the 
possibility of the planar cyclohexane ring, pro- 
posed by Langseth and Bak, as does the electron 
diffraction investigation of gaseous cyclohexane 
by Pauling and Brockway.‘ Beach and his 
co-workers’ have found staggered configurations 
for ethylene chloride, chlorobromide, and _ bro- 
mide, and 1,1,2-trichloroethane as we have for 
the 2,3-dibromobutanes and propylene bromide. 
Klug reported crystalline ethylene iodide’ to 
have a irans and thus staggered configuration. 

Our electron diffraction photographs® of gase- 
ous symmetrical tetrachloroethane were obtained 


TABLE I. 








So Xx. N. So 


1.61 11.85 
2.63 12.74 
3.68 13.64 
4.74 14.38 
5.54 15.54 
6.14 16.72 
6.77 17.97 
7.68 18.70 
8.47 19.54 
9.43 20.55 
10.54 21.85 





—13 











2 M. de Hemptinne and C. Velghe, Physica 5, 958 (1938). 
3a. O. Hassel and H. Kringstad, Tidskr. Kemi. Berg- 
vesen, 10, 128 (1930). b. R. G. Dickinson and C. Bilicke, 


J. Am. Chem. Soc. 50, 764 (1928). c. E. Holméy and 
O. Hassel, Zeits. f. physik. Chemie B16, 234 (1932). 

4L. Pauling and L. O. Brockway, J. Am. Chem. Soc. 
59, 1223 (1937). 

5a. J. Y. Beach and K. J. Palmer, J. Chem. Phys. 6, 
639 (1938). b. J. Y. Beach and A, Turkevich, J. Am. 
Chem. Soc. 61, 303 (1939). c. A. Turkevich and J. Y. 
Beach, tbid. 61, 3127 (1939). 

6a. D. P. Stevenson and V. Schomaker, J. Am. Chem. 
Soc. 61, 3173 (1939). b. V. Schomaker and D. P. Stevenson, 
unpublished observations. 

7H. P. Klug, Zeits. f. Krist. 90, 495 (1935). 

8L. O. Brockway, Rev. Mod. Phys. 8, 231 (1936). 
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Fic. 1. The radial distribution function for symmetrical 
tetrachloroethane. The solid vertical lines give the 
interatomic distances to be expected for a molecule with 
bond distances and angles as in ethyl! chloride and methyl- 
ene chloride. The lengths of the lines indicate the relative 
importance of the distances in the scattering formula for 
the more symmetrical (C2,) model. For the less sym- 
metrical staggered model (C2) the 3.25A line would become 
three times as long as the 4.31A line. The dotted lines 
correspond to the orientation-dependent distances in the 
symmetric (Coy) eclipsed model; for the C2 eclipsed model 
the longer of these distances becomes three times as 
important as the shorter one. The arrows relate the 
corresponding orientation-dependent distances of the two 
types of models, 


from a sample (b.p. 142.8+0.1°) which was part 
of the middle fraction obtained on distilling a 
large quantity of Eastman symmetrical tetra- 
chloroethane through a 25-cm helix packed 
column with controllable reflux ratio. The photo- 
graphs were taken with the high temperature 
nozzle. The temperature of the boiler was 
varied from 70° to 120°. The wave-length of the 
electrons was 0.0615A (ao for gold, 4.070A) and 
the camera distance was 10.85 cm. 

The photographs showed eleven rings whose 
measured s values (s= (47/) sin (#/2)) are given 
in Table I along with the coefficients C, of the 
radial distribution function.!° Characteristic 
features of the diffraction pattern are the rela- 
tively weak third and fifth maxima, the broad 
appearance of the eighth and the shelf-like 
nature of the tenth maximum. 

The radial distribution function shown in 
Fig. 1 has maxima at 1.74, 2.86, 3.31 and 4.32A. 
The first of these (1.74A) corresponds to the 

*L. O. Brockway and K. J. Palmer, J. Am. Chem. Soc. 
59, 2181 (1937). 

104. L. Pauling and L. O. Brockway, J. Am. Chem. Soc. 


57, 2684 (1935). b. Verner Schomaker, American Chemical 
Society Meeting, Baltimore, Maryland, April, 1939. 


VERNER SHOMAKER AND D. P. 








STEVENSON 


carbon-chlorine bond distance (1.76A), and the 
second (2.86A) to an unresolved average of the 
nonbonded carbon-chlorine distance (2.74A in 
ethyl chloride"*) and the distance between 
chlorine atoms bound to the same carbon atom 
(2.92A in methylene chloride”). The positions 
of the other two maxima at 3.31A and 4.32A 
correspond to the chlorine-chlorine distances 
(3.25A and 4.32A) which would be expected if 
the bond distances and bond angles are the 
same as the corresponding distances and angles 
in ethyl chloride and methylene chloride and if 
the molecule has one (or both) of the ‘‘staggered”’ 
relative orientation of the two —CHCLz2 groups 
about the C—C bond. Moreover, it may be seen 
from Fig. 1 that these peaks are in complete 
disagreement with the distances to be anticipated 
for the ‘‘opposed’”’ configuration proposed by 
Langseth and Bernstein.'» The peaks of the 
radial distribution function which depend on 
the relative orientation of the two ends of the 
molecule are quite narrow. It therefore seems 
likely that the molecules are fairly rigid. In 
particular, the sharpness of the 3.31A peak 
indicates that the amplitude of libration is less 
than that (15—20°) which has been estimated for 
ethylene chloride, etc.s® The approximate 3 : 2 
ratio of the heights of the 3.31A and 4.32A 
peaks of the radial distribution function suggest 
that the gas at the temperatures of our experi- 
ments consist of a roughly equimolal mixture of 
the two staggered isomers. We believe that the 
complete treatment of our data by the correlation 
method will provide a more precise estimate of 
the relative abundance of the two isomers as 
well as the answer to the rather large discrepancy 
between the position of the 3.31A peak of the 
radial distribution function and the value calcu- 
lated from our model. Our tentative conclusions 
regarding the relative abundance and rigidity of 
the isomers are in good agreement with those of 
Smyth and McAlpine” from dipole moment 
measurements on the gas. 

We are indebted to Professor Linus Pauling 
for helpful counsel about the matters discussed 
in this paper. 

4a. J. Y. Beach and D. P. Stevenson, J. Am.’ Chem. 
Soc. 61, 2643 (1939). b. L. E. Sutton and L. O. Brockway, 
thid. 57, 479 (1935). 


2C, P. Smyth and K. B. McAlpine, J. Am. Chem. Soc. 
57, 979 (1935). 








id the 
of the 
4A in 
tween 
atom 
sitions 
4.32A 
tances 
ted if 
e the 
angles 
and if 
ered”’ 
TOUps 
e seen 
plete 
pated 
d by 
f the 
id on 
f the 
seems 


1. In 
peak 
s less 
d for 
$:2 
.32A 
gest 
(peri- 
ire of 
t the 
ation 
te of 
rS as 
ancy 
' the 
alcu- 
sions 
-y of 
se of 
nent 


ling 
ssed 


hem. 


way’, 


Soc. 


AUGUST, 1940° 


JOURNAL OF CHEMICAL PHYSICS 


VOLUME 8 


Thermal Diffusion Separation of Radioactive and Ordinary Hydrogen Isotopes 


G. T. SEABorG, A. C. WAHL AND J. W. KENNEDY 
Department of Chemistry, University of California, Berkeley, California 


(Received June 11, 1940) 


T seems reasonable that in the separation of 

gaseous isotopes by the thermal diffusion 
method of Clusius and Dickel'? the separation 
factor* should be independent of the relative 
concentrations of the two components. The 
simplified equation of Enskog* indicates that the 
separation factor for the elementary thermal 
diffusion process is independent of concentration. 
However, it is not impossible that a variation 
in S.F. with mole fraction could arise from the 
process producing multiplication of the unit 
separation by convective action; also, the results 
on the separation of the chlorine isotopes re- 
ported by Clusius and Dickel* indicate a varia- 
tion in S.F. amounting to a factor of about two 
or three between the values at high mole fractions 
(70 percent to 40 percent Cl**) and at low mole 
fractions (12 percent to 0.6 percent Cl*®).7 

To test the behavior of a Clusius and Dickel 
column at extremely small mole fractions we 
have employed a radioactive isotopic species, 
namely the long-lived radioactive H* of Alvarez 
and Cornog.4 A column was constructed 7.5 
meters in length and with an inside diameter of 
1.0 cm. A water-cooling jacket served to keep 
this tube at room temperature; provision was 
made for electrically heating an axial, molybde- 
num wire (0.40 mm diameter). A 5-liter gas 
reservoir was connected to the top of the 
column, and samples were withdrawn from the 
bottom. ; 

In a set of preliminary experiments we studied 


'K. Clusius and G. Dickel, Naturwiss. 26, 546 (1938); 
27, 148 and 487 (1939). 

°K. Clusius and G. Dickel, Zeits. f. physik. Chemie 
B44, 397 and 451 (1939). 

*For this discussion we use the separation factor 
(S.F.) =(Ni/N2)pottom/(Ni/Ne2)top, where N; and N2 are 
the mole fractions of the isotopes. 

3 See reference 2, page 427. 

} An algebraic error in the derivation of an equation in 
their paper (reference 2, Eq. (26), p. 430) may have led 
them to exaggerate the magnitude of this variation (Fig. 8, 
p. 471). In addition the experimental values were neces- 
sarily determined in a region where the effect of errors in 
measurement are magnified. 

(1939) W. Alvarez and R. Cornog, Phys. Rev. 56, 613 
9). 


the separation of deuterium-hydrogen mixtures. 
With the wire at 250°C and the gas at 1 atmos- 
phere pressure, equilibrium was reached after 
3 days and thermal conductivity analyses of 
samples removed showed 18 percent deuterium 
in the reservoir and 87 percent deuterium at the 
bottom of the column. These figures give S.F. 
=30. (With the wire at higher temperatures 
S.F. could not be accurately measured because 
the deuterium content of the samples withdrawn 
from the bottom was too large.) 

A mixture of deuterium and radioactive H?* 
(mole fraction ~10-") was introduced into the 
same column under identical conditions of 
temperature and pressure as above. The samples 
withdrawn were analyzed for H* content by 
introduction into a Geiger-Miiller counter. At 
equilibrium we found S.F.=5.0. When account 
is taken of the mass dependence of S.F. according 
to the equations* of Enskog or Chapman it is 
seen that the effectiveness of the column for 
separating hydrogen isotopes is substantially 
the same at moderate or extremely small mole 
fractions. 

We have performed exploratory experiments 
concerning the dependence of S.F. on the 
temperature of the hot wire, using deuterium- 
radioactive H* mixtures for convenience in 
analysis. S.F. increased with increasing tempera- 
ture, reaching the value 9 at the highest temper- 
ature available, 800°C. (To prevent vibration 
of the wire at this temperature spacers in the 
column were required even though direct current 
was used for heating.) It should be possible to 
increase further this separation by finding the 
optimum values of temperature, pressure and 
tube radius. Another way to increase S.F. is to 
use several such columns in series.” It is easily 
seen that if the S.F. is independent of the mole 
fraction, as is demonstrated approximately by 
the present experiments, then the S.F. for several 
columns in series will be given by the product of 
the individual S.F.’s. We have actually connected 
in series with the column described above a 
second column, which we have previously used 
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in separating® the chlorine isotopes, and obtained 
the expected S.F. although a longer time was 
required to reach equilibrium (with our coupling 
arrangement about 10 times as long as for a 
single column). 

It would be interesting to know how well the 
Clusius and Dickel column would separate 
deuterium-hydrogen mixtures at the higher 
temperature (800°C). If the temperature de- 
pendence of this S.F. is the same as that of the 
deuterium-radioactive H* mixture, at 800° we 
would expect S.F.=54. However, calculations 
based on the Enskog equation suggest a much 
higher value. (The same calculations give the 


6 toy Kennedy and G. T. Seaborg, Phys. Rev. 57, 843 


(1940). 


WAHL AND J. W. 


KENNEDY 


temperature dependence of the deuterium- 
radioactive H® separation within a factor of two 
for the region investigated.) A value of 54 for 
S.F., which we for these reasons consider a lower 
limit, would imply that, starting from a large 
reservoir of the natural isotopic mixture, one 
column of the type described would produce 
1.1 percent deuterium, two such columns would 
give 37 percent deuterium, and three columns 
could be used to deliver 97 percent deuterium. 

We are indebted to Professor E. O. Lawrence 
and other members of the Radiation Laboratory, 
in particular Dr. R. Cornog, for the sample of 
radioactive hydrogen. Also it is a pleasure to 
thank Dr. Melvin Calvin, who made available 
to us his thermal conductivity apparatus. 
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LETTERS TO THE EDITOR 


This section will accept reports of new work, provided these 
are terse and contain few figures, and especially few halftone 
cuts. The Editorial Board will not hold itself responsible for 
opinions expressed by the correspondents. Contributions to 
this section should not exceed 600 words in length and must 


Silicate Glasses 


Volumes, refractions and dispersions of well-annealed 
silicate glasses are simple additive functions of the compo- 
sitions, when these properties are calculated for the 
quantity of glass which contains one gram atom of oxygen 
and the compositions are expressed as numbers of atoms 
(Nw) of the “metallic” elements per atom of oxygen. 
Specifically, 

Vo=k+bgitcsiNsit2'cuNm, 
(m,—1) Vo= Rox = Zam,r.Nm, 
and = (m,, — y,) Vo= Do, a4, = 2 (Gm, r, — 2, dq) Nau. 


k is a small constant, depending apparently on the 
annealing technique. The constants bs; and cgi have 
different values for different ranges of Ngi, as shown by 
the changes of slope of the line in the accompanying figure 
for NazO—SiO, glasses. The same constants are applicable 
over the same ranges of Ng; values, regardless of the kinds 
and amounts of other metallic elements present, insofar 
as can be determined from the data at present available. 
The different slopes indicate a d!fference in the closeness 


Vo,% ° Yo" * - 670%, 





Na,0- SiO, 


* MOREY AND 
© FINN AND CO-WORK! 
@ TURNER AND CO-WORKERS 


Fic. 1. 
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reach the office of the Managing Editor not later than the 15th 
of the month preceding that of the issue in which the letter 
is to appear. No proof will be sent to the authors. The usual 
publication charge ($3.00 per page) will not be made and no 
reprints will be furnished free. 


of the packing of the atoms in the silicon-oxygen framework 
within different Ng; ranges. 

The break at Ng,;=0.435 may indicate the limit, below 
which all the oxygen atoms are approximately in a close- 
packed arrangement and above which some of them are 
more loosely packed. Ng;=0.40, where another break 
occurs, is the composition required for (SizOs~~), sheets 
or other structures in which 3 of the 4 oxygens around 
each Si are bonded also to another Si. The break at about 
0.345 is not very well defined by present data. Actually 
there may be a gradual curvature in this region or there 
may be two or more breaks, e.g. at 0.333, corresponding 
to (SiO;~~), chains or rings (in which 2 of the 4 oxygens 
surrounding each Si form bridges to other silicons), and at 
0.364, corresponding to (Si,O;:~*), double chains. 

Satisfactory data for glasses having Ngi less than about 
0.27 are lacking. One may predict a definite break at 0.25, 
however, since below that concentration some of the 
oxygens must be bonded to no silicons. 

Details of this work will be published in the Journal of 
the Optical Soctety of America. 


Maurice L. HuGGINs 
Kodak Research Laboratories, 
Rochester, New York, 
July 1, 1940. 





Elasticity of Crystals 


Recent theoretical discussions! of the elasticity of 
crystals, particularly of the rigidity, in the neighborhood 
of the melting point have brought out the scantiness of 
relevant experimental material. Some preliminary meas- 
urements of the rigidity of polycrystalline aluminum, to 
650°C, may consequently be of interest. A dynamical 
method? gave the resonant frequency of torsional vibrations 
of a right circular cylinder, approximately 20 cm long. 
The frequencies (of the lowest mode) varied from about 
8000 cycles per second at room temperature, to about 
5100 at 650°. The melting point of this aluminum (com- 
mercial 2S grade) as determined with the same Chromel- 
Alumel thermocouple which was used for the other 
temperature measurements, was 663°, using the calibration 
tables given by the manufacturer. In Table I are given 
the ratios of frequencies at even temperatures to the 
frequency at 0°C and the corresponding ratios of rigidity, 
computed from the relations, V=2FL and G= V%p, where 
V is the velocity of propagation of torsional waves, F the 
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TABLE I. Relative rigidity of aluminum. 








F /Fo 


1.000 
0.973 
.942 
.907 
-865 
816 
.756 
.750 .553 
.743 .542 
-731 .524 
-712 498 
648 412 
melting point 


G /Go 





1.000 
0.945 
-883 
-816 
-740 
-657 
.562 








fundamental frequency of vibration of a cylindrical rod of 
length L, p the density, and G the rigidity. Corrections for 
thermal expansion have been applied, according to the 
data of Hidnert. 

These measurements show a considerably smaller falling 
off of rigidity than do those of Koch and Dannecker,? 
also carried out on polycrystalline aluminum, to 580°; the 
method of the torsional pendulum used by these writers 
probably involved strains beyond the elastic limit at the 
higher temperatures. 

It is evidently not possible to say from these results 
whether the rigidity vanishes at the melting point or 
retains a finite small value. The interpretation is somewhat 
complicated by the fact that the measurements deal with 
an aggregate of crystals; it seems unlikely, however, that 
an aggregate would remain stable for shearing strains after 
the crystals of which it was composed had become unstable 
for shear about any single crystal direction. 

The shape of the curve giving the variation of rigidity 
with temperature is perceptibly different from the shape 
of the sample curves calculated by Born: The observed 
variation is smaller at low temperatures, and the decrease 
in the neighborhood of the melting point more precipitate. 
Born’s curves were calculated, however, for the case of 
central forces between atoms. 

FRANCIS BIRCH 


DENNISON BANCROFT 
Harvard University, 
Cambridge, Massachusetts, 
June 24, 1940. 


1 Max Born, J. Chem. Phys. 7, 591 (1939); L. Brillouin, Phys. Rev. 
54, 916 (1938). 

2? Francis Birch, J. App. Phys. 8, 129 (1937). 

3 Koch and Dannecker, Ann. d. Physik 47, 197 (1915). 





Lattice Summations for Hexagonal Close-Packed Crystals* 


The potential energy of interaction between the atoms 
or ions of a crystal may be expressed by a set of terms 
which have the exponential or inverse power form of 
attraction or repulsion. The work required to remove a 
single particle from the crystal is determined by summing 
up the interaction between that particle and all the other 
particles of the crystal. This summation for the inverse 
power terms has been made by Lennard-Jones for the 
cubic crystal.1 The summations necessary for the determi- 
nation of the elastic constants of a cubic crystal have been 
made by Herzfeld and Goeppert-Mayer? and by May.® 
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For the hexagonal crystal, the method used by Lennard- 
Jones could not easily be adopted, so we have used a more 
tedious and less elegant method. The summations were 
made directly to ten times the distance between nearest 
neighbors omitting half the contribution of the outside 
shell. Then an integration was made from this outside 
shell to infinity. The crystal considered was the hexagonal 
close-packed crystal with ideal axial ratio (1.633). Sum- 
mations were made over x—y planes for the inverse sixth, 
eighth, tenth, and twelfth power of the distance between 
particles. The threefold axis of the crystal was chosen as 
the z axis. From symmetry considerations and the fact 
that r?=x?+y?+2?, the following relations between the 
sums entering the elastic constants can be established: 
(1) Lxtr— 2-2 = Ly"? = h( Lr" — Ser), 
(2) Dxtr-2 4 = Ly'r-*-4§ = 3 Ixy n- 4 
(3) Dxty2r-2-4 = (Lr — 2222-24 Tatr24), 
(4) © Dxtstr—a4 x Dytety—-4 = (D4 "9-8 — Laty—n-4), 
(5) All sums containing odd powers of any one coordinate 
are zero. Consequently, only the lattice sums <r’, 
D22r-"-2, and LYz'r-"-4 have to be evaluated in order to 
obtain the contributions of the interaction potential r~" 


to all the elastic constants. 
The summations were made for »=6, and n=8. 


TABLE I. Latlice sums for hexagonal close-packed crystals. 








n=8 

12.8028 
4.270 
4.266 
2.804 
2.650 
0.883 
0.733 

12.1323 


n=6 


=r-" 14.4549 
4.821 
4.817 
3.114 
2.973 
0.991 
0.854 

12.3124 


Dzr-n-2 
Dx2r-n-2 
DVziy-n-4 
Dxir-n-4 
Dx2y2r-n-4 
Dx222x-n- 
Sr-n-4 














Since the lattice sums obtained have almost the same 
value as for the face-centered cubic crystal, a direct 
comparison of the two lattice types was made as a numer- 
ical check. If one of the threefold axes of the cubic crystal 
is chosen as the z axis then the first and second x—y planes 
of the two types of crystal can be made to coincide. The 
next three planes differ in the two types and the following 
three differ only by a rotation about the z axis, which does 
not affect the summation. A summation was made of the 
inverse sixth and eighth powers over the first three planes 
on which they differed. The difference beyond these planes 
is inappreciable. The difference between the two types 
was found to be exactly the same as the difference taken 
over the entire crystal as found by us and by Lennard- 
Jones (AZr~*=0.0010 and AXr~*=0.0009). 


BROTHER GABRIEL KANE 
Manhattan College, 
MARIA GOEPPERT- MAYER 
Columbia University, 
New York, New York, 
July 11, 1940. 


* Publication assisted by the Ernest Kempton Adams Fund to! 
Physical Research of Columbia University. 
Lennard-Jones and A. E. Ingham, Proc. Roy. Soc. A107, 63 


.F. Herzfeld and M. Goeppert- Mayer, Phys. Rev. 46, 995 (1934)- 
3A. May, Phys. Rev. 52, 339 (1937). 
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An X-Ray Investigation of the Phase Transitions of 
Sodium Palmitate 


In the course of some phase transition studies carried 
out with the collaboration of Lever Brothers Company, 
a temperature-controlled x-ray powder camera! was 
employed as a means of detecting the transitions of 
sodium palmitate. As a result of this investigation evidence 
is provided for the existence of seven phases, with their 
corresponding six transitions, between room temperature 
and the condition of isotropic liquid. 

The sodium palmitate used in this experiment was 
prepared by neutralizing an alcoholic solution of Eastman 
Kodak Company palmitic acid with alcoholic sodium 
hydroxide. Alcohol and residual water were removed by 
evacuation at an elevated temperature. A hot alcohol 
solution of the sodium palmitate was exactly neutral to 
phenolphthalein. The anhydrous sodium palmitate was in 
the form of small fibrous chunks ranging in size from 
extremely fine powder to pieces approximately 1-2 mm 
on the long dimension. The material was not unlike 
satin-spar in appearance. 

A small amount of the crystalline solid was loaded and 
sealed into a fine glass capillary tube. All of the x-ray 
photographs were taken under identical conditions of 
exposure and development. The exposure time was two 
hours, in each case being preceded by 15 minutes of 
heating under the temperature conditions at which the 
exposure was to be made. These photographs indicated 
transitions at 67°, 117°, 138°, 205°, 257°, and 290°C. 

The transitions in sodium palmitate have already been 
investigated by R. D. Vold and M. J. Vold,? who employed 
a dilatometric method. They reported five successive 
phase changes. The 67° transition found here by x-ray 
methods was not reported by the Volds. The revised list 


TABLE I. 








TRANSITION TEMPERATURE* 





X-Ray METHOD DILATOMETER 
(CHESLEY) (VOLD) 





Curd fiber | 
Curd fiber II 
Subwaxy 
Waxy 
Subneat 


Neat 


Isotropic liquid | 








ian Since x-ray photographs were taken at 5-degree intervals, the 
sbove temperatures are correct to +2.5°C 
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of successive phases and their corresponding transitions, 
using the Vold nomenclature, is shown in Table I. 

It may be recalled that no single method provides an 
infallible criterion of a transition. The dilatometric 
method fails if the volume change at a particular transition 
is extremely small. The data listed by the Volds proves 
that the volume change from the phase Curd fiber I to the 
phase Curd fiber II is indeed very small. Nevertheless, 
the data? plotted in their Fig. 1 do show a small volume 
discontinuity in the region near 70°C, This they either 
overlooked or failed to recognize as being characteristic of 
a transition. The x-ray method proves quite sensitive in 
detecting the 67° transition. 

FRANK G. CHESLEY 

Mineralogical Laboratory, 

Massachusetts Institute of Technology, 


Cambridge, Massachusetts, 
June 26, 1940. 


1 A full discussion of this instrument will appear elsewhere. 
2 R. D. Vold and M. J. Vold, J. Am. Chem. Soc. 64, 808 (1939). 





Raman Spectrum of Aluminum Bromide 


The Raman spectrum of a very pure sample of aluminum 
bromide, kindly supplied by Dr. H. Pines of the Universal 
Oil Products Company, has been photographed. The 
material was kept in the liquid state at 100° by the heat 
from the mercury arc light sources. The results are as 
follows: 67.0(5b)(e, +k), 79.3(5)(e, +k), 112.8(5)(te, +h), 
140.3(5)(-te, +), 185(0)(e), 208.2(10)(+e, i, +h), 
223(1)(+e), 407?(Ovd)(e), 488(1vd)(e). It seems likely that 
407 and 488 are not fundamentals. 185 may be the combi- 
nation of 113 and 67, and 223 may be the combination of 
140 and 79. 

Two plausible structures for Al,Bre are the halogen 
bridge model with symmetry D2, adopted by Palmer and 
Elliott! and the ethane-like model with symmetry Ds, 
Ds, or Dsa. The maximum number of frequencies active 
in the Raman effect permitted by the selection rules is 9 
for the former model, while this number is 10, 9, or 6 for 
D3, Ds, or Dsa, respectively. According to Stitt? the 
number of distinct Raman-active fundamentals for Ds 
with free rotation is 6. Since the symmetry Dy (plane of 
reflection perpendicular to the threefold axis) is not very 
probable for Al2Bre, 6 Raman-active fundamentals might 
be expected for the ethane-like structure and 9 for the 
halogen bridge structure. While a clean-cut distinction 
between these alternatives cannot yet be made, it is 
probably significant that (1) we have observed 5 strong 
Raman lines both with Al,Bre and with Alo(CHs)s for which 
a bridge structure is presumably impossible; (2) the Raman 
spectrum of C2Cl¢ consists of 5 lines;? and (3) six funda- 
mentals of SizCle are observed in its Raman spectrum’. 

The study of this point will be continued by observation 
of the Raman and infra-red absorption spectra of some 
halogen-methyl compounds of aluminum. 

E. J. ROSENBAUM 

Department of Chemistry, 

University of Chicago, 


Chicago, Illinois, 
June 20, 1940. 


1 Palmer and Elliott, J. Am. Chem. Soc. 60, 1852 (1938). 
2 Stitt, J. Chem. Phys. 5, 90 (1937). 
3 Heidenreich, Zeits. f. Physik 97, 277 (1935). 
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The Determination of the Proton-Attracting Properties 
of Liquids* 

In a recent article! Gordy and Stanford have shown 
that a plot of the shift in wave-length of the OD vibrational 
band of the CH;OD molecule due to hydrogen bonding 
with a variety of amines of known basicity is a linear 
function of the logarithm of the ionization constant K; 
of the base. By a bold extrapolation they have then 
estimated the base strengths of a number of weakly basic 
oxygen compounds, but have apparently overlooked the 
fact that a value for the basicity of one of these, aceto- 
phenone, has previously been obtained by an entirely 
different method.? This depends upon the determination 
of the range of sulfuric acid-water compositions required 
for the change in ultraviolet absorption which accompanies 
the conversion of the base to its conjugate acid. Gordy 
and Stanford’s value of K; for acetophenone is 1X 10~*, 
that obtained from the data of Flexser, Hammett and 
Dingwall is 0.9X10-*°. Further, somewhat less direct 
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evidence of agreement of the two methods is also available. 
The Gordy and Stanford value for methyl and benzyl 
benzoates is 5X10-*; the Flexser, Hammett and Dingwall 
value for benzoic acid is 6X 10-*. The Gordy and Stanford 
value for nitrobenzene is 4X 10-*5; one may estimate from 
the incomplete ionization of the substance in nearly 
anhydrous sulfuric acid* that the value of K; is in the range 
from 10~* to 10~*5, 

The agreement between these two entirely different 
approaches to the determination of the strengths of very 
weak bases strikingly confirms the validity of both. 

Louis P. HAMMETT 

Department of Chemistry, 

Columbia University, 


New York, New York, 
July 13, 1940. 


* Publication assisted by the Ernest Kempton Adams Fund for 
Physical Research of Columbia University. 

1 Gordy and Stanford, J. Chem. Phys. 8, 170 (1940). 

— Hammett and Dingwall, J. Am. Chem. Soc. 57, 2103 
(1935). 

3 Treffers and Hammett, J. Am. Chem. Soc. 59, 1708 (1937). 





